Chaotic Pendulum. Experiment CP. University of Florida Department of Physics PHY4803L Advanced Physics Laboratory. Objective. References.
|
|
- Wesley Lane
- 5 years ago
- Views:
Transcription
1 Chaotic Pendulum Experiment CP University of Florida Department of Physics PHY4803L Advanced Physics Laboratory Objective Both regular and chaotic motion for a springpendulum system will be measured and analyzed. You will learn in your investigations how to recognize systems for which chaotic behavior can be expected, how to display data arising from these systems, how to quantify features of the dynamical variables involved, and how to perform computer simulations and comparisons with experimental data. References Henry D. I. Abarbanel, The analysis of observed chaotic data in physical systems, Reviews of Modern Physics, 65, (1993). Gregory L. Baker, Jerry P. Gollub, Chaotic Dynamics: An Introduction, 2nd ed., Cambridge University Press, New York, Francis C. Moon, Chaotic and Fractal Dynamics, John Wiley and Sons, New York, Introduction The harmonic oscillator is a paradigm of predictability. The resonance response of a harmonic oscillator can be exploited to make clocks accurate to 1 part in Harmonic oscillations are normally observed in systems operating near equilibrium and governed by linear or Hooke s law restoring forces. Nonetheless, for small enough excursions from equilibrium, systems governed by nonlinear restoring forces can also display harmonic oscillations. A good example is the pendulum. When its motion stays near equilibrium (at θ = 0), the nonlinear restoring torque (proportional to sin θ) behaves like a linear restoring torque (proportional to θ) and the pendulum executes simple harmonic motion. However, when driven strongly, the excursions from equilibrium can grow into regions where the nonlinearity becomes important. The predictable pattern of repeating oscillations might give way to chaos nonrepeating motion characterized by a particular kind of unpredictability. Chaotic systems are completely deterministic. For example, our apparatus is governed by equations easily derived from F = ma. The unpredictability arises not from any inherent randomness but rather from an extreme sensitivity to initial conditions. Harmonic oscillators ultimately settle on some well defined state of motion. Momentarily perturbed from this steady state motion, they will return to it according to a predictable transient time constant. In contrast, perturb a nonlinear oscillator and the difference between the CP 1
2 CP 2 Advanced Physics Laboratory perturbed and unperturbed states of motion can grow exponentially. You will investigate just such a nonlinear oscillator. Illustrated in Fig. 1, this springpendulum apparatus can be configured to display both chaotic and non-chaotic behavior. The non-chaotic data sets will be graphed and analyzed using traditional techniques. The chaotic data sets will be graphed in Poincaré sections and quantified by fractal dimensions and Lyapunov exponents. SXOOH\ D[OH GLVN θ SHQGXOXP PDVV VWULQJ VSULQJV Exercise 1 The applicability of the superposition principle is the distinguishing property between linear and nonlinear systems. Show that solutions to the pendulum equation, d 2 θ/dt 2 = (g/l) sin θ, obey the superposition principle only if one can assume sin θ = θ throughout the motion. That is, if θ 1 (t) and θ 2 (t) are solutions, show that θ 1 + θ 2 is a solution when sin θ = θ, but not otherwise. Apparatus The apparatus consists of a physical pendulum coupled to a drive motor using a pair of springs, a pulley, and string. The physical pendulum consists of an aluminum disk, a removable pendulum mass, a pulley, and an axle. The axle and pulley are those of the Pasco Rotary Motion Sensor which has a rotary encoder for measuring the pendulum angle. A computer equipped with a National Instruments PCI 6601 timer/counter board controls the drive motor and reads the rotary encoder. The computer is also used to save, display, and analyze the resulting data. The drive system consists of a stepper motor, shaft, and controller. A uniform pulse train is sent by the computer to the controller which sequentially energizes the various motor windings to effect a single angular step of the motor shaft for each pulse. The stepper motor VWHSSHU PRWRU Q φ $ VWULQJ VKDIW WXQLQJ SHJ Figure 1: The chaotic pendulum apparatus. Not shown: electronics, supports, and the rotary encoder housing the axle. requires 200 pulses to complete one revolution. As the stepper motor rotates, the driven end of the spring oscillates, keeping the pendulum in motion. For measuring the pendulum angle, the rotary encoder emits logic pulses (1440/rev) on two separate phonojacks for clockwise and counterclockwise rotations. The computer is capable of direct up-down counting with these signals. The count, directly proportional to the rotation angle, is read with every stepper motor pulse and saved into an array. Because
3 Chaotic Pendulum CP 3 the same pulses initiate a reading of the angular count and step the motor, the pendulum angle and the drive angle are simultaneously determined. The raw counter readings, spaced 1/200th of the drive period apart in time, are scaled and numerically smoothed and differentiated using Savitsky-Golay filtering. 1 The filters are equivalent to a separate least squares fit to a polynomial for each data point with the fitting region symmetrically surrounding the point. They are used to determine the angle θ, the angular velocity ω = dθ/dt, and the angular acceleration α = d 2 θ/dt 2. Equation of Motion As shown in Fig. 1, the pendulum angle θ is defined from the vertical with θ = 0 for the straight-up, or inverted, position. Clockwise angles are taken positive. Also illustrated on Fig. 1, the drive phase φ is taken as the angle of the motor shaft from vertical. The uniform pulse train sent to the motor has an adjustable rate and will cause the shaft to rotate at a constant angular velocity Ω so that φ = Ωt (1) where the shaft is assumed to be straight up at t = 0. Then, the displacement d of the driven end of the spring can be taken to be 2 d = A cos Ωt (2) 1 Abraham Savitsky and Marcel J. E. Golay, Smoothing and differentiation of data by simplified least squares procedures, Analytical Chemistry, 36 pp , With the distance from the center of the drive motor to the guide hole given by L, the driven end of the spring moves according to d = A cos φ A2 4L (1 cos 2φ). The offset term, A 2 /4L, is typically no larger than 2 mm and the term proportional to cos 2φ causes a small amplitude driving oscillation at an angular frequency of 2Ω which will be ignored in the analysis. Figure 2: Near scale drawing of the potential V (θ) for our system. The disk, pulley and axle have a moment of inertia I 0. They are well balanced and cause virtually no torque in a gravitational field. The pendulum mass m placed a distance l from the rotation axis produces a gravitational torque mgl sin θ. It also increases the moment of inertia to I = I 0 + I m + ml 2, where I m is its principal moment of inertia. The string between the two springs rides over the pulley without slipping, and the relatively weak springs stay elongated at all times. The springs/strings/pulley arrangement is adjusted as described later so that the springs have equal stretch (provide equal tension) when θ = 0 and d = 0. In this inverted configuration, the pendulum would have no torque acting on it, but it would be in an unstable equilibrium, ready to fall to the left or right. Exercise 2 (a) Show that the torque arising from the two springs and from gravity acting on the pendulum mass is given by τ c = 2kr 2 θ + krd + mgl sin θ (3) where k is the force constant for each spring and r is the pulley radius. (b) Show that for d = 0, (i.e., A = 0, which means the pendulum is undriven) the potential energy is given by V = kr 2 θ 2 + mgl cos θ (4)
4 CP 4 Advanced Physics Laboratory (c) Show that for mgl 2kr 2, V (θ) has one minimum at θ = 0; and that for mgl > 2kr 2, V (θ) has a maximum at θ = 0 and two equally deep minima at ±θ e where sin θ e θ e = 2kr2 mgl (5) As shown in Fig. 2, our system potential has a double well shape. Dissipation The torque τ c (Eq. 3) arises from the conservative forces of gravity and Hooke s law. There are also nonconservative torques due to eddy current damping and axle friction through which energy is continually drained or dissipated from the system. If this energy is not continually replenished, for example, by the stepper motor drive system, the motion slows and ultimately stops. A damping torque proportional to and opposite the angular velocity is created by induced eddy-currents when a neodymium magnet is brought near the rotating aluminum disk. This torque can be modeled as b ω. Friction, arising largely in the axle bearings, causes an additional damping torque of constant magnitude b that is always opposite the angular velocity. This axle friction torque can be modeled b sgn ω, where sgn ω = ω/ ω. Thus, the net nonconservative torque is given by τ f = b ω b sgn ω (6) This form for the damping torque will be verified by fits of data taken from the apparatus. It will also be used in simulations. The (probably small) effect of static friction as the pendulum momentarily stops each time ω goes through zero has not been explored. Occasionally, the entire axle friction term will be dropped in this writeup for readability or whenever analytic solutions are desired. (Including axle friction makes the equation of motion nonlinear and the solutions are not expressible in closed form.) Exercise 3 Show that the equation of motion for the apparatus can be written as an autonomous 3 set of three first order differential equations. θ = ω (7) ω = Γω Γ sgn ω κθ (8) + µ sin θ + ɛ cos φ φ = Ω (9) where a dot over a variable represents its time derivative and Γ = b/i (10) Linear Dynamics Γ = b /I (11) κ = 2kr 2 /I (12) µ = mgl/i (13) ɛ = kra/i (14) Despite being nonlinear, the system can still behave like a harmonic oscillator as long as the total energy remains small compared to the height of the potential barrier at θ = 0. Exercise 4 (a) Show that a Taylor expansion 4 of V (Eq. 4) about either minimum θ 0 = ±θ e gives V = kr 2 θ 2 e + mgl cos θ e (15) (2kr2 mgl cos θ e )(θ θ 0 ) 2 + higher order terms 3 Autonomous means that time does not appear explicitly on the right sides of the equations. 4 The first three terms are given by f(x) = f(x 0 ) + f (x 0 )(x x 0 ) f (x 0 )(x x 0 ) 2, where primes indicate first and second derivatives with respect to x.
5 Chaotic Pendulum CP 5 Why is the first order term in (θ θ 0 ) missing? (b) Show that keeping only up to the quadratic term, (and taking Γ = 0), the equations of motion can be written θ + Γ θ + Ω 2 0θ = ɛ cos Ωt (16) where a double dot over a variable represents a second derivative with respect to time, θ = θ θ 0 is the angular displacement from equilibrium, and Ω 0 = κ µ cos θ e (17) is the natural resonance frequency. As discussed below, Eq. 16 is the standard driven harmonic oscillator equation. Although it is required for chaotic motion, the pendulum mass m can be removed, turning the apparatus into a good example of a linear oscillator. The apparatus is not perfectly linear because of the nonzero Γ. The equation of motion with m = 0 (and Γ = 0) can be expressed θ + Γ θ + Ω 2 0θ = ɛ cos Ωt (18) where Ω 2 0 = κ, and Γ, κ, and ɛ are given by Eqs. 10, 12 and 14 (with I = I 0 ). Exercise 5 (a) Show that for ɛ = 0 (undriven system), Eq. 18 has an underdamped (Γ/2 < Ω 0 ) solution of the form where θ = Ce 1 2 Γt cos(ω t + δ) (19) Ω = Ω 2 0 Γ 2 /4 (20) is the oscillation frequency (pulled slightly below the natural resonance frequency Ω 0 for weak damping). Although you are not asked to do so here, C and δ can be determined from the initial conditions on θ and θ at t = 0. (b) Show that for a nonzero ɛ (driven system), the solution to Eq. 18 is the sum of Eq. 19 (called the homogeneous or transient solution) and a particular or steady state solution given by where and C = θ = C cos(ωt + δ ) (21) ɛ (Ω 2 0 Ω 2 ) 2 + Ω 2 Γ 2 (22) tan δ = ΩΓ Ω 2 Ω 2 0 Non-linear Dynamics (23) For chaotic motion, there are no analytic solutions, no closed form expressions for θ vs. t like Eqs. 19 and 21. So how can the measurements be analyzed? The following few sections present some of the techniques for representing, analyzing, and comparing chaotic data sets. Phase space The Poincaré-Bendixson theorem states that at least three time dependent (dynamical) variables are required for a system to display chaotic behavior. Our system is at this minimum as demonstrated by the three equations of motion, Eqs. 7-9, one for each variable. Consider carefully how θ, ω, and φ are the only time-dependent variables in the system, how they are independent from one another, and how they completely describe the state of the system at any particular point in time. θ Angle of the pendulum. ω How fast and in what direction the pendulum is moving. φ Angle of the drive shaft.
6 CP 6 Advanced Physics Laboratory The triplet of values (θ, ω, φ) can be considered a point u in a three dimensional coordinate system called phase space. u = θ ω φ (24) ω θ φ As the three variables change continuously in time, the phase point describing the system moves continuously through phase space. The path through phase space is called a trajectory and a plot of a trajectory is called a phase plot. The trajectory is predicted by Eqs. 7-9 which can be expressed compactly using the vector notation where u = G(u) (25) u = θ ω φ (26) and G(u) is a 3-vector function of the three dynamic variables. G(u) = For our system G θ (θ, ω, φ) G ω (θ, ω, φ) G φ (θ, ω, φ) (27) G(u) = (28) ω Γω Γ sgn ω κθ + µ sin θ + ɛ cos φ Ω For any point u in phase space, the equation u = G(u) predicts the value of u a short time later. For a small enough time step τ, the future phase can be predicted in terms of the present phase via the definition of the derivative, u = (u(t + τ) u(t))/τ, as u(t + τ) = u(t) + τ G(u(t)) (29) Figure 3: The three-dimensional phase space for the driven nonlinear oscillator. Each rectangle represents a θ-ω plane at a different drive phase φ. A darkened plane is shown for one value of φ. To determine u over a longer time interval, this equation would be iterated using a small enough time step τ. With certain refinements to improve accuracy and speed, this is how computer simulations and differential equation solvers work. The phase space coordinate system can be represented as a continuous set of θ-ω planes encircling a common axis as shown in Fig. 3. Each plane is in one-to-one correspondence with the drive phase φ. One turn of the drive system can be thought of as forcing one orbit through the planes. C.Q. 1 For undriven motion, φ is not a dynamical variable for the system and phase space becomes two-dimensional a single θ-ω plane. Show that the trajectory for undamped, undriven harmonic oscillations lies on an ellipse in the θ-ω phase plane. What relative axes scaling is needed to make the ellipse circular? What determines the direction of motion around the ellipse? What determines the time it takes to go around the ellipse? How do the trajectories depend on the initial conditions?
7 Chaotic Pendulum CP 7 Hand sketch at least two different ellipses on the same θ-ω plane. C.Q. 2 How does the addition of damping change the trajectories above? How do the trajectories depend on the initial conditions? Hand sketch at least two different trajectories on the same θ-ω plane. Where on the phase plane do all trajectories, no matter the initial conditions, ultimately end up? C.Q. 1 demonstrates that oscillations at constant amplitude have a circular trajectory in an appropriately-scaled θ-ω phase plane, while C.Q. 2 demonstrates that the trajectory spirals into the point (0, 0) when damping is present. What do trajectories look like in a driven system, for which phase space is three dimensional like that of Fig. 3? Figure 4(a) shows a plot of θ vs. t for a driven harmonic oscillator. The fast oscillations at the beginning of the curve correspond to the transient behavior of Eq. 19. These fast oscillations are superimposed on the slower steady state oscillations corresponding to Eq. 21. In plot (b), this motion is shown as a three dimensional phase space trajectory. Starting at the endpoint on the left, the fast transient oscillations become the corkscrew motion about the circular pattern, on which the trajectory repeatedly cycles and representing the steady state oscillations. The steady state orbit is not quite circular, but rides on the surface of a torus as shown in plot (c). With proper scaling, a cross section of the torus would have the twodimensional circular trajectory of C.Q. 1. But now, while the trajectory cycles around the circle in the θ-ω planes, it simultaneously cycles through the drive phase φ around the long direction of the torus. Of course, other trajectories are also possible. For example, the natural frequency Ω 0 might be closer to or lower than the drive frequency Ω rather than 20 times higher as in Figs. 4(a) and (b). Or, the damping might be weaker or stronger. These variations require a change in the system parameters. The only variations possible in a given system depend on how the motion is started and determine whether the transient motion starts off large or small. In contrast, the steady state trajectory does not depend on initial conditions it is an invariant of the system. It is an invariant because the amplitude and phase relations (Eqs. 22 and 23) specify a single closed curve through phase space. Because the phase point repeatedly cycles around a closed curve, the motion is called a limit cycle. The motion repeats after a single drive period because the steady state oscillation frequency within the θ-ω planes is the same as the drive frequency through them. Thus, the closed curve is a single loop and called a period-1 limit cycle. Not possible with a linear oscillator, our apparatus can display motion that repeats only after the drive motor has gone around five (or more) times. This would be called a period-5 (or higher) limit cycle. Attractors The single invariant curve representing the steady state motion of a driven harmonic oscillator is called an attractor for that (three dimensional) system. The final ending point (0, 0) for the damped undriven motion of C.Q. 2 is an attractor for that (two dimensional) system. Attractors are the ultimate fate of the trajectories after transients have decayed away. Nonlinear oscillators also have attractors. For our system, the type and number of attractors that occur depend on the drive frequency Ω, its amplitude ɛ, and the other system parameters κ, µ, Γ, and Γ. At low drive amplitudes, only two steady state motions are possible near harmonic os-
8 CP 8 Advanced Physics Laboratory D E F Figure 4: (a) θ vs. t for a linear oscillator. (b) The same data represented as a phase space trajectory. (c) The dark curve represents the steady state trajectory (attractor) for a linear oscillator. It lies on a torus. cillations centered about one or the other equilibrium point. Each of these motions corresponds to a period-1 limit cycle attractor. An interesting question not relevant to linear oscillators is to determine which of the two attractors the motion will end on. The initial conditions do, in fact, unambiguously determine which attractor the motion will end on. Thus, a basin of attraction can be associated with each attractor, and is defined as those phase points that lead to motion on that attractor. For nonlinear systems, the shapes of these basins can be quite complex. At larger drive amplitudes, chaotic motion can occur for which the trajectory never repeats. The phase point cycles through the full set of θ-ω planes once per drive period yet never hits the same point twice. Moreover, an attractor still describes the ultimate motion. As with the linear oscillator, arbitrary initial conditions cause motion off the attractor, but after this transient motion decays away, the system phase point ultimately moves on a subset of points in phase space. The subset is the attractor, and determining its shape and properties are important parts of the analysis. One trivial property is that it must lie within a bounded three dimensional volume of phase space. Under most operating conditions the bounds for θ are θ < 2π (the pendulum never goes over the top twice in the same direction). The range of φ is also bounded; its extents can be taken from 0 to 2π. And ω must be bounded as well. The kinetic energy Iω 2 /2 will not grossly exceed the largest possible potential energy k(rθ) 2 with θ 2π. Just as a point in phase space represents the state of a system, a volume in phase space represents a set of states. For every value of (θ, ω, φ) in that volume, the relation u = G(u) provides unique values for the derivatives of each variable. The derivatives specify a unique flow for each trajectory point so that as the trajectories evolve, they can converge toward one another and they can diverge away from one another but they can never cross. A crossing would require two different values for u at the same value of u, forbidden under the assumption that the pendulum is governed by u = G(u). As outlined in the appendix on the Lyapunov exponent algorithm, one can prove that any finite volume of phase points simultaneously evolving according to Eq. 25 (with Eq. 28) will exponentially decay to zero with a rate constant given by Γ. Thus, consider a set of initial conditions filling a finite vol-
9 Chaotic Pendulum CP 9 ume of phase space that completely covers the attractor. As the transients decay away, the phase points must all come to lie on the attractor while at the same time the volume of the attractor must exponentially decay to zero. This is possible if the attractor consists of any number of limit cycles because a limit cycle is a finite length curve through phase space and thus has zero volume. However, limit cycles are not the only possibility. In fact, chaotic motion never repeats. Consequently, the attractor on which it lies can not be a limit cycle and must be a curve with an infinite length. Thus, the chaotic attractor must be a set of points with an infinite length but zero volume. How can this be? A two dimensional surface of arbitrary shape in phase space is a set of points with this property, but it will turn out that our attractor is not any kind of oddly shaped area in phase space. You will determine that the dimensionality of our attractors is between that of an area (a two dimensional object) and a volume (a three dimensional object). Chaotic attractors apparently have strange shapes with fractional dimensions (fractals). Fractal attractors are called strange attractors. Poincaré sections We will only be able to measure a subset of the complete attractor a very long trajectory. Overnight runs of 50,000 drive cycles are typical for our apparatus. Such a long trajectory is nearly enough to completely represent the attractor and we will often not distinguish between the trajectory and the attractor. However, you should always keep in mind that the full attractor has infinitely more points than any finite trajectory. A long trajectory is difficult to display on a two-dimensional medium such as a computer screen or a piece of paper. One way to achieve Figure 5: A Poincaré section from the apparatus of Fig. 1. Figure 6: A sequence of Poincaré sections arranged to illustrate the complete three dimensional attractor in the phase space coordinate system of Fig. 3 a compact representation of a very long trajectory is with a Poincaré surface of section (Poincaré section for short) so named because it is a cross section of the full threedimensional trajectory. In a Poincaré section, θ and ω are recorded as a single point each time the trajectory pierces through the θ-ω plane corresponding to one value of the drive phase φ. Thus, once per drive period, a new (θ, ω) point is added to the plot as the trajectory repeatedly pierces that particular plane. A Poincaré section from our apparatus is shown in Fig. 5. A collection of Poincaré sec- ω θ
10 CP 10 Advanced Physics Laboratory tions at sequential drive phases φ is a representation of the complete attractor. Figure 6 shows a three dimensional arrangement of 20 Poincaré sections illustrating the attractor. Several sequences of 40 Poincaré sections can be seen as movie loops on the laboratory web page. Fractal dimension The attractor for chaotic motion is called strange because it has a strange shape called a fractal. A fractal is a self-similar object, meaning it looks similar on ever decreasing scales. One invariant of a strange attractor is its fractal dimension. Since our attractor is dimensionally the same as an extruded Poincaré section, its fractal dimension is one plus the fractal dimension of any of its Poincaré sections. Consequently, it is the fractal dimension of a Poincaré section that will be sought. While there are several different varieties of fractal dimensions, the most common is the capacity dimension. The box-counting algorithm used to determine the capacity dimension also demonstrates the self similarity of a fractal Poincaré section. The algorithm for a two dimensional fractal is as follows. Determine the bounds in the plane over which the fractal extends. Divide that area into M rm grid boxes (r = 1 is fine but not required). Count the number of grid boxes N covering any part of the fractal. Repeat the procedure as M is increased, making the grid boxes smaller and smaller. The slope of the log N vs. log M graph (in the limit as M ) is the capacity dimension. The procedure above is for the idealized case of an infinite trajectory a Poincaré section from the complete attractor. It will need to be modified for an experimental Poincaré section obtained from a finite trajectory. The phrase covering any part of the fractal is modified to having at least one phase point. The phrase in the limit as M must be dropped because if the size of boxes become too small, N will saturate at the total number of phase points in the Poincaré section. Thus, one must simply check that the graph of log N vs. log M appears to be a straight line before N saturates. For our 50,000-point Poincaré sections, r = 1/2 (twice as many grid lines in θ as in ω) and M in a range from 10 to 300 or so works well. C.Q. 3 The ordinary geometric dimensions of a point, line, area, and volume are 0, 1, 2, and 3, respectively. (a) Explain how the boxcounting algorithm would give a dimension of 0, 1, and 2, respectively, for a very large number of phase points distributed uniformly at a point, along a line, and over an area in the θ- ω plane. (b) When there are two or more such objects, e.g., points and lines or lines and areas, explain why the algorithm should give the dimension of the highest dimensional object in the set. Lyapunov exponents Fractal attractors are one hallmark of chaotic behavior. Extreme sensitivity to initial conditions is another. It is typified by an exponential growth in the phase-space separation of two trajectories which start out very near one another. Actually, at each point u on the attractor there are three different directions (global Lyapunov vectors) ê i that characterize three different exponential time dependencies that occur. Any phase point u near u having a small difference δu = u u (30) that is parallel to one of these directions, evolves with a single exponential time dependence characterized by one of three Lyapunov
11 Chaotic Pendulum CP 11 exponents. In other words, any initial separation δu(0) in the direction of a particular Lyapunov vector ê i evolves according to δu(t) = δu(0) e λ it (31) where λ i is the corresponding Lyapunov exponent for the particular direction ê i. The λ i can be positive, negative, or zero. A negative value corresponds to a direction in which the separation decays over time. A zero value corresponds to a direction in which the separation remains constant. And for chaotic motion, at least one of the λ i must be positive corresponding to a direction in which the separation grows. No matter how small, any two initial conditions will always have some component along the expanding direction, and the trajectories evolving therefrom will always diverge. Small random forces and torques acting on the pendulum (e.g., from the pickup of vibrations in the pendulum supports and from small variations in friction), cause small perturbations to the phase space trajectory. These perturbations then also evolve according to the system dynamics and are called dynamical noise. Extremely close phase points, which might take a long time to separate in a noiseless system, do so more quickly in the presence of dynamical noise. A study of the effects of dynamical noise would make an interesting advanced project. However, the dynamical noise in our apparatus is small enough that many predictions made assuming it is negligible will still be valid. One of the three Lyapunov directions and its corresponding exponent is essentially trivial and easily visualized. Imagine two identical trajectories differing only in that one leads the other by some small time interval δt. For small enough δt, Eq. 25 implies their phases will differ by δu = G δt, and although δu changes as the two systems evolve (because G depends on u), there is no average increase or decrease with time. Thus, this propagation direction has a corresponding Lyapunov exponent of zero. All dynamical systems governed by an autonomous set of differential equations will have such a direction and at least one Lyapunov exponent will be zero. Because the motion in the propagation direction is basically trivial, the behavior of nearby points on the same Poincaré section provide all the essential information. The general behavior is illustrated in the two Poincaré sections shown in Fig. 7. Consider the single phase point u n marked by the cursor in the section on the left. The highlighted points centered around u n are a set of nearest neighbors u j found by searching through all points on this Poincaré section and selecting those inside an elliptical region satisfying ( ) 2 ( ) 2 δθj δωj + < 1 (32) θ r where δθ j = θ n θ j, δω j = ω n ω j. The quantities θ r and ω r specify the major axes of the elliptical neighborhood (0.15 rad and 0.30 rad/s, respectively, for the section on the left in Fig. 7). For the point u n and the set of points u j, a corresponding point u n+1 and set u j+1 are also thereby determined as those points one time step later. These points are shown in the section on the right for a time step τ equal to 1/5th of the drive period. From these two sets are constructed two sets of deviations δu n = u n u j and δu n+1 = u n+1 u j+1 one time step forward. For small values of δu, the transformation between the δu in each neighborhood is predicted to be relatively simple rotations and/or contractions and/or expansions. While the rotations can not be discerned in Fig. 7, a contracting and expanding direction are evident. (The scaling is the same ω r
12 CP 12 Advanced Physics Laboratory Figure 7: A small area of a Poincaré section showing how an initial set of δu n around the cursor in the left plot move and change shape when propagated 1/5 period. The scale ranges are the same in the two figures and the evolved set δu n+1 is the elongated and flattened patch around the cursor in the right figure. in the two sections.) The basic transformation can be modeled as a 2 2 matrix multiplication ( ) ( ) ( ) δθn+1 b11 b = 12 δθn (33) δω n+1 b 21 b 22 δω n with the matrix elements b ij depending on where u is located on the attractor. A least squares fit of the data (the δu n and δu n+1 ) to Eq. 33 provides best estimates of the transformation matrix elements. 5 The matrix elements are then used to determine two local Lyapunov exponents which must be averaged along a long trajectory to determine the global Lyapunov exponents characteristic of the attractor as a whole. From an arbitrary starting point, this procedure finding the nearest neighbors, determining the δu n and δu n+1, performing the fit, and determining the two local Lyapunov exponents is repeated at sequential points along a trajectory. The local Lyapunov exponents are monitored and 5 As shown in the Appendix on Lyapunov exponents a modified version of Eq. 33 is actually fit to the data. the calculation continues until their average is deemed to have stabilized. A more complete description of the algorithm and its theoretical underpinnings is provided in the Lyapunov Exponent addendum that can be found on the laboratory web site. Fit of the angular acceleration The Savitsky-Golay filters can be used to calculate the angular acceleration at every point along a trajectory. Modeling the values of α (= ω) according to Eq. 8, a linear regression can be used to experimentally determine Γ, Γ, κ, µ, and ɛ. For the fit, it is important to take into account the differences between the measured angles and those of the theory. For example, θ = θ m +δθ would be used to describe the fact that the computer value θ m may be offset from that of the theory θ. Similarly, φ = φ m + δφ. Lastly, the possibility that the spring equilibrium may not be set exactly at θ = 0 should be taken into account; it is taken to occur at
13 Chaotic Pendulum CP 13 some unknown angle θ 0. With these offsets in mind, Eq. 8 becomes α = Γω Γ sgn ω κ(θ m + δθ θ 0 ) (34) +µ sin(θ m + δθ) + ɛ cos(φ m + δφ) This equation can be put in a form suitable for regression analysis using the trigonometric identities: sin(a + b) = sin a cos b + sin b cos a and cos(a + b) = cos a cos b sin a sin b. C.Q. 4 (a) Show that the regression expression can be written: α = C[1] + G[ω] + G [sgn ω] + K[θ m ] +M s [sin θ m ] + M c [cos θ m ] (35) +E c [cos φ m ] + E s [sin φ m ] where each regression term is in brackets and where the coefficients of each term are given by C = κ(δθ θ 0 ) (36) G = Γ G = Γ K = κ M s = µ cos δθ M c = µ sin δθ E c = ɛ cos δφ E s = ɛ sin δφ (b) After obtaining these coefficients from the fit, they are used to determine the parameters. Show that this becomes: Γ = G (37) Γ = G κ = K µ = Ms 2 + Mc 2 δθ = tan 1 M c ɛ = M s E 2 c + E 2 s δφ = tan 1 E s E c θ 0 = tan 1 M c M s C K While the Savitsky-Golay filters can be quite accurate when fitting continuous functions, the discontinuous behavior of axle friction (it changes sign at ω = 0) presents a slight problem. Because the filters average over a finite time, they smooth the discontinuity over a range that depends on the number of points used in the filter. For the recommended 33- point filters, the axle friction smoothing appears to be over a range of ±3 rad/s. To prevent systematic errors, the fits should exclude data points in this range. LabVIEW Programs The following sections describe the use of the data acquisition and analysis programs available for this experiment. All analysis programs are limited to three-dimensional attractors like those for the system described here. The analyses typically begin by converting the data sets to a discrete set of two-dimensional Poincaré sections equally spaced along the trivial (drive phase) dynamical variable. The analysis programs are designed to work with several kinds of data sets that can be generated in our lab. In the programs they are labeled: Poincaré: A text-type spreadsheet file whose first two columns contain a single Poincaré section. I16: A binary file containing a single row of 16-bit signed integers. This is the correct file type for the chaotic pendulum. datalog: A binary file containing multiple rows of double precision numbers.
14 CP 14 Advanced Physics Laboratory This is the type of file made by the simulation programs. N I16: A binary file containing multiple rows of 16-bit signed integers. This is the type of file made when collecting data for a driven electronic circuit. The first two indices in the var. order control then determine which variables are used in the Poincaré sections. All of the programs are started by clicking on the Run button on the toolbar. Different analysis programs may run with only particular kinds of data sets. The File type tab control on the front panel of the program will have a tab for each allowed type and the type desired must be selected before hitting the Run button. For all but the DisplayI16 program for real time display of experimental data, various data set parameters (described next) must also be selected before running the program. For the I16 file type there is a delta t control giving the time between each measured integer. It is used with the Savitsky-Golay filtering to determine properly scaled derivatives. For the chaotic pendulum it should be set to 1/200 of the drive period. For both the I16 and N I16 file types there is a delta y control (array for the latter) which gives the physical step size for each change of 1 in the 16-bit integer. For the chaotic pendulum, it should be set to 2π/1440, the default setting. Savitsky-Golay filtering is used to determine θ, ω = θ, and α = θ values from the raw rotary encoder readings. 6 This filtering is equivalent to performing an unweighted, linear least squares fit of the readings to a polynomial at each point. The control for filter averaging (pts) is used to set the number of consecutive data points used (33 is recommended) and the polynomial order is set by the polynomial order 6 See the laboratory chaos web page for addendum on the theory of Savitsky-Golay filters. control (a 4th order or quartic polynomial is recommended). The programs may have from and to controls and/or a Slice control. The from and to controls determine which drive cycles in the data set will be analyzed. For a data set with 50,000 drive cycles, choosing the from and to to be 5001 and would cause the program to use the ten thousand cycles starting from cycle Programs requiring multiple Poincaré sections from the attractor will have a sections/cycle control. These programs use a set of n equally-spaced Poincaré sections, where n is the setting for this control. Programs that use only a single Poincaré section will have an additional slice control to select which of the n Poincaré sections the program will use. Both kinds of programs must also know the total number of data points per drive cycle and will have a readings/cycle control, which should be set for 200 for experimental chaotic data (I16 files) or 40 for simulations (default-sized datalog files). Acquire chaotic pendulum data Raw data from the rotary encoder (1440 pulses per revolution) are collected and saved to a file by the Acquire chaotic pendulum data program. This program steps the (200 pulse-perrevolution) drive motor at a constant rate and triggers a reading of the rotary encoder with each motor pulse. The drive frequency (in rotations per second) is set by the stepper motor frequency control. After the Run button in the toolbar is pressed, the control is updated with the nearest frequency the system was able to reach. This program will ask for a file name in which it will store the data. Navigate to the My Documents area to store the data. Unless this is a run you are sure you will want to save, keep
15 Chaotic Pendulum CP 15 the default name Data.bin and click OK when it asks if it can be overwritten. You can always rename it afterward if you decide it will be needed later. The reset button can be pressed while taking data to restart saving data at the beginning of the file. However, previous data remains part of the file until overwritten. You can change the drive frequency two ways. The automatic method repeatedly changes the frequency by a specific amount while continuously taking data and it should only be used for independent study. At the bottom of the front panel there is a button to turn the frequency sweep feature on and off and a control for the frequency change per step (negative values are OK). Each step is applied after a fixed number of pulses to the stepper motor (step delay). The default value gives a sweep of Hz/cycle. With the frequency sweep off, the frequency can be changed manually and on the fly while data collection continues. To do so, simply change the value in the stepper motor frequency control and then click on change frequency button. Keep in mind that unless the drive frequency is constant for the entire run, any programs using the delta t control (described later) will not have correct time or frequency scales. All information about the run, such as drive frequency, the magnet damping distance and the drive shaft length, should be written in your lab notebook. The angular position of the drive shaft and of the rotary encoder when the toolbar Run button is pressed are taken to be φ = 0 and θ = 0, respectively, in the data analysis programs. Consequently, if you want φ = 0 and θ = 0 to be close to the straight up positions, you will have to set the drive shaft and hold the pendulum in these positions as you press the Run button. To set the drive shaft, turn the motor controller off, adjust the shaft by hand to the vertical position, and turn the motor controller back on before pressing the Run button. However, these zeroing procedures are hardly ever relevant. Just keep in mind that in any analysis, the drive shaft and pendulum angle for any data set will generally have an unknown offset. The offsets are irrelevant for most analyses, and when they are relevant, for example in the fits to the angular acceleration, they are included as parameters. For collecting data on undriven motion, set the drive frequency to 1.0 Hz. This will establish the rotary encoder data collection rate at 200 readings per second. Set the drive shaft for A = 0 and shut off the motor controller to prevent the motor from turning. DisplayI16 Use the DisplayI16 program to look at the data while it is being collected, to look at short data sets, or to make ASCII files of the data suitable for spreadsheets. Before hitting the toolbar Run button, enter values for the readings/cycle (200), delta y (2π/1440) and delta t (T/200, where T is the drive period). The history length (cycles) control (described shortly) must also be selected before running the program. After hitting the Run button, the program asks for a data file. Specify the file name given for the Acquire chaotic pendulum data program Data.bin or any previously collected binary data files. However, depending on the file and buffer size, only points from the end of the data set may be available to look at. The topmost raw readings graph displayed on the front panel is the raw count from the rotary encoder vs. the cycle number. The history length (cycles) control sets the number of drive cycles N that will be available for use by this program. As the data set grows, only the most recent N cycles will be available. The button labeled real time updating should be set
16 CP 16 Advanced Physics Laboratory to on to update the buffer with new data as it is being collected. New data is then added to the end of the buffer as old data is removed from the beginning. Click real time updating to off to stop reading in new data and to make it easier to select data for further processing. When you later click it back on, all unread data will then be read in. Set the cursors on the raw readings graph to cover the desired data and click on the update graphs button to the left of this graph to display the selected data on the other three graphs on the front panel or to save it to a spreadsheet file. The leading cursor determines the first data point processed, while the sections/cycle determine the subsequent data points processed until such point would be after the trailing cursor. (For example, at 40 sections/cycle, every fifth point is processed.) Each of the other three graphs and the save x, v, a to spreadsheet have a button next to them labeled update or don t update. For the graphs, the state of these buttons determines whether the graphed data is updated according to the new settings or left unchanged. For the save θωα to spreadsheet the update/don t update determines whether the processed data will be saved to a text file readable by most spreadsheet programs. The format %.7e can be changed. Change the 7 to the number of digits to the right of the decimal point. The graph directly below the raw readings is a plot of x (θ), v (ω), and/or a (α) vs. t, as chosen by the check boxes to the left of the graph. The time interval between points is calculated assuming the raw acquisition rate is constant 200 times the (starting, saved) motor frequency. Thus, if frequency scanning were used, the time scales on these graphs will be incorrect. To the right of the vs. t graphs is the Fourier power spectrum (associated with x vs. t). Below the vs. t graphs is a graph of ω vs. θ. Set the Readings/rev to 1 and the plot style to points/no line to get a Poincaré section. Set the Readings/rev anywhere from and change the plot style to line/no points for a phase plot. (To set the plot style click in the legend box at the upper left of the graph. The Common Plots options let you set whether points, lines, or both are displayed.) Chaotic pendulum fit The algorithm is described in the theory section. The program needs to be told the data for the fit. You do not want to fit too big a data set. It will bog down the computer. Fitting about ten thousand drive cycles at 10 sections per cycle works well and goes reasonably fast. You can try fitting data from the beginning, middle, and end of the data set (with the from and to controls) to see if the parameter estimates are stable throughout the run. The parameters displayed on the right of the front panel are those of Eq. 34. However, the displayed covariance matrix is that of the fitting parameters of Eq. 35. It is a matrix of the sample covariances s 2 ij (diagonal elements are the square of the estimated parameter uncertainties), with the rows and columns of the matrix in the following order: G, K, M s, M c, E s, E c, C, G. The fit residuals are the differences between α as determined by the Savitsky-Golay filters and the value obtained through Eq. 8 using the fitted parameter values. These can be plot versus θ, ω, or φ, and the excluded points can be included in this plot or not. Box counting dimension The algorithm is described in the theory section. Recall that the fractal dimension it will find is that of one Poincaré section from the attractor as selected with the Slice control. After clicking Run on the toolbar, the program determines the chosen Poincaré section and then
17 Chaotic Pendulum CP 17 starts histogramming it. The histogramming consists of setting up a phase space grid and determining the number of measured phase points in each grid box. The overall phase space area it will grid up is determined by the minimum and maximum values of θ and ω found in the data. Thus, you should autoscale the graph of the Poincaré section to make sure it does not contain any phase points far from the main attractor which would cause the program to cover an inappropriately large area. The program repeatedly performs the histogramming as it decreases the grid box size (increases the number of grid boxes). The grids are determined by the array labeled grid size M. For each entry in this array, a phase space grid is formed with M θ-intervals and rm ω-intervals, with the ω-factor (r) as given in the control for this parameter (0.5 is the default). Thus the first default grid (M = 10) is a 10 5 grid and thus has 50 grid boxes. The final M = 400 value produces a grid of 80,000 boxes. The Probability density shown below the Poincaré section is an intensity graph with the z-value indicated by the color and is proportional to the number of phase points per unit area in each grid box. Looking at this graph, you will see the grid size decrease as the program proceeds. As each histogram is determined, the cells occupied N indicator displays the number of grid boxes in which at least one phase point was found. This data is plotted on a log N (ordinate) vs. log M (abscissa) plot. As described in the theory section, there are good reasons why the value of N for the smallest and largest values of M may not behave as the fractal scaling would predict. Use the cursors to select the starting and ending points of an appropriate region to use in the fit to the straight line. Then click on the compute slope button. The slope of this line is the capacity dimension. Besides the capacity dimension, two common fractal dimensions are the information dimension and correlation dimension. In all three cases, the dimension is the slope of a fitted line with log M as the abscissa. For the capacity dimension, each grid box contributes 1 (if there is a phase point in it) or 0 (if it is empty), and the ordinate is the log of the sum of the contributions. For the information dimension, each grid box contributes f i log f i, where f i is the fraction of the phase points in the grid box i, and the graph ordinate is the sum of this quantity over all grid boxes. For the correlation dimension, each grid box contributes fi 2, and the ordinate is the log of the sum of the contributions from all grid boxes. These other two dimensions can be selected by the Dimension type control. Lyapunov The algorithm is described in the theory section and more fully in the Lyapunov Exponent Addendum on the laboratory web site. This program reads the data and makes a set of equally spaced Poincaré sections as determined by the Readings/rev control. (It always uses slice 0 as the first Poincaré section.) Ten sections (i.e., slices) have been found to be enough. No check has been made for the minimum number of drive cycles needed, but 40,000-50,000 is enough to get good results. Use more if you have them, but the computer may start to get bogged down if you try to use too many slices, or too many drive cycles. The size of the ellipse for selecting δu n around each point u n along the trajectory is determined by the θ-radius and ω-radius controls (θ r and ω r in Eq. 32). They should be set to produce a near-circular area and be as small as possible while big enough that the area will nearly always have more than some
Addendum: Lyapunov Exponent Calculation
Addendum: Lyapunov Exponent Calculation Experiment CP-Ly Equations of Motion The system phase u is written in vector notation as θ u = (1) φ The equation of motion can be expressed compactly using the
More informationMechanical Resonance and Chaos
Mechanical Resonance and Chaos You will use the apparatus in Figure 1 to investigate regimes of increasing complexity. Figure 1. The rotary pendulum (from DeSerio, www.phys.ufl.edu/courses/phy483l/group_iv/chaos/chaos.pdf).
More informationUpdated 2013 (Mathematica Version) M1.1. Lab M1: The Simple Pendulum
Updated 2013 (Mathematica Version) M1.1 Introduction. Lab M1: The Simple Pendulum The simple pendulum is a favorite introductory exercise because Galileo's experiments on pendulums in the early 1600s are
More informationChaotic pendulum: The complete attractor
Chaotic pendulum: The complete attractor Robert DeSerio a) Department of Physics, University of Florida, Gainesville, Florida 32611 Received 16 April 2002; accepted 11 October 2002 A commercial chaotic
More informationPhysics 106b: Lecture 7 25 January, 2018
Physics 106b: Lecture 7 25 January, 2018 Hamiltonian Chaos: Introduction Integrable Systems We start with systems that do not exhibit chaos, but instead have simple periodic motion (like the SHO) with
More informationChaotic Motion (One or two weights)
Page 1 of 8 Chaotic Motion (One or two weights) Exercises I through IV form the one-weight experiment. Exercises V through VII, completed after Exercises I-IV, add one weight more. This challenging experiment
More informationLab 11 - Free, Damped, and Forced Oscillations
Lab 11 Free, Damped, and Forced Oscillations L11-1 Name Date Partners Lab 11 - Free, Damped, and Forced Oscillations OBJECTIVES To understand the free oscillations of a mass and spring. To understand how
More informationMAS212 Assignment #2: The damped driven pendulum
MAS Assignment #: The damped driven pendulum Sam Dolan (January 8 Introduction In this assignment we study the motion of a rigid pendulum of length l and mass m, shown in Fig., using both analytical and
More informationLab M1: The Simple Pendulum
Spring 2003 M1.1 Introduction. Lab M1: The Simple Pendulum The simple pendulum is a favorite introductory exercise because Galileo's experiments on pendulums in the early 1600s are usually regarded as
More informationPH 120 Project # 2: Pendulum and chaos
PH 120 Project # 2: Pendulum and chaos Due: Friday, January 16, 2004 In PH109, you studied a simple pendulum, which is an effectively massless rod of length l that is fixed at one end with a small mass
More informationCHAOS -SOME BASIC CONCEPTS
CHAOS -SOME BASIC CONCEPTS Anders Ekberg INTRODUCTION This report is my exam of the "Chaos-part" of the course STOCHASTIC VIBRATIONS. I m by no means any expert in the area and may well have misunderstood
More informationPreLab 2 - Simple Harmonic Motion: Pendulum (adapted from PASCO- PS-2826 Manual)
Musical Acoustics Lab, C. Bertulani, 2012 PreLab 2 - Simple Harmonic Motion: Pendulum (adapted from PASCO- PS-2826 Manual) A body is said to be in a position of stable equilibrium if, after displacement
More informationInverted Pendulum System
Introduction Inverted Pendulum System This lab experiment consists of two experimental procedures, each with sub parts. Experiment 1 is used to determine the system parameters needed to implement a controller.
More informationChapter 14 (Oscillations) Key concept: Downloaded from
Chapter 14 (Oscillations) Multiple Choice Questions Single Correct Answer Type Q1. The displacement of a particle is represented by the equation. The motion of the particle is (a) simple harmonic with
More informationChaotic motion. Phys 750 Lecture 9
Chaotic motion Phys 750 Lecture 9 Finite-difference equations Finite difference equation approximates a differential equation as an iterative map (x n+1,v n+1 )=M[(x n,v n )] Evolution from time t =0to
More informationOscillatory Motion. Simple pendulum: linear Hooke s Law restoring force for small angular deviations. Oscillatory solution
Oscillatory Motion Simple pendulum: linear Hooke s Law restoring force for small angular deviations d 2 θ dt 2 = g l θ θ l Oscillatory solution θ(t) =θ 0 sin(ωt + φ) F with characteristic angular frequency
More informationTransitioning to Chaos in a Simple Mechanical Oscillator
Transitioning to Chaos in a Simple Mechanical Oscillator Hwan Bae Physics Department, The College of Wooster, Wooster, Ohio 69, USA (Dated: May 9, 8) We vary the magnetic damping, driver frequency, and
More informationOscillatory Motion. Simple pendulum: linear Hooke s Law restoring force for small angular deviations. small angle approximation. Oscillatory solution
Oscillatory Motion Simple pendulum: linear Hooke s Law restoring force for small angular deviations d 2 θ dt 2 = g l θ small angle approximation θ l Oscillatory solution θ(t) =θ 0 sin(ωt + φ) F with characteristic
More informationDriven Harmonic Oscillator
Driven Harmonic Oscillator Physics 6B Lab Experiment 1 APPARATUS Computer and interface Mechanical vibrator and spring holder Stands, etc. to hold vibrator Motion sensor C-209 spring Weight holder and
More informationExperiment 5. Simple Harmonic Motion
Reading and Problems: Chapters 7,8 Problems 7., 8. Experiment 5 Simple Harmonic Motion Goals. To understand the properties of an oscillating system governed by Hooke s Law.. To study the effects of friction
More informationThe Pendulum. The purpose of this tab is to predict the motion of various pendulums and compare these predictions with experimental observations.
The Pendulum Introduction: The purpose of this tab is to predict the motion of various pendulums and compare these predictions with experimental observations. Equipment: Simple pendulum made from string
More informationExperiment IV. To find the velocity of waves on a string by measuring the wavelength and frequency of standing waves.
Experiment IV The Vibrating String I. Purpose: To find the velocity of waves on a string by measuring the wavelength and frequency of standing waves. II. References: Serway and Jewett, 6th Ed., Vol., Chap.
More informationPhysics 2310 Lab #3 Driven Harmonic Oscillator
Physics 2310 Lab #3 Driven Harmonic Oscillator M. Pierce (adapted from a lab by the UCLA Physics & Astronomy Department) Objective: The objective of this experiment is to characterize the behavior of a
More informationLab Partner(s) TA Initials (on completion) EXPERIMENT 7: ANGULAR KINEMATICS AND TORQUE
TA name Lab section Date TA Initials (on completion) Name UW Student ID # Lab Partner(s) EXPERIMENT 7: ANGULAR KINEMATICS AND TORQUE 117 Textbook Reference: Walker, Chapter 10-1,2, Chapter 11-1,3 SYNOPSIS
More informationEXPERIMENT 7: ANGULAR KINEMATICS AND TORQUE (V_3)
TA name Lab section Date TA Initials (on completion) Name UW Student ID # Lab Partner(s) EXPERIMENT 7: ANGULAR KINEMATICS AND TORQUE (V_3) 121 Textbook Reference: Knight, Chapter 13.1-3, 6. SYNOPSIS In
More informationChaotic motion. Phys 420/580 Lecture 10
Chaotic motion Phys 420/580 Lecture 10 Finite-difference equations Finite difference equation approximates a differential equation as an iterative map (x n+1,v n+1 )=M[(x n,v n )] Evolution from time t
More informationThe Torsion Pendulum (One or two weights)
The Torsion Pendulum (One or two weights) Exercises I through V form the one-weight experiment. Exercises VI and VII, completed after Exercises I -V, add one weight more. Preparatory Questions: 1. The
More informationName Date: Course number: MAKE SURE TA & TI STAMPS EVERY PAGE BEFORE YOU START. Grade: EXPERIMENT 4
Laboratory Section: Last Revised on June 18, 2018 Partners Names: Grade: EXPERIMENT 4 Moment of Inertia & Oscillations 0 Pre-Laboratory Work [20 pts] 1 a) In Section 31, describe briefly the steps you
More informationChaotic Motion (One or two weights)
Page 1 of 8 Chaotic Motion (One or two weights) Exercises I through IV form the one-weight experiment. Exercises V through VII, completed after Exercises I-IV, add one weight more. This challenging experiment
More informationLab 1: Damped, Driven Harmonic Oscillator
1 Introduction Lab 1: Damped, Driven Harmonic Oscillator The purpose of this experiment is to study the resonant properties of a driven, damped harmonic oscillator. This type of motion is characteristic
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department Physics 8.01T Fall Term 2004 Experiment 06: Work, Energy and the Harmonic Oscillator Purpose of the Experiment: In this experiment you allow a cart
More informationLab 1: damped, driven harmonic oscillator
Lab 1: damped, driven harmonic oscillator 1 Introduction The purpose of this experiment is to study the resonant properties of a driven, damped harmonic oscillator. This type of motion is characteristic
More informationAP Physics C Mechanics Objectives
AP Physics C Mechanics Objectives I. KINEMATICS A. Motion in One Dimension 1. The relationships among position, velocity and acceleration a. Given a graph of position vs. time, identify or sketch a graph
More informationSimple Harmonic Motion - MBL
Simple Harmonic Motion - MBL In this experiment you will use a pendulum to investigate different aspects of simple harmonic motion. You will first examine qualitatively the period of a pendulum, as well
More informationLaboratory handout 5 Mode shapes and resonance
laboratory handouts, me 34 82 Laboratory handout 5 Mode shapes and resonance In this handout, material and assignments marked as optional can be skipped when preparing for the lab, but may provide a useful
More informationRotational Motion. 1 Purpose. 2 Theory 2.1 Equation of Motion for a Rotating Rigid Body
Rotational Motion Equipment: Capstone, rotary motion sensor mounted on 80 cm rod and heavy duty bench clamp (PASCO ME-9472), string with loop at one end and small white bead at the other end (125 cm bead
More informationTorsion Spring Oscillator with Dry Friction
Torsion Spring Oscillator with Dry Friction Manual Eugene Butikov Annotation. The manual includes a description of the simulated physical system and a summary of the relevant theoretical material for students
More informationChapter 16: Oscillatory Motion and Waves. Simple Harmonic Motion (SHM)
Chapter 6: Oscillatory Motion and Waves Hooke s Law (revisited) F = - k x Tthe elastic potential energy of a stretched or compressed spring is PE elastic = kx / Spring-block Note: To consider the potential
More informationDamped harmonic motion
Damped harmonic motion March 3, 016 Harmonic motion is studied in the presence of a damping force proportional to the velocity. The complex method is introduced, and the different cases of under-damping,
More informationNonlinear Autonomous Systems of Differential
Chapter 4 Nonlinear Autonomous Systems of Differential Equations 4.0 The Phase Plane: Linear Systems 4.0.1 Introduction Consider a system of the form x = A(x), (4.0.1) where A is independent of t. Such
More informationGood Vibes: Introduction to Oscillations
Chapter 14 Solutions Good Vibes: Introduction to Oscillations Description: Several conceptual and qualitative questions related to main characteristics of simple harmonic motion: amplitude, displacement,
More informationThe Coupled Pendulum Experiment
The Coupled Pendulum Experiment In this lab you will briefly study the motion of a simple pendulum, after which you will couple two pendulums and study the properties of this system. 1. Introduction to
More informationForced Oscillations in a Linear System Problems
Forced Oscillations in a Linear System Problems Summary of the Principal Formulas The differential equation of forced oscillations for the kinematic excitation: ϕ + 2γ ϕ + ω 2 0ϕ = ω 2 0φ 0 sin ωt. Steady-state
More informationOscillatory Motion SHM
Chapter 15 Oscillatory Motion SHM Dr. Armen Kocharian Periodic Motion Periodic motion is motion of an object that regularly repeats The object returns to a given position after a fixed time interval A
More informationUnit 7: Oscillations
Text: Chapter 15 Unit 7: Oscillations NAME: Problems (p. 405-412) #1: 1, 7, 13, 17, 24, 26, 28, 32, 35 (simple harmonic motion, springs) #2: 45, 46, 49, 51, 75 (pendulums) Vocabulary: simple harmonic motion,
More informationOscillations. PHYS 101 Previous Exam Problems CHAPTER. Simple harmonic motion Mass-spring system Energy in SHM Pendulums
PHYS 101 Previous Exam Problems CHAPTER 15 Oscillations Simple harmonic motion Mass-spring system Energy in SHM Pendulums 1. The displacement of a particle oscillating along the x axis is given as a function
More informationIntroduction. Pre-Lab Questions: Physics 1CL PERIODIC MOTION - PART II Fall 2009
Introduction This is the second of two labs on simple harmonic motion (SHM). In the first lab you studied elastic forces and elastic energy, and you measured the net force on a pendulum bob held at an
More informationChapter 15 Periodic Motion
Chapter 15 Periodic Motion Slide 1-1 Chapter 15 Periodic Motion Concepts Slide 1-2 Section 15.1: Periodic motion and energy Section Goals You will learn to Define the concepts of periodic motion, vibration,
More information2: SIMPLE HARMONIC MOTION
2: SIMPLE HARMONIC MOTION Motion of a mass hanging from a spring If you hang a mass from a spring, stretch it slightly, and let go, the mass will go up and down over and over again. That is, you will get
More informationChapter 14 Periodic Motion
Chapter 14 Periodic Motion 1 Describing Oscillation First, we want to describe the kinematical and dynamical quantities associated with Simple Harmonic Motion (SHM), for example, x, v x, a x, and F x.
More informationLaboratory handouts, ME 340
Laboratory handouts, ME 340 This document contains summary theory, solved exercises, prelab assignments, lab instructions, and report assignments for Lab 4. 2014-2016 Harry Dankowicz, unless otherwise
More informationPHYSICS 1 Simple Harmonic Motion
Advanced Placement PHYSICS 1 Simple Harmonic Motion Student 014-015 What I Absolutely Have to Know to Survive the AP* Exam Whenever the acceleration of an object is proportional to its displacement and
More informationPHY 123 Lab 6 - Angular Momentum
1 PHY 123 Lab 6 - Angular Momentum (updated 10/17/13) The purpose of this lab is to study torque, moment of inertia, angular acceleration and the conservation of angular momentum. If you need the.pdf version
More informationChapter 1, Section 1.2, Example 9 (page 13) and Exercise 29 (page 15). Use the Uniqueness Tool. Select the option ẋ = x
Use of Tools from Interactive Differential Equations with the texts Fundamentals of Differential Equations, 5th edition and Fundamentals of Differential Equations and Boundary Value Problems, 3rd edition
More informationSHM Simple Harmonic Motion revised May 23, 2017
SHM Simple Harmonic Motion revised May 3, 017 Learning Objectives: During this lab, you will 1. communicate scientific results in writing.. estimate the uncertainty in a quantity that is calculated from
More informationName (please print): UW ID# score last first
Name (please print): UW ID# score last first Question I. (20 pts) Projectile motion A ball of mass 0.3 kg is thrown at an angle of 30 o above the horizontal. Ignore air resistance. It hits the ground 100
More informationRaymond A. Serway Chris Vuille. Chapter Thirteen. Vibrations and Waves
Raymond A. Serway Chris Vuille Chapter Thirteen Vibrations and Waves Periodic Motion and Waves Periodic motion is one of the most important kinds of physical behavior Will include a closer look at Hooke
More informationExamining the Effect of Additional Driving Forces on a Mechanical Oscillator
Examining the Effect of Additional Driving Forces on a Mechanical Oscillator Nate Stone Physics Department, The College of Wooster, Wooster, Ohio 44691, USA (Dated: May 7, 2015) In this experiment, the
More informationBig Idea 4: Interactions between systems can result in changes in those systems. Essential Knowledge 4.D.1: Torque, angular velocity, angular
Unit 7: Rotational Motion (angular kinematics, dynamics, momentum & energy) Name: Big Idea 3: The interactions of an object with other objects can be described by forces. Essential Knowledge 3.F.1: Only
More informationMotion on a linear air track
Motion on a linear air track Introduction During the early part of the 17 th century, Galileo experimentally examined the concept of acceleration. One of his goals was to learn more about freely falling
More informationPhysics Mechanics. Lecture 32 Oscillations II
Physics 170 - Mechanics Lecture 32 Oscillations II Gravitational Potential Energy A plot of the gravitational potential energy U g looks like this: Energy Conservation Total mechanical energy of an object
More information2: SIMPLE HARMONIC MOTION
2: SIMPLE HARMONIC MOTION Motion of a Mass Hanging from a Spring If you hang a mass from a spring, stretch it slightly, and let go, the mass will go up and down over and over again. That is, you will get
More informationRotational Motion. Figure 1: Torsional harmonic oscillator. The locations of the rotor and fiber are indicated.
Rotational Motion 1 Purpose The main purpose of this laboratory is to familiarize you with the use of the Torsional Harmonic Oscillator (THO) that will be the subject of the final lab of the course on
More informationPhysics 326 Lab 6 10/18/04 DAMPED SIMPLE HARMONIC MOTION
DAMPED SIMPLE HARMONIC MOTION PURPOSE To understand the relationships between force, acceleration, velocity, position, and period of a mass undergoing simple harmonic motion and to determine the effect
More informationPhysics 101 Final Exam Problem Guide
Physics 101 Final Exam Problem Guide Liam Brown, Physics 101 Tutor C.Liam.Brown@gmail.com General Advice Focus on one step at a time don t try to imagine the whole solution at once. Draw a lot of diagrams:
More informationPhysics 106a, Caltech 4 December, Lecture 18: Examples on Rigid Body Dynamics. Rotating rectangle. Heavy symmetric top
Physics 106a, Caltech 4 December, 2018 Lecture 18: Examples on Rigid Body Dynamics I go through a number of examples illustrating the methods of solving rigid body dynamics. In most cases, the problem
More informationLab 10 - Harmonic Motion and the Pendulum
Lab 10 Harmonic Motion and the Pendulum L10-1 Name Date Partners Lab 10 - Harmonic Motion and the Pendulum L (measured from the suspension point to the center of mass) Groove marking the center of mass
More informationL = 1 2 a(q) q2 V (q).
Physics 3550, Fall 2011 Motion near equilibrium - Small Oscillations Relevant Sections in Text: 5.1 5.6 Motion near equilibrium 1 degree of freedom One of the most important situations in physics is motion
More information4. Complex Oscillations
4. Complex Oscillations The most common use of complex numbers in physics is for analyzing oscillations and waves. We will illustrate this with a simple but crucially important model, the damped harmonic
More informationSIMPLE HARMONIC MOTION
SIMPLE HARMONIC MOTION PURPOSE The purpose of this experiment is to investigate simple harmonic motion. We will determine the elastic spring constant of a spring first and then study small vertical oscillations
More information5-Sep-15 PHYS101-2 GRAPHING
GRAPHING Objectives 1- To plot and analyze a graph manually and using Microsoft Excel. 2- To find constants from a nonlinear relation. Exercise 1 - Using Excel to plot a graph Suppose you have measured
More informationIntroduction. Pre-Lab Questions: Physics 1CL PERIODIC MOTION - PART II Spring 2009
Introduction This is the second of two labs on simple harmonic motion (SHM). In the first lab you studied elastic forces and elastic energy, and you measured the net force on a pendulum bob held at an
More informationPHYS2330 Intermediate Mechanics Fall Final Exam Tuesday, 21 Dec 2010
Name: PHYS2330 Intermediate Mechanics Fall 2010 Final Exam Tuesday, 21 Dec 2010 This exam has two parts. Part I has 20 multiple choice questions, worth two points each. Part II consists of six relatively
More informationMechanics IV: Oscillations
Mechanics IV: Oscillations Chapter 4 of Morin covers oscillations, including damped and driven oscillators in detail. Also see chapter 10 of Kleppner and Kolenkow. For more on normal modes, see any book
More informationChapter 14 Oscillations
Chapter 14 Oscillations Chapter Goal: To understand systems that oscillate with simple harmonic motion. Slide 14-2 Chapter 14 Preview Slide 14-3 Chapter 14 Preview Slide 14-4 Chapter 14 Preview Slide 14-5
More informationPHYSICS. Chapter 15 Lecture FOR SCIENTISTS AND ENGINEERS A STRATEGIC APPROACH 4/E RANDALL D. KNIGHT Pearson Education, Inc.
PHYSICS FOR SCIENTISTS AND ENGINEERS A STRATEGIC APPROACH 4/E Chapter 15 Lecture RANDALL D. KNIGHT Chapter 15 Oscillations IN THIS CHAPTER, you will learn about systems that oscillate in simple harmonic
More informationPhysics I: Oscillations and Waves Prof. S. Bharadwaj Department of Physics and Meteorology. Indian Institute of Technology, Kharagpur
Physics I: Oscillations and Waves Prof. S. Bharadwaj Department of Physics and Meteorology Indian Institute of Technology, Kharagpur Lecture No 03 Damped Oscillator II We were discussing, the damped oscillator
More informationBrown University Physics 0030 Physics Department Lab 5
Oscillatory Motion Experiment 1: Oscillations of a spring As described in the text, a system of a simple spring exhibits simple harmonic motion and is a good introduction to a study of oscillations, which
More informationComputer Problems for Taylor Series and Series Convergence
Computer Problems for Taylor Series and Series Convergence The two problems below are a set; the first should be done without a computer and the second is a computer-based follow up. 1. The drawing below
More informationOscillator Homework Problems
Oscillator Homework Problems Michael Fowler 3//7 1 Dimensional exercises: use dimensions to find a characteristic time for an undamped simple harmonic oscillator, and a pendulum Why does the dimensional
More informationThe phenomenon: complex motion, unusual geometry
Part I The phenomenon: complex motion, unusual geometry Chapter 1 Chaotic motion 1.1 What is chaos? Certain long-lasting, sustained motion repeats itself exactly, periodically. Examples from everyday life
More informationRotational motion problems
Rotational motion problems. (Massive pulley) Masses m and m 2 are connected by a string that runs over a pulley of radius R and moment of inertia I. Find the acceleration of the two masses, as well as
More informationStepping Motors. Chapter 11 L E L F L D
Chapter 11 Stepping Motors In the synchronous motor, the combination of sinusoidally distributed windings and sinusoidally time varying current produces a smoothly rotating magnetic field. We can eliminate
More informationHB Coupled Pendulums Lab Coupled Pendulums
HB 04-19-00 Coupled Pendulums Lab 1 1 Coupled Pendulums Equipment Rotary Motion sensors mounted on a horizontal rod, vertical rods to hold horizontal rod, bench clamps to hold the vertical rods, rod clamps
More informationPHY 123 Lab 9 Simple Harmonic Motion
PHY 123 Lab 9 Simple Harmonic Motion (updated 11/17/16) The purpose of this lab is to study simple harmonic motion of a system consisting of a mass attached to a spring. You will establish the relationship
More informationCHAPTER 7: OSCILLATORY MOTION REQUIRES A SET OF CONDITIONS
CHAPTER 7: OSCILLATORY MOTION REQUIRES A SET OF CONDITIONS 7.1 Period and Frequency Anything that vibrates or repeats its motion regularly is said to have oscillatory motion (sometimes called harmonic
More informationME scope Application Note 28
App Note 8 www.vibetech.com 3/7/17 ME scope Application Note 8 Mathematics of a Mass-Spring-Damper System INTRODUCTION In this note, the capabilities of ME scope will be used to build a model of the mass-spring-damper
More informationHuman Arm. 1 Purpose. 2 Theory. 2.1 Equation of Motion for a Rotating Rigid Body
Human Arm Equipment: Capstone, Human Arm Model, 45 cm rod, sensor mounting clamp, sensor mounting studs, 2 cord locks, non elastic cord, elastic cord, two blue pasport force sensors, large table clamps,
More informationLecture V : Oscillatory motion and spectral analysis
Lecture V : Oscillatory motion and spectral analysis I. IDEAL PENDULUM AND STABILITY ANALYSIS Let us remind ourselves of the equation of motion for the pendulum. Remembering that the external torque applied
More informationLinear Second-Order Differential Equations LINEAR SECOND-ORDER DIFFERENTIAL EQUATIONS
11.11 LINEAR SECOND-ORDER DIFFERENTIAL EQUATIONS A Spring with Friction: Damped Oscillations The differential equation, which we used to describe the motion of a spring, disregards friction. But there
More informationLectures Chapter 10 (Cutnell & Johnson, Physics 7 th edition)
PH 201-4A spring 2007 Simple Harmonic Motion Lectures 24-25 Chapter 10 (Cutnell & Johnson, Physics 7 th edition) 1 The Ideal Spring Springs are objects that exhibit elastic behavior. It will return back
More informationPhysics 1050 Experiment 6. Moment of Inertia
Physics 1050 Moment of Inertia Prelab uestions These questions need to be completed before entering the lab. Please show all workings. Prelab 1 Sketch a graph of torque vs angular acceleration. Normal
More informationThe University of Hong Kong Department of Physics. Physics Laboratory PHYS3350 Classical Mechanics Experiment No The Physical Pendulum Name:
The University of Hong Kong Department of Physics Physics Laboratory PHYS3350 Classical Mechanics Experiment No. 3350-2 The Physical Pendulum Name: University No: Introduction One of the practical uses
More informationPhysics C: Mechanics
Physics C: Mechanics 2013 2014 PISCATAWAY TOWNSHIP SCHOOLS COURSE SYLLABUS Mr. Rohan Gokhale rgokhale@pway.org www.piscatawayschools.org/phs Brief Description of Course The AP Physics course is a full
More informationInvestigating the Relationship Between Cavendish Temperature Fluctuation and Torsional Oscillation
Investigating the Relationship Between Cavendish Temperature Fluctuation and Torsional Oscillation John Grasel 4 March 2010 Abstract The Cavendish apparatus measures the gravitational attraction between
More informationCircular Motion and Centripetal Force
[For International Campus Lab ONLY] Objective Measure the centripetal force with the radius, mass, and speed of a particle in uniform circular motion. Theory ----------------------------- Reference --------------------------
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department. Experiment 03: Work and Energy
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department Physics 8.01 Fall Term 2010 Experiment 03: Work and Energy Purpose of the Experiment: In this experiment you allow a cart to roll down an inclined
More informationPHYSICS 211 LAB #8: Periodic Motion
PHYSICS 211 LAB #8: Periodic Motion A Lab Consisting of 6 Activities Name: Section: TA: Date: Lab Partners: Circle the name of the person to whose report your group printouts will be attached. Individual
More informationRotational Motion. Variable Translational Motion Rotational Motion Position x θ Velocity v dx/dt ω dθ/dt Acceleration a dv/dt α dω/dt
Team: Rotational Motion Rotational motion is everywhere. When you push a door, it rotates. When you pedal a bike, the wheel rotates. When you start an engine, many parts rotate. Electrons rotate in an
More informationt = g = 10 m/s 2 = 2 s T = 2π g
Annotated Answers to the 1984 AP Physics C Mechanics Multiple Choice 1. D. Torque is the rotational analogue of force; F net = ma corresponds to τ net = Iα. 2. C. The horizontal speed does not affect the
More information