Class 4: More Pendulum results
|
|
- Maude Warner
- 6 years ago
- Views:
Transcription
1 Class 4: More Pendulum results The pendulum is a mechanical problem with a long and interesting history, from Galileo s first ansatz that the period was independent of the amplitude based on watching priests swinging incense burners, to the modern study of chaotic systems, the pendulum has been there. Pendulum clocks became the very model of deterministic science after Newton, when analogies were made to a clockwork universe that obeyed the laws of physics and was completely determined by its initial conditions. From a paradigm of determinism to a classic example of unpredictable chaos, the pendulum is an ideal system to learn some basic dynamics. We calculated the formal analytic solution last time and found that, although useful for proving various properties of pendulum dynamics, we need to the numerical solution to build physical intuition for the system. Pendulum orbits: Numerical solution Time series of angle and angular velocity We plotted the analytic solution last time and, at least in one situation, we got gaps in the plot, probably due to insufficient sampling (or possibly no-unique behavior of the Elliptic functions in some regions). Here's a numerical solution of the pendulum equation of motion θθ + ωω 0 sinθθ = 0 with ωω 0 =. θ0 =.5; s = NDSolve[{θ''[t] + Sin[θ[t]] == 0, θ[0] θ0, θ'[0] == 0}, θ, {t, 0, 30}]; Plot {Evaluate[{θ[t]} /. s], θ0 Cos[t]}, {t, 0, 30}, PlotStyle {Thick, Automatic}, Frame True, FrameLabel {"t", "θ"}, RotateLabel False, AspectRatio GoldenRatio θθ t The blue curve is the actual pendulum - the Jacobi elliptic function - and the orange curve is the cosine (harmonic oscillator) oscillation. You can see that the shape of the curves are similar but that the frequencies are different. In fact, nonlinear oscillators like the pendulum usually have amplitudedependent frequencies or frequency-dependent amplitudes, depending on how you look at it. This
2 dependent frequencies or frequency-dependent amplitudes, depending on how you look at it. This plot is for an initial displacement angle of θθ 0 =.5 radians, almost 90 degrees. If we increase the initial angle to 3 radians (about 7 degrees) we get: θ0 = 3.0; s = NDSolve[{θ''[t] + Sin[θ[t]] == 0, θ[0] θ0, θ'[0] == 0}, θ, {t, 0, 50}]; Plot {Evaluate[{θ[t]} /. s], θ0 Cos[t]}, {t, 0, 50}, PlotStyle {Thick, Automatic}, Frame True, FrameLabel {"t", "θ"}, RotateLabel False, AspectRatio GoldenRatio 3 θθ t This is the case that gave problems plotting from the analytic solution. No problem numerically. For this larger amplitude the harmonic oscillator curve has the same frequency it did for the smaller amplitude, but the pendulum frequency is now much lower. The shape of the pendulum curve even looks different than the cosine function -- it's curvature at the maxima no longer looks sinusoidal, but is more rounded. You're seeing the difference between a harmonic (linear) oscillation and one kind of nonlinear oscillation. Internalize it. Before moving on, it's worth looking at the pendulum motion for a small initial angle, say 0.07 radians, or about. θ0 = 0.07; s = NDSolve[{θ''[t] + Sin[θ[t]] == 0, θ[0] θ0, θ'[0] == 0}, θ, {t, 0, 30}]; Plot {Evaluate[{θ[t]} /. s], θ0 Cos[t]}, {t, 0, 30}, PlotStyle {{Thick, Dashed}, Automatic}, Frame True, FrameLabel {"t", "θ"}, RotateLabel False, AspectRatio GoldenRatio For this small angle, the pendulum motion is virtually identical to the cosine function (yes, there are curves plotted here, they just overlap and look like one). You can easily verify this theoretically, by solving the equation of motion (4) for small angles (use the first term in a Taylor series for the sine function) - you'll find the small angle oscillation is sinusoidal with fixed frequency ωω 0. Producing a Phase Portrait for the pendulum Phase space plots, like the one we obtained last class, show us how the dynamical variables are related. But, for the pendulum, we already know how they are related: by the energy integral! We don't need a numerical solution to get a phase plot, just use the conservation of energy equation : m L θθ + m ωω0 L cosθθ = E
3 3 Solving this equation for θθ + gives a relation that represents a curve in θθ-θθ + space. Each curve will be associated with a different value of the total energy E. Before we plot this curve, let's try to use our physical intuition to figure out what we expect the result to look like. For example, for small angles, θθ <<, we know the pendulum motion will be a simple oscillation - i.e. swinging back and forth - also called libration. If we start the pendulum with zero angular velocity and initial angle θθ 0 we expect it to traverse a curve not unlike the phase space orbit above: it will move downward and to the left until it gets to θθ = 0, where its angular velocity will reach a maximum. After that it will swing over to negative angles until it reaches its maximum angle θθ 0 with zero θθ +, then returning to θθ = 0 with a positive θθ +, and finally reaching its initial position with zero angular velocity - only to repeat the whole operation ad infinitum. On a θθ-θθ + graph, this would look like a closed curve. Let's let Mathematica plot it exactly for us: define En = E / m L so, and let's plot for a sample frequency ωω 0 =, so the energy integral becomes θθ + cosθθ = E m L = En, or, converting to Mathematica language: θdot = Cos[θ] + En, Cos[θ] + En ; As noted above, small oscillations will have negative energies, so let's choose En = 0.5 for an example: c = Plot θdot /. En.5, {θ, π /, π / }, PlotPoints 00, Frame True, FrameLabel "θ", "θ ", RotateLabel False, AspectRatio GoldenRatio θθ θθ This nearly circular phase space orbit for E / m L = 0.5 has an amplitude of (can you prove this?). Now, imagine a different initial condition, say starting the pendulum bob at θθ 0 = ππ, with positive initial velocity θθ + 0 > 0. What happens? A little thought should tell you that the pendulum would now rotate, since it has kinetic energy at the top of its swing so when it reaches that maximum height again, it retains that energy and will keep moving. What does this rotation look like on a phase-space plot? It would start at θθ = ππ then increase upward and leftward until it reaches θθ = 0, again with a maximum angular velocity. Then it would swing into the negative angle range until it gets to θθ = ππ again, where it would have angular velocity equal to its initial velocity - its minimum angular speed - and then it would keep going, repeating this pattern for larger and larger negative angles. In Mathematica, it would look like this:
4 4 c = Plot θdot /. En 3, {θ, 5 π, π}, PlotPoints 00, PlotStyle {Purple}, Frame True, FrameLabel "θ", "θ ", RotateLabel False, AspectRatio GoldenRatio 3 θθ θθ The upper curve is the one we just described, starting with a positive angular velocity. The lower curve is the corresponding trajectory for a negative initial angular velocity at the same energy. Combining these plots for two different initial conditions on the same axes gives: Show[c, c, PlotRange All] 3 θθ θθ So, we can easily see the difference between libration (swinging) and rotation by looking at the graph. In fact, we can guess that there should be a boundary curve that divides regions of libration from regions of rotation: a curve that's usually called the separatrix. Here it is for the example case we're plotting:
5 5 c3 = Plot θdot /. En, {θ, 5 π, π}, PlotPoints 00, PlotStyle {Red}, Frame True, FrameLabel "θ", "θ ", RotateLabel False, AspectRatio GoldenRatio ; Show[c, c, c3, PlotRange All] 3 θθ θθ What does the separatrix curve represent physically for the pendulum? It's an idealized orbit that can't really be achieved with a real pendulum: for example, start it at the top of its swing, θθ = ππ, with zero velocity - if it could then somehow be infinitesimally perturbed (with no change in energy) so it would swing, it would swing around and come back to the top and stop (and it would take an infinite time to do it)! That's the separatrix orbit. These three orbit types are characteristic of the pendulum - all orbits of a plane pendulum are of one of these types. Let's plot some other orbits (i.e. different energies) to get a global picture of pendulum behavior.
6 6 θdotp = θdot /.En.5; θdotp = θdot /.En.; θdotp3 = θdot /.En ; θdotp4 = θdot /.En ; θdotp5 = θdot /.En 3; θdotp6 = θdot /.En 4; Plot {θdotp, θdotp, θdotp3, θdotp4, θdotp5, θdotp6}, {θ, 3 π, 3 π}, Frame True, PlotPoints 00, FrameLabel "θ", "θ ", RotateLabel False, PlotLabel "Pendulum Phase Portrait" 3 Pendulum Phase Portrait θθ θθ This plot is called a phase portrait for the pendulum. Does this agree with your calculation of the turning points in Homework? From the phase portrait you can determine all the possible types of behavior of this system with two variables (θθ and θθ + ). Phase portraits can be an extremely useful tool for initial investigations of a new nonlinear system. Note that the phase portrait is really just a way of visualizing the energy integral (each curve is a solution to the energy integral with a different constant energy), so we see that the energy integral really contains all the information we need to understand the motion! Using a contour plot to visualize the Phase Portrait for the pendulum Here's another, probably easier, way of viewing the phase portrait: the energy integral can be thought of as a two-dimensional function, energy E θθ, θθ + as a function of the two dynamical variables, θθ and θθ +. If you plot a contour map of the energy function, the contours will be at constant values of energy and therefore will represent the phase space orbits at each energy. Rewriting the energy integral:
7 7 energy = θdot Cos[θ]; ContourPlot[energy, {θ, 0, 0}, {θdot, 3.5, 3.5}, ImageSize Small] This default plot can be made more useful with options. For example, let's get rid of the color shading between contours. Also, one advantage of the contour plot is that we can label each orbit (contour line) with its energy: ContourPlot energy, {θ, 0, 0}, {θdot, 3.5, 3.5}, ContourShading False, ContourLabels All, AspectRatio Automatic, FrameLabel "θ", "θ ", RotateLabel False, PlotLabel "Pendulum Phase Portrait" If you don't like the default contour energy values, you can select your own. If you want tweak further, you can specify the significant figures and number of decimal places for the labels. Let s also clean up the plot by eliminating those gaps at multiples of θθ = ππ, and by labeling the θθ-axis in multiples of ππ:
8 8 θmin = 0; θmax = 0; tickmin = Ceiling θmin π ; tickmax = Floor[θmax / π]; pendpp = ContourPlot energy, {θ, θmin, θmax}, {θdot, 3.5, 3.5}, PlotPoints 50, ContourShading False, Contours {0.53, 0.0,.0,.5, 5.0}, ContourStyle Blue, ContourLabels (Text[NumberForm[#3, {, }], {#, #}] &), FrameTicks {Table[tickmin π + k π, {k, 0, tickmax tickmin}], Automatic}, AspectRatio Automatic, FrameLabel "θ", "θ ", RotateLabel False, PlotLabel "Pendulum Phase Portrait", ImageSize Large The NumberForm command tells Mathematica to label variable 3 (energy) with significant figures and decimal places. I also added more sampling points with the PlotPoints command to clear up that "gap" in the separatrix curve; the extra sampling looks nicer, but does take more compute time. Finally, I changed the θθ axis to be in multiples of ππ radians using the Ticks command. Compare the pendulum phase portrait to that of the simple harmonic oscillator I mentioned above that if you approximate the sine function near θθ = 0 with the first term in its Taylor series, that you end up with the equation of motion for a simple harmonic oscillator: θθ + ωω 0 θθ 0 which (as you no doubt learned in PHY 0...) has a simple sinusoidal solution θθ(t) = A sinωω 0 t + B cosωω 0 t. It's not hard to show that the energy integral for this equation is θθ+ + ωω 0 θθ = E m L Let's plot the harmonic oscillator phase portrait using this energy equation. Calling energy = E / m L and taking ωω 0 = :
9 9 energysho = θdot + θ ; shopp = ContourPlot energysho, {θ, 4, 4}, {θdot, 3.5, 3.5}, PlotPoints 00, ContourShading False, Contours {0.0,.0,.5, 5.0, 7.5}, ContourStyle Red, AspectRatio Automatic, FrameLabel "θ", "θ ", RotateLabel False, PlotLabel "Simple Harmonic Oscillator Phase Portrait" ; Show[pendPP, shopp, PlotRange {{4, 4}, Automatic}, PlotLabel "Pendulum (blue), Harmonic oscillator (red)"] You can see that the oscillator phase portrait is just a set of concentric circles representing oscillations: no rotational motion and no separatrix. Thus it is only a good approximation for very small angles, that is it's valid locally, near the origin (0, 0) of the phase plot. Using the small angle Taylor expansion to first order is sometimes called linearizing the equation of motion locally. The global behavior, shown in the full pendulum phase plot is quite different. As we shall see, we can often learn some important information by linearizing locally, but we need to be careful that we don't extrapolate the local behavior to the global system, which may be totally different.
10 0 Higher order approximation to the pendulum equation You looked at the third order approximation to the pendulum equation in Class. Now lets look at higher order approximations. 3th order approximation Here s an example of the 3 th order approximation to the pendulum force (keeping 3 terms in the Taylor series): order3 = Normal[Series[Sin[θ[t]], {θ[t], 0, 3}]] EOM = θ''[t] + ω0 (order3) En = Integrate[EOM θ'[t], t] θ[t] θ[t]3 6 + θ[t]5 0 θ[t] θ[t]9 θ[t] θ[t] ω0 θ[t] θ[t]3 6 + θ[t]5 0 ω0 θ[t] 4 ω0 θ[t] ω0 θ[t] 6 ω0 θ[t] ω0 θ[t] θ[t] θ[t]9 θ[t] θ[t] θ [t] ω0 θ[t] ω0 θ[t] θ [t] Prepare to plot by getting separatrix locations. One way to do that is to locate the fixed points of the motion: where the acceleration and the velocity vanish. The acceleration here is proportional to the series for sinθθ, i.e. what I called oder3 above. Let s set that equal zero and find the real solutions only: Solve[order3 0, θ[t], Reals]; N[%] {{θ[t] 0.}, {θ[t] }, {θ[t] 3.46}, {θ[t] 3.46}, {θ[t] }} We now have positive real roots (and two negative). Get the energy of the first real root: En = En /. ω0 /. {θ[t] θ, θ'[t] ω} Es = En /. {ω 0, θ 3.46} θ. θ4 4 + θ6 θ θ0 θ θ ω Get the energy of the second real root: Es = En /. {ω 0, θ } Here are the plots:
11 approx3 = ContourPlot En, {θ, 8, 8}, {ω, 3, 3}, PlotPoints 50, ContourShading False, Contours {0.0, Es, 0.50,.00,.50, Es,.50, 3.0, 3.5, 4.0, 4.5, 5.0}, ContourStyle Blue, FrameLabel θ, θ, RotateLabel False, AspectRatio Automatic, ImageSize Large ; Epend = ω + ( Cos[θ]); pend = ContourPlot Epend, {θ, 8, 8}, {ω, 3, 3}, PlotPoints 50, ContourShading False, Contours {0.0, Es, 0.50,.00,.50, Es,.50, 3.0, 3.5, 4.0, 4.5, 5.0}, ContourLabels True, ContourStyle Red, FrameLabel θ, θ, RotateLabel False, AspectRatio Automatic, ImageSize Large ; Show[approx3, pend, PlotRange All, PlotLabel "Full pendulum (red) and 3th order (blue)"] We get pretty good agreement out to an angle of about 5 radians, but the separatrix structure is still inaccurate. As mentioned in class, note that for every other order (5 th, 9 th, 3 th, etc.) the unphysical orbits that blow up are not present. 5th order approximation Here s an example of the 5 th order approximation (keeping 5 terms in the Taylor series): order5 = Normal[Series[Sin[θ[t]], {θ[t], 0, 5}]] EOM = θ''[t] + ω0 (order5); En = Integrate[EOM θ'[t], t]; θ[t] θ[t]3 6 + θ[t]5 0 θ[t] θ[t]9 θ[t] θ[t]3 θ[t] Prepare to plot by getting fixed point locations (real solutions only) & real root energies:
12 Solve[order5 0, θ[t], Reals]; N[%] En = En /. ω0 /. {θ[t] θ, θ'[t] ω}; Es = En /. {ω 0, θ 3.459} Es = En /. {ω 0, θ } {{θ[t] 0.}, {θ[t] 7.057}, {θ[t] }, {θ[t] 3.459}, {θ[t] 3.459}, {θ[t] }, {θ[t] 7.057}} Here are the plots: approx5 = ContourPlot En, {θ, 8, 8}, {ω, 3, 3}, PlotPoints 50, ContourShading False, Contours {Es, 0.0, 0.06, 0.50,.00,.50, Es,.50, 3.0, 3.5, 4.0, 4.5, 5.0}, ContourStyle Blue, FrameLabel θ, θ, RotateLabel False, AspectRatio Automatic, ImageSize Large ; pend = ContourPlot Epend, {θ, 8, 8}, {ω, 3, 3}, PlotPoints 50, ContourShading False, Contours {Es, 0.0, 0.06, 0.50,.00,.50, Es,.50, 3.0, 3.5, 4.0, 4.5, 5.0}, ContourLabels True, ContourStyle Red, FrameLabel θ, θ, RotateLabel False, AspectRatio Automatic, ImageSize Large, PlotLabel "Full Pendulum: phase portrait" ; Show[approx5, pend, PlotRange All, PlotLabel "Full pendulum (red) and 5th order approximation (blue)"] For this approximation the infinite orbits return at right and left. We get pretty good agreement out to an angle of a little over 5 radians, and now there are two unstable positive fixed points, but the second separatrix is not connected to the first! 7th order approximation Here s an example of the 7 th order approximation (keeping 7 terms in the Taylor series):c
13 3 order7 = Normal[Series[Sin[θ[t]], {θ[t], 0, 7}]] EOM = θ''[t] + ω0 (order7); En = Integrate[EOM θ'[t], t]; Solve[order7 0, θ[t], Reals]; N[%] En = En /. ω0 /. {θ[t] θ, θ'[t] ω}; Es = En /. {ω 0, θ 3.459}; Es = En /. {ω 0, θ 6.839}; Es3 = En /. {ω 0, θ 9.494}; θ[t] θ[t]3 + θ[t]5 6 0 θ[t] θ[t] θ[t] θ[t] θ[t] θ[t]9 θ[t] θ[t] θ[t] θ[t] θ[t] {{θ[t] 0.}, {θ[t] 9.494}, {θ[t] 6.839}, {θ[t] 3.459}, {θ[t] 3.459}, {θ[t] 6.839}, {θ[t] 9.494}} approx7 = ContourPlot En, {θ, 3, 3}, {ω, 3, 3}, PlotPoints 50, ContourShading False, Contours {Es, Es, Es3, 0.50,.00,.5,.50, 3.0, 3.5, 4.0, 4.5, 5.0}, ContourStyle Blue, FrameLabel θ, θ, RotateLabel False, AspectRatio Automatic, ImageSize Large ; pend = ContourPlot Epend, {θ, 3, 3}, {ω, 3, 3}, PlotPoints 50, ContourShading False, Contours {Es, Es, Es3, 0.50,.00,.5,.50, 3.0, 3.5, 4.0, 4.5, 5.0}, ContourLabels True, ContourStyle Red, FrameLabel θ, θ, RotateLabel False, AspectRatio Automatic, ImageSize Large ; Show[approx7, pend, PlotRange All, PlotLabel "Full pendulum (red) and 7th order approximation (blue)"] OK, things are looking up! We get good agreement out to an angle of about 0 radians, including the separatrix with multiple fixed points, after which the divergent orbits show up. Conclusion: expanding the force in a Taylor series really does get more and more accurate as the order of the approximation increases -- even though the actual terms are minuscule (e.g. 7 order term has coefficient 0 8 ).
14 4 Fixed points of the pendulum equation The phase portrait for the pendulum has some obvious points of interest: the points where the separatrix appears to intersect itself, for example (θθ = ( n + )ππ, i.e. odd multiples of ππ) and the centers of the oscillatory orbit regions (θθ = nππ, i.e. even multiples of ππ). A few moments thought and you ll realize that if you start the pendulum at any of those points, nothing happens! They re called fixed points or equilibrium points, and are defined by zero acceleration and velocity. From the equation of motion θθ + ωω 0 sinθθ = 0 with ωω 0 =, we see that zero acceleration implies sinθθ = 0. So the fixed points occur at all angles θθ * satisfying θθ * = sin 0, or fixed point position: θθ * = n ππ, for any integer n We can determine what the motion is like near the fixed point by linearizing about the fixed point, i.e. expanding the force to first order, dropping higher order terms: F(θθ) = F(θθ * ) + F θθ * (θθ θθ* ) = F θθ * (θθ θθ* ) where the second equality is true because the Force must vanish at any fixed point. The next step is to guess an exponential solution θθ θθ * = A e αα t and solve for the constant αα: any positive real value implies the motion leaves the fixed point (motion is unstable), and imaginary αα implies oscillation about the fixed point (a type of stable motion). Let s try the pendulum fixed points: F = sinθθ (for ωω 0 = ) so the linearized equation of motion is θθ = (cosθθ) * (θθ θθ * ). Fixed point θθ * = 0: In this case, the linearized equation becomes θθ = θθ and trial solution θθ = A e αα t gives αα = so αα = ±i, which implies oscillatory motion about the fixed point, a type of stable motion. Fixed point θθ * = ππ: In this case, the linearized equation becomes θθ = +θθ and trial solution θθ = A e αα t gives αα = + so αα = ±, which implies exponential motion, with the positive exponent giving unstable motion. It should be clear that all the fixed points at multiples of ππ will be the same as the zero fixed point, and the same for thee odd multiples of ππ: they ll be the same as the θθ * = ππ fixed point. Hence: Pendulum fixed points: Fixed point θ * = n π θ * = ( n + ) π Stability Stable Unstable where n is an integer. This agrees with the phase portrait above: there are closed loop orbits surrounding the stable fixed points, implying oscillation, while most orbits leave the unstable points, implying instability. Note that you will have both oscillatory and rotating orbits in the neighborhood of the unstable fixed points.
15 5
Nonlinear Oscillators: Free Response
20 Nonlinear Oscillators: Free Response Tools Used in Lab 20 Pendulums To the Instructor: This lab is just an introduction to the nonlinear phase portraits, but the connection between phase portraits and
More informationContents. Contents. Contents
Physics 121 for Majors Class 18 Linear Harmonic Last Class We saw how motion in a circle is mathematically similar to motion in a straight line. We learned that there is a centripetal acceleration (and
More informationAP Physics C Mechanics
1 AP Physics C Mechanics Simple Harmonic Motion 2015 12 05 www.njctl.org 2 Table of Contents Click on the topic to go to that section Spring and a Block Energy of SHM SHM and UCM Simple and Physical Pendulums
More informationOSCILLATIONS ABOUT EQUILIBRIUM
OSCILLATIONS ABOUT EQUILIBRIUM Chapter 13 Units of Chapter 13 Periodic Motion Simple Harmonic Motion Connections between Uniform Circular Motion and Simple Harmonic Motion The Period of a Mass on a Spring
More informationMass on a Horizontal Spring
Course- B.Sc. Applied Physical Science (Computer Science) Year- IInd, Sem- IVth Subject Physics Paper- XIVth, Electromagnetic Theory Lecture No. 22, Simple Harmonic Motion Introduction Hello friends in
More informationChapter 11 Vibrations and Waves
Chapter 11 Vibrations and Waves If an object vibrates or oscillates back and forth over the same path, each cycle taking the same amount of time, the motion is called periodic. The mass and spring system
More informationMITOCW 6. Standing Waves Part I
MITOCW 6. Standing Waves Part I The following content is provided under a Creative Commons license. Your support will help MIT OpenCourseWare continue to offer high quality educational resources for free.
More informationLab 10: Harmonic Motion and the Pendulum
Lab 10 Harmonic Motion and the Pendulum 119 Name Date Partners Lab 10: Harmonic Motion and the Pendulum OVERVIEW A body is said to be in a position of stable equilibrium if, after displacement in any direction,
More informationChapter 13 Oscillations about Equilibrium. Copyright 2010 Pearson Education, Inc.
Chapter 13 Oscillations about Equilibrium Periodic Motion Units of Chapter 13 Simple Harmonic Motion Connections between Uniform Circular Motion and Simple Harmonic Motion The Period of a Mass on a Spring
More informationLecture 1 Notes: 06 / 27. The first part of this class will primarily cover oscillating systems (harmonic oscillators and waves).
Lecture 1 Notes: 06 / 27 The first part of this class will primarily cover oscillating systems (harmonic oscillators and waves). These systems are very common in nature - a system displaced from equilibrium
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 informationEnergy in a Simple Harmonic Oscillator. Class 30. Simple Harmonic Motion
Simple Harmonic Motion Class 30 Here is a simulation of a mass hanging from a spring. This is a case of stable equilibrium in which there is a large extension in which the restoring force is linear in
More informationVibratory Motion -- Conceptual Solutions
Vibratory Motion Vibratory Motion -- Conceptual Solutions 1.) An ideal spring attached to a mass m =.3 kg provides a force equal to -kx, where k = 47.33 nt/m is the spring's spring constant and x denotes
More informationThe Nonlinear Pendulum
The Nonlinear Pendulum - Pádraig Ó Conbhuí - 08531749 TP Monday 1. Abstract This experiment was performed to examine the effects that linearizing equations has on the accuracy of results and to find ways
More informationLAB 10: HARMONIC MOTION AND THE PENDULUM
163 Name Date Partners LAB 10: HARMONIC MOION AND HE PENDULUM Galileo reportedly began his study of the pendulum in 1581 while watching this chandelier swing in Pisa, Italy OVERVIEW A body is said to be
More informationWhy are Discrete Maps Sufficient?
Why are Discrete Maps Sufficient? Why do dynamical systems specialists study maps of the form x n+ 1 = f ( xn), (time is discrete) when much of the world around us evolves continuously, and is thus well
More informationMATH 415, WEEKS 7 & 8: Conservative and Hamiltonian Systems, Non-linear Pendulum
MATH 415, WEEKS 7 & 8: Conservative and Hamiltonian Systems, Non-linear Pendulum Reconsider the following example from last week: dx dt = x y dy dt = x2 y. We were able to determine many qualitative features
More informationDepending on your Mathematica installation, you might see the solution curve corresponding to q
Problem 5 (a) Clear@pendD; pend@t_, x0_, y0_d := x@td, y@td
More informationFor a rigid body that is constrained to rotate about a fixed axis, the gravitational torque about the axis is
Experiment 14 The Physical Pendulum The period of oscillation of a physical pendulum is found to a high degree of accuracy by two methods: theory and experiment. The values are then compared. Theory For
More informationAP Physics 1. April 11, Simple Harmonic Motion. Table of Contents. Period. SHM and Circular Motion
AP Physics 1 2016-07-20 www.njctl.org Table of Contents Click on the topic to go to that section Period and Frequency SHM and UCM Spring Pendulum Simple Pendulum Sinusoidal Nature of SHM Period and Frequency
More informationMITOCW MITRES18_005S10_DiffEqnsMotion_300k_512kb-mp4
MITOCW MITRES18_005S10_DiffEqnsMotion_300k_512kb-mp4 PROFESSOR: OK, this lecture, this day, is differential equations day. I just feel even though these are not on the BC exams, that we've got everything
More informationChapter 4. Oscillatory Motion. 4.1 The Important Stuff Simple Harmonic Motion
Chapter 4 Oscillatory Motion 4.1 The Important Stuff 4.1.1 Simple Harmonic Motion In this chapter we consider systems which have a motion which repeats itself in time, that is, it is periodic. In particular
More informationResonance and response
Chapter 2 Resonance and response Last updated September 20, 2008 In this section of the course we begin with a very simple system a mass hanging from a spring and see how some remarkable ideas emerge.
More information!T = 2# T = 2! " The velocity and acceleration of the object are found by taking the first and second derivative of the position:
A pendulum swinging back and forth or a mass oscillating on a spring are two examples of (SHM.) SHM occurs any time the position of an object as a function of time can be represented by a sine wave. We
More informationLab 11. Spring-Mass Oscillations
Lab 11. Spring-Mass Oscillations Goals To determine experimentally whether the supplied spring obeys Hooke s law, and if so, to calculate its spring constant. To find a solution to the differential equation
More informationPhysics 235 Chapter 4. Chapter 4 Non-Linear Oscillations and Chaos
Chapter 4 Non-Linear Oscillations and Chaos Non-Linear Differential Equations Up to now we have considered differential equations with terms that are proportional to the acceleration, the velocity, and
More informationMaps and differential equations
Maps and differential equations Marc R. Roussel November 8, 2005 Maps are algebraic rules for computing the next state of dynamical systems in discrete time. Differential equations and maps have a number
More informationLAB 10 - HARMONIC MOTION AND THE PENDULUM
L10-1 Name Date Partners LAB 10 - HARMONIC MOION AND HE PENDULUM θ L Groove marking the center of mass Photogate s = 0 s F tan mg θ OVERVIEW Figure 1 A body is said to be in a position of stable equilibrium
More informationChapter 14 Oscillations. Copyright 2009 Pearson Education, Inc.
Chapter 14 Oscillations Oscillations of a Spring Simple Harmonic Motion Energy in the Simple Harmonic Oscillator Simple Harmonic Motion Related to Uniform Circular Motion The Simple Pendulum The Physical
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 informationDaba Meshesha Gusu and O.Chandra Sekhara Reddy 1
International Journal of Basic and Applied Sciences Vol. 4. No. 1 2015. Pp.22-27 Copyright by CRDEEP. All Rights Reserved. Full Length Research Paper Solutions of Non Linear Ordinary Differential Equations
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 informationPhysics 351, Spring 2015, Homework #5. Due at start of class, Friday, February 20, 2015 Course info is at positron.hep.upenn.
Physics 351, Spring 2015, Homework #5. Due at start of class, Friday, February 20, 2015 Course info is at positron.hep.upenn.edu/p351 When you finish this homework, remember to visit the feedback page
More informationLinear and Nonlinear Oscillators (Lecture 2)
Linear and Nonlinear Oscillators (Lecture 2) January 25, 2016 7/441 Lecture outline A simple model of a linear oscillator lies in the foundation of many physical phenomena in accelerator dynamics. A typical
More informationLab 12. Spring-Mass Oscillations
Lab 12. Spring-Mass Oscillations Goals To determine experimentally whether the supplied spring obeys Hooke s law, and if so, to calculate its spring constant. To determine the spring constant by another
More informationspring mass equilibrium position +v max
Lecture 20 Oscillations (Chapter 11) Review of Simple Harmonic Motion Parameters Graphical Representation of SHM Review of mass-spring pendulum periods Let s review Simple Harmonic Motion. Recall we used
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 informationChaos and Liapunov exponents
PHYS347 INTRODUCTION TO NONLINEAR PHYSICS - 2/22 Chaos and Liapunov exponents Definition of chaos In the lectures we followed Strogatz and defined chaos as aperiodic long-term behaviour in a deterministic
More informationChapter 11 Parametric Equations, Polar Curves, and Conic Sections
Chapter 11 Parametric Equations, Polar Curves, and Conic Sections ü 11.1 Parametric Equations Students should read Sections 11.1-11. of Rogawski's Calculus [1] for a detailed discussion of the material
More informationTheoretical physics. Deterministic chaos in classical physics. Martin Scholtz
Theoretical physics Deterministic chaos in classical physics Martin Scholtz scholtzzz@gmail.com Fundamental physical theories and role of classical mechanics. Intuitive characteristics of chaos. Newton
More informationChap. 15: Simple Harmonic Motion
Chap. 15: Simple Harmonic Motion Announcements: CAPA is due next Tuesday and next Friday. Web page: http://www.colorado.edu/physics/phys1110/phys1110_sp12/ Examples of periodic motion vibrating guitar
More informationSimple Harmonic Motion
1. Object Simple Harmonic Motion To determine the period of motion of objects that are executing simple harmonic motion and to check the theoretical prediction of such periods. 2. Apparatus Assorted weights
More informationP321(b), Assignement 1
P31(b), Assignement 1 1 Exercise 3.1 (Fetter and Walecka) a) The problem is that of a point mass rotating along a circle of radius a, rotating with a constant angular velocity Ω. Generally, 3 coordinates
More informationProblem Set 1 October 30, 2017
1. e π can be calculated from e x = x n. But that s not a good approximation method. The n! reason is that π is not small compared to 1. If you want O(0.1) accuracy, then the last term you need to include
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 informationCHAPTER 9 -- VIBRATORY MOTION QUESTION SOLUTIONS
Solutions--Ch. 9 (Vibratory Motion) CHAPTER 9 -- VIBRATORY MOTION QUESTION SOLUTIONS 9.1) An ideal spring attached to a mass m =.3 kg provides a force equal to -kx, where k = 47.33 nt/m is the spring's
More informationDifferential Equations
Electricity and Magnetism I (P331) M. R. Shepherd October 14, 2008 Differential Equations The purpose of this note is to provide some supplementary background on differential equations. The problems discussed
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 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 informationTopic 5 Notes Jeremy Orloff. 5 Homogeneous, linear, constant coefficient differential equations
Topic 5 Notes Jeremy Orloff 5 Homogeneous, linear, constant coefficient differential equations 5.1 Goals 1. Be able to solve homogeneous constant coefficient linear differential equations using the method
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 informationLecture 2: Series Expansion
Lecture 2: Series Expansion Key points Maclaurin series : For, Taylor expansion: Maple commands series convert taylor, Student[NumericalAnalysis][Taylor] 1 Maclaurin series of elementary functions for
More informationLab 4 Numerical simulation of a crane
Lab 4 Numerical simulation of a crane Agenda Time 10 min Item Review agenda Introduce the crane problem 95 min Lab activity I ll try to give you a 5- minute warning before the end of the lab period to
More informationGeneral Physics I Spring Oscillations
General Physics I Spring 2011 Oscillations 1 Oscillations A quantity is said to exhibit oscillations if it varies with time about an equilibrium or reference value in a repetitive fashion. Oscillations
More informationPhysics 342 Lecture 23. Radial Separation. Lecture 23. Physics 342 Quantum Mechanics I
Physics 342 Lecture 23 Radial Separation Lecture 23 Physics 342 Quantum Mechanics I Friday, March 26th, 2010 We begin our spherical solutions with the simplest possible case zero potential. Aside from
More informationThe dynamics of the Forced Damped Pendulum. John Hubbard Cornell University and Université de Provence
The dynamics of the Forced Damped Pendulum John Hubbard Cornell University and Université de Provence Three ways to view the pendulum Three ways to view the pendulum 1. As a physical object Three ways
More informationSimple harmonic motion the motion of springs is a very important topic in physics.
Chapter 11 Potential and Kinetic Energy Together: Simple Harmonic Motion In This Chapter Using Hooke s law Working with simple harmonic motion Calculating simple harmonic motion velcoity Finding simple
More informationQM and Angular Momentum
Chapter 5 QM and Angular Momentum 5. Angular Momentum Operators In your Introductory Quantum Mechanics (QM) course you learned about the basic properties of low spin systems. Here we want to review that
More informationA plane autonomous system is a pair of simultaneous first-order differential equations,
Chapter 11 Phase-Plane Techniques 11.1 Plane Autonomous Systems A plane autonomous system is a pair of simultaneous first-order differential equations, ẋ = f(x, y), ẏ = g(x, y). This system has an equilibrium
More informationLecture 15 - Orbit Problems
Lecture 15 - Orbit Problems A Puzzle... The ellipse shown below has one focus at the origin and its major axis lies along the x-axis. The ellipse has a semimajor axis of length a and a semi-minor axis
More information28. Pendulum phase portrait Draw the phase portrait for the pendulum (supported by an inextensible rod)
28. Pendulum phase portrait Draw the phase portrait for the pendulum (supported by an inextensible rod) θ + ω 2 sin θ = 0. Indicate the stable equilibrium points as well as the unstable equilibrium points.
More informationHow To Pump A Swing. Tareq Ahmed Mokhiemer Research Assistant, Physics Department.
How To Pump A Swing Tareq Ahmed Mokhiemer Research Assistant, Physics Department Abstract The analysis of pumping a swing has been done in the standing mode and seated mode in different ways using Mathematica
More informationCalculus I Homework: The Derivatives of Polynomials and Exponential Functions Page 1
Calculus I Homework: The Derivatives of Polynomials and Exponential Functions Page 1 Questions Example Differentiate the function y = ae v + b v + c v 2. Example Differentiate the function y = A + B x
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 informationThe object of this experiment is to study systems undergoing simple harmonic motion.
Chapter 9 Simple Harmonic Motion 9.1 Purpose The object of this experiment is to study systems undergoing simple harmonic motion. 9.2 Introduction This experiment will develop your ability to perform calculations
More informationPhysics 8 Monday, December 4, 2017
Physics 8 Monday, December 4, 2017 HW12 due Friday. Grace will do a review session Dec 12 or 13. When? I will do a review session: afternoon Dec 17? Evening Dec 18? Wednesday, I will hand out the practice
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 informationFebruary 13, Option 9 Overview. Mind Map
Option 9 Overview Mind Map Return tests - will discuss Wed..1.1 J.1: #1def,2,3,6,7 (Sequences) 1. Develop and understand basic ideas about sequences. J.2: #1,3,4,6 (Monotonic convergence) A quick review:
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 informationChapter 13 Lecture. Essential University Physics Richard Wolfson 2 nd Edition. Oscillatory Motion Pearson Education, Inc.
Chapter 13 Lecture Essential University Physics Richard Wolfson nd Edition Oscillatory Motion Slide 13-1 In this lecture you ll learn To describe the conditions under which oscillatory motion occurs To
More informationMITOCW ocw f99-lec23_300k
MITOCW ocw-18.06-f99-lec23_300k -- and lift-off on differential equations. So, this section is about how to solve a system of first order, first derivative, constant coefficient linear equations. And if
More informationFigure 1: Doing work on a block by pushing it across the floor.
Work Let s imagine I have a block which I m pushing across the floor, shown in Figure 1. If I m moving the block at constant velocity, then I know that I have to apply a force to compensate the effects
More informationMITOCW 5. Traveling Waves without Damping
MITOCW 5. Traveling Waves without Damping The following content is provided under a Creative Commons license. Your support will help MIT OpenCourseWare continue to offer high-quality educational resources
More informationComputational Physics (6810): Session 8
Computational Physics (6810): Session 8 Dick Furnstahl Nuclear Theory Group OSU Physics Department February 24, 2014 Differential equation solving Session 7 Preview Session 8 Stuff Solving differential
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 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 informationDifferentiation of Parametric Space Curves. Goals: Velocity in parametric curves Acceleration in parametric curves
Block #2: Differentiation of Parametric Space Curves Goals: Velocity in parametric curves Acceleration in parametric curves 1 Displacement in Parametric Curves - 1 Displacement in Parametric Curves Using
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 informationMechanics Oscillations Simple Harmonic Motion
Mechanics Oscillations Simple Harmonic Motion Lana Sheridan De Anza College Dec 3, 2018 Last time gravity Newton s universal law of gravitation gravitational field gravitational potential energy Overview
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 informationTIphysics.com. Physics. Pendulum Explorations ID: By Irina Lyublinskaya
Pendulum Explorations ID: 17 By Irina Lyublinskaya Time required 90 minutes Topic: Circular and Simple Harmonic Motion Explore what factors affect the period of pendulum oscillations. Measure the period
More informationThe Not-so-simple Pendulum: Balancing a Pencil on its Point Peter Lynch, UCD, Dublin, May 2014
The Not-so-simple Pendulum: Balancing a Pencil on its Point Peter Lynch, UCD, Dublin, May 204 arxiv:406.25v [nlin.si] 4 Jun 204 ABSTRACT. Does quantum mechanics matter at everyday scales? We generally
More informationLab 12: Periodic Motion
Lab 12: Periodic Motion Objectives: To devise an experiment to test variables that might affect the period of a pendulum To carry out an experiment testing variables that might affect the period of a pendulum,
More informationSolution Set Five. 2 Problem #2: The Pendulum of Doom Equation of Motion Simple Pendulum Limit Visualization...
: Solution Set Five Northwestern University, Classical Mechanics Classical Mechanics, Third Ed.- Goldstein November 4, 2015 Contents 1 Problem #1: Coupled Mass Oscillator System. 2 1.1 Normal Modes.......................................
More informationx = B sin ( t ) HARMONIC MOTIONS SINE WAVES AND SIMPLE HARMONIC MOTION Here s a nice simple fraction: y = sin (x) Differentiate = cos (x)
SINE WAVES AND SIMPLE HARMONIC MOTION Here s a nice simple fraction: y = sin (x) HARMONIC MOTIONS dy Differentiate = cos (x) dx So sin (x) has a stationary value whenever cos (x) = 0. 3 5 7 That s when
More informationChapter 14 Oscillations. Copyright 2009 Pearson Education, Inc.
Chapter 14 Oscillations 14-1 Oscillations of a Spring If an object vibrates or oscillates back and forth over the same path, each cycle taking the same amount of time, the motion is called periodic. The
More informationSimple and Physical Pendulums Challenge Problem Solutions
Simple and Physical Pendulums Challenge Problem Solutions Problem 1 Solutions: For this problem, the answers to parts a) through d) will rely on an analysis of the pendulum motion. There are two conventional
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 informationMechanics Cycle 1 Chapter 12. Chapter 12. Forces Causing Curved Motion
Chapter 1 Forces Causing Curved Motion A Force Must be Applied to Change Direction Coordinates, Angles, Angular Velocity, and Angular Acceleration Centripetal Acceleration and Tangential Acceleration Along
More informationThe Harmonic Oscillator
The Harmonic Oscillator Math 4: Ordinary Differential Equations Chris Meyer May 3, 008 Introduction The harmonic oscillator is a common model used in physics because of the wide range of problems it can
More informationClassification of Phase Portraits at Equilibria for u (t) = f( u(t))
Classification of Phase Portraits at Equilibria for u t = f ut Transfer of Local Linearized Phase Portrait Transfer of Local Linearized Stability How to Classify Linear Equilibria Justification of the
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 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 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 informationMITOCW R11. Double Pendulum System
MITOCW R11. Double Pendulum System The following content is provided under a Creative Commons license. Your support will help MIT OpenCourseWare continue to offer high quality educational resources for
More information11/17/10. Chapter 14. Oscillations. Chapter 14. Oscillations Topics: Simple Harmonic Motion. Simple Harmonic Motion
11/17/10 Chapter 14. Oscillations This striking computergenerated image demonstrates an important type of motion: oscillatory motion. Examples of oscillatory motion include a car bouncing up and down,
More informationLaboratory Instruction-Record Pages
Laboratory Instruction-Record Pages The Driven Pendulum Geology 200 - Evolutionary Systems James Madison University Lynn S. Fichter and Steven J. Baedke Brief History of Swinging Many things in this universe
More informationChapter 14 Oscillations
Chapter 14 Oscillations If an object vibrates or oscillates back and forth over the same path, each cycle taking the same amount of time, the motion is called periodic. The mass and spring system is a
More informationHow do the physical aspects of the oscillators affect the Period?
LAST NAME FIRST NAME DATE 10.4 The Pendulum & Spring Mass Oscillator Conceptual Questions 10, 11, 12, 13 page 314 Problems 40 page 317 How do the physical aspects of the oscillators affect the Period?
More informationNonlinear dynamics & chaos BECS
Nonlinear dynamics & chaos BECS-114.7151 Phase portraits Focus: nonlinear systems in two dimensions General form of a vector field on the phase plane: Vector notation: Phase portraits Solution x(t) describes
More information