# L = 1 2 a(q) q2 V (q).

Save this PDF as:
Size: px
Start display at page:

## Transcription

1 Physics 3550, Fall 2011 Motion near equilibrium - Small Oscillations Relevant Sections in Text: Motion near equilibrium 1 degree of freedom One of the most important situations in physics is motion near equilibrium. From a pendulum in a clock, to a bridge, to photons in quantum electrodynamics, the harmonic approximation for motion near equilibrium comes into play. We will now touch on some of the key ideas. You will recall that a Newtonian particle has an equilibrium point r 0 at locations where the force on the particle vanishes F ( r 0 ) = 0. In terms of a potential energy functions V ( r), an equilibrium point is a critical point of the potential energy: V ( r 0 ) = 0. So, let us begin by considering the motion of a one-dimensional system near a critical point of the potential energy. While only the simplest systems can be reduced to one degree of freedom, this case is still particularly important because as we shall soon see multi-dimensional systems (in the vicinity of stable equilibrium) can be reduced to multiple copies of the one dimensional case. We assume the Lagrangian is of the form: L = 1 2 a(q) q2 V (q). Note we assume (for now) that the system is closed, i.e., energy is conserved. Let q 0 be a critical point of the potential, so that V (q 0 ) = 0. Let us approximate the motion by assuming that x := q q 0 is small. We expand the Lagrangian in a Taylor series in x, keeping only the first non-trivial terms (exercise): L(x, ẋ) = 1 2 mẋ2 1 2 kx Here we have set m = a(q 0 ) and k = V (q 0 ). Of course, the kinetic energy must be positive, so that m > 0. Note that we also have dropped the irrelevant constant V (q 0 ). If desired we can adjust the zero point of potential energy so that this constant is zero. The EL equations are (in the domain of validity of the approximation) mẍ = kx. 1

2 This equation is easily integrated: x = Re(Ae iωt ), when k 0. x = at + b, if k = 0, where ω = ± km and A is a complex constant encoding the two real integration constants, which can be fixed by initial conditions. Taking the real part of the exponential expresses the solution in terms of a cosine. An alternative form of the exponential solution corresponds to using a sine: x = Im(Ae iωt ), when k 0. If k > 0, then q 0 is a point of stable equilibrium, and we get harmonic motion. In particular, if x is small initially and the initial velocity is sufficiently small, then x(t) remains small (exercise), so that our approximation is self-consistent. On the other hand, if k 0, then the motion of the particle need not maintain our approximation of small x; our approximation is not self-consistent and must be abandoned after a very short time. Of course, the case k < 0 corresponds to unstable equilibrium, for which a small perturbation leads to a rapid motion away from equilibrium. If k = 0, then the critical point q 0 is not a maximum or minimum but is a saddle point ( neutral equilibrium ); our approximation again becomes invalid, although the time scale for this is larger than the case of unstable equilibrium. Thus in a neighborhood of a point of stable equilibrium it is consistent to make the harmonic approximation to the potential and kinetic energies. In the harmonic approximation the motion of the system is mathematically the same as that of a simple harmonic oscillator. A very familiar example of all of this is the planar pendulum of mass m and length l for which q = θ is the deflection from a vertical position. We have V (θ) = mgl(1 cos θ), and a(θ) = 1 2 ml2. The Lagrangian in the harmonic approximation near equilibrium at θ = 0 is (exercise) L = 1 2 ml2 θ2 mgl 2 θ2, so that the harmonic motion has mass ml 2 and angular frequency ω = g l. (Exercise: what happens to motion near θ = π?) 2

3 Example Consider a mass m which is constrained to move on a straight line. The mass is bound to a fixed point by harmonic force with potential energy V = 2 1K(r R)2, where K is a constant, r is the distance of the particle to the fixed point. The distance from the point to the line is l > R. A mechanical model of this system is a mass sliding on a straight track; the mass being connected to a fixed point by a spring. Our goal is to find the stable equilibrium position(s) and compute the frequency of small oscillations about the equilibrium. Evidently, r 2 = x 2 + l 2, where x is the position of the particle along the given line, with x = 0 the location at distance l from the center of force. Using x as the generalized coordinate, the Lagrangian for this system is (exercise) L = 1 2 mẋ2 1 2 K( x 2 + l 2 R) 2. You can check that x = 0 is the only point of stable equilibrium. (As an exercise you can check that the point x = 0 is a point of unstable equilibrium if we assume R l.) You can expand the potential to second order to find that (exercise) k = K(l R), l so that K(l R) ω = lm. Note that K(l R) is the force needed to move a particle from r = R to the point r = l (exercise) in the presence of the potential V. Example We return to our example of a plane pendulum of length l with horizontally moving point of support. The mass m 1 (the point of support) has position x, and the angular deflection of the pendulum (mass m 2 ) is denoted by φ. The kinetic energy is T = 1 2 (m 1 + m 2 )ẋ m 2(l 2 φ2 + 2lẋ φ cos φ), and the potential energy is V = m 2 gl cos φ. We have already seen that we can use conservation of the momentum conjugate to x to effect a Lagrangian reduction which eliminates the x degree of freedom. In detail, the conservation law we need is P = (m 1 + m 2 )ẋ + m 2 l φ cos φ = constant. 3

4 This conservation law is just the conservation of the x component of the center of mass of the system (exercise). Let us compute the frequency of small oscillations of the pendulum in the reference frame in which P = 0, so that (exercise) (m 1 + m 2 )x + m 2 l sin φ = constant. In this reference frame, which is the rest frame of the x component of the center of mass, the reduced kinetic energy is (exercise) T = 1 2 m 2l 2 φ2 and the reduced potential energy is (exercise) ( ) m 1 2 cos 2 φ, m 1 + m 2 V = m 2 gl cos φ. Clearly φ = 0 is an equilibrium point (exercise). Expanding in powers of φ we get, in the harmonic approximation, T 1 m 1 m 2 l 2 2 m 1 + m φ2, 2 V 1 2 m 2glφ 2. It is now straightforward to see that the frequency of small oscillations is (exercise) ω = g(m 1 + m 2 ). m 1 l Of course, in a generic inertial reference frame the motion of the system in the harmonic approximation is small oscillations at the above frequency superimposed with a uniform translation of the center of mass along the x direction. As a sanity check, let us consider the limit in which m 2 m 1 << 1, i.e., the mass at the point of support is becoming large. Physically, we expect that the point of support moves with a uniform translation along x. In the rest frame of the point of support (which is now the approximate center of mass) we have a traditional plane pendulum problem. In this limit we get (exercise) g ω l, as expected. 4

5 Damped Oscillations Often times motion near equilibrium of a realistic system includes dissipation of energy into an environment. We will consider the case where this dissipation is proportional to velocity. There are other possibilities, of course. The equation of motion is mẍ + bẋ + kx = 0, with b > 0 controlling the dissipation. You can see that the effect of dissipation is always a deceleration. A simple example where such dissipation occurs is in an LRC electrical circuit. If q(t) is the charge on the capacitor at time t, we have L q + R q + 1 C q = 0. If we divide the equation of motion for x(t) by m, we have ẍ + 2βẋ + ω 2 x = 0, where β = 2m b. This is a linear ODE with constant coefficients. A standard approach to solving such an equation is to try to find 2 independent solutions of the form x(t) = Ce rt. Substituting this into the differential equation is straightforward and gives an equation for r: r 2 + 2βr + ω 2 = 0, which is easy to solve using the quadratic formula: r = β ± β 2 ω 2. Our two independent solutions are now x(t) = C 1 e βt+ β 2 ω 2t + C 2 e βt+ β 2 ω 2t. The constants C 1 and C 2 correspond to initial position and velocity. If β = 0, we see harmonic motion at frequency ω. If β < ω, we get damped, oscillatory motion. If β ω we get damped motion. You will explore this all a bit more in the homework. Forced Oscillations A typical scenario in which small oscillations is relevant is where one has a system in stable equilibrium which is subjected to an external force F which moves the system 5

6 from equilibrium. We allow this force to be time dependent. This introduces an additional potential energy term V 1 to the quadratically approximated Lagrangian given by We thus consider the Lagrangian V 1 = F (t) x. L = 1 2 mẋ2 1 2 kx2 + xf (t), where k > 0. The equations of motion are that of a harmonic oscillator subjected to an external, time varying force F (t) (exercise): ẍ + ω 2 x = F (t) m. Of course, we must assume that F remains suitably small so that the solutions do not violate the approximation needed for their validity. This inhomogeneous differential equation can be directly integrated as follows. Define ξ(t) = ẋ(t) + iωx(t). In terms of this complex variable the equation of motion takes the form (exercise) ξ iωξ = F (t) m. You can easily see that, when F = 0, the solution is of the form Ae iωt, where A is a constant. So try a solution of the form ξ(t) = A(t)e iωt. Plugging into the ODE, we find that A satisfies which has solution (exercise) A(t) = Ȧ = e iωt F (t) m, dt 1 m F (t)e iωt + constant. Putting this all together, we see that the general solution of the forced oscillator equation is (exercise) x(t) = 1 ω Im(ξ(t)) = 1 ω Im [e iωt (B + 6 ] 1 m F (t)e iωt dt),

7 where B is an arbitrary complex constant. Notice that the solution for x(t) is the sum of two terms: x(t) = x H (t) + x I (t), where and x H (t) = Im( B ω eiωt ), x I (t) = 1 ] [e ω Im iωt 1 m F (t)e iωt dt You will recognize x H (t) as the general solution to the un-forced harmonic oscillator equation. You can easily check that x I (t) is a solution to the forced harmonic oscillator equation. What we have here is a very general kind of result, valid when the equation being solved is linear, inhomogeneous. The general solution of a linear, inhomogeneous equation is the sum of any particular solution and the general solution of the associated homogeneous equation. As an example of all this, suppose that F (t) = f cos(γt). Then, provided γ 2 ω 2, the general solution is of the form (exercise) x(t) = a cos(ωt + α) + f m(ω 2 γ 2 ) cos(γt), where a and α are real constants. We see that in this case the motion is a superposition of two oscillations at the two frequencies ω and γ inherent in the problem. The relative importance of the forced oscillation component depends, of course, on the the size of f, but also on the relative magnitudes of ω and γ. When ω γ the forced oscillation amplitude diverges and our form of the solution given above becomes invalid; this situation is called resonance. To get the correct solution in this case we set γ = ω in our integral expression or the general solution. We then get a solution of the form (exercise) x(t) = a cos(ωt + α) + f t sin ωt, 2mω where, again, a and α are constants. Note the linear growth in t, which eventually destroys the harmonic approximation. We can get a useful picture of the behavior of the system near resonance by an approximation scheme. Let γ = ω + ɛ, 7

8 where ɛ << ω. Write the general solution for x(t) (off resonance) in the complex form (exercise): ξ(t) (A + Be iɛt )e iωt. Over one period, 2π ω, the amplitude C = A + Beiɛt changes very little. Thus the motion is approximately that of free oscillation with a slowly varying amplitude. In particular, the amplitude is of the form (exercise) C = a 2 + b 2 + 2ab cos(ɛt + φ), where A = ae iα, B = be iβ, and φ = β α. Thus the amplitude varies (slowly) between the values a + b and a b. The oscillatory behavior is said to exhibit beats. Typically, a general force F (t) can be Fourier analyzed into sinusoidal components. Likewise, we can Fourier analyze the solution x(t). We can view the above example as illustrating the behavior of a typical Fourier component. The general motion of the system is then a superposition of motions such as given above (exercise). Finally, let us note that since the Lagrangian for a system executing forced oscillations is explicitly time dependent (provided df dt 0), there will be no conservation of energy for the oscillator. This should not surprise you, since the oscillator is clearly exchanging energy with its environment. We can compute the energy transferred during a time interval (t 1, t 2 ) by noting that the oscillator energy can be written as and then using our explicit formula for ξ(t), E = 1 2 m(ẋ2 + ω 2 x 2 ) = 1 2 m ξ 2, t ξ = e iωt 1 (B + 0 m F (t)e iωt dt) to compute the energy at time t. For example, let us suppose that the system is at equilibrium before t = 0, a force acts for a period of time after t = 0 after which the force is zero again. Then B = 0 and the change in the oscillator energy can be written as (exercise) E = 1 2m F (t)e iωt 2 dt. Thus, the energy transfer is controlled by the absolute value of the Fourier component of the force with frequency ω. If the time during which the force acts is small compared to 1 ω, then eiωt is approximately constant in the integral, and hence E 1 2m 2 F (t) dt. 8

9 Here the change in energy is controlled solely by the impulse imparted by the force since the time scale is so short that no appreciable change in potential energy occurs while the force acts. In the limit where F (t) = fδ(t t 0 ), this approximation becomes exact (exercise). Finally, we remark that one can treat the forced oscillator with damping included in a similar way to what was done above. See your text for details. More than one degree of freedom Here we consider a couple of important and illustrative examples of oscillations with more than one degree of freedom. In general if V = V (q 1, q 2,..., q n ) is the potential energy function for a Newtonian system with n degrees of freedom. Then equilibrium q 0 is defined by V q i (q 0) = 0, i = 1, 2,..., n. and the small oscillation approximation involves the Taylor series of the potential energy to second order. Ignoring the zeroth order (constant) term, the first order term vanishes, and we have near equilibrium V 1 n 2 V 2 q i,j=1 i q j (q 0)x i x j, x i = q i q0 i. So, in this approximation the potential energy becomes a second order polynomial in the displacement from equilibrium. We will briefly explore a couple of standard situations which arise for this approximation. Isotropic and anisotropic oscillators We consider a system with two degrees of freedom (x, y) whose Lagrangian near equilibrium takes the form L = 1 2 m(ẋ2 + ẏ 2 ) 1 2 k xx k yy 2, k x, k y constants This corresponds to a matrix of second derivatives of the form 2 ( ) V q i q j (q kx 0 0) = 0 k. y 9

10 This system is known as a two-dimensional* anisotropic oscillator. An example where such a Lagrangian would be used is when studying a marble rolling near the bottom of a bowl which is sharply curved in one direction and less sharply curved in an orthogonal direction. Another example of such a Lagrangian (when generalized to 3-d) arises in condensed matter physics when studying the motion of an atom when displaced from equilibrium in a crystal which is not isotropic. The equations of motion for this Lagrangian are easily seen to be just two, decoupled harmonic oscillator equations: The solutions are mẍ + k x x = 0, mÿ + k y y = 0. x(t) = A cos(ω x t + α), y(t) = B cos(ω y t + β), where ω x = kx ky m, ω y = m and A, B, α, β are constants which can be determined by initial conditions. The motion in each direction is simple (harmonic), but the combined motion can be quite complex. To be sure, with initial position and velocity in one of the directions set to zero, the motion is just harmonic oscillations in the other direction. But, in the generic case, where ω x and ω y are two random constants and the initial conditions are not fine tuned in any way, the particle wanders in the plane with a motion which never repeats itself! It can be shown that if ω x and ω y are rationally related ω x /ω y is a rational number then the motion is periodic. If ω x = ω y then we speak of an isotropic oscillator. Now the motion is much more regular. Indeed you can check that, with ω x = ω y the motion is generically an ellipse. See your text for details. The isotropic oscillator provides an example of a rotationally invariant central force; we shall revisit the isotropic oscillator when we study the 2-body central force problem. Coupled oscillators Let us consider a Lagrangian of the form L = 1 2 m(ẋ2 + ẏ 2 ) 1 2 mω2 0 (x2 + y 2 ) 1 2 mα2 (x y) 2 = 1 ] [ẋ 2 m 2 + ẏ 2 (ω0 2 + α2 )(x 2 + y 2 ) + 2α 2 xy * I think you can easily generalize everything I will do here to more than two degrees of freedom. 10

11 This system can be viewed as two identical one-dimensional harmonic oscillators (natural frequency ω 0 ) with a coupling by a harmonic force (natural frequency α) (exercise). So, a toy model for this system could be two pendulums connected by a spring, or two masses each connected by identical springs to two walls and inter-connected by a different spring. In any case, this is the Lagrangian for motion near equilibrium whenever the potential energy is such that: 2 V q i q j (q 0) = The equations of motion can be written as: ( ) ( A B m(ω 2 = 0 + α 2 ) mα 2 B A mα 2 m(ω0 2 + α2 ) ẍ + (ω α2 )x α 2 y = 0, ÿ + (ω α2 )y α 2 x = 0. As you can see, the spring characterized by the frequency α couples the two onedimensional oscillators. How to solve these coupled, linear, homogeneous, second-order, constant coefficent, ordinary differential equations? A simple, linear change of coordinates will do the trick. Set q 1 = x + y,, q 2 = x y. Add and subtract the two equations of motion to get q 1 + ω 2 0 q 1 = 0, ). and q 2 + (ω α2 )q 2 = 0. The change of variables has decoupled the two ODEs and, moreover, has yield two 1-d harmonic oscillator equations! So, in term of the variables q 1 and q 2 the motion mathematically that of an anisotropic oscillator two decoupled simple harmonic oscillations as we treated in the previous subsection. We have the general solution q 1 (t) = C 1 cos(ω 0 t + β 1 ), q 2 (t) = C 2 cos( ω α2 t + β 2 ), or x(t) = 1 ) (C 1 cos(ω 1 t + β 1 ) + C 2 cos(ω 2 t + β 2 ) 2 y(t) = 1 ) (C 1 cos(ω 1 t + β 1 ) C 2 cos(ω 2 t + β 2 ). 2 It can be shown that this type of result the motion can be viewed as a combination of independent harmonic oscillations arises whenever one approximates motion near 11

12 equilibrium as we have done here, for any number of degrees of freedom! The harmonic motions are called normal modes of vibration and the associated frequencies are the characteristic frequencies. Evidently, the normal mode at frequency ω 0 corresponds to the two masses moving exactly in phase, so the coupling does not come into play, while the oscillation at frequency ω α2 has the two masses moving exactly out of phase. The general motion of the system is a superposition of these two normal modes; the particular superposition which occurs is determined by initial conditions. In particular, if we choose the initial conditions just right we can have motion which is purely one or the other normal mode. For example, normal mode q 2 is excited by itself if one picks initial x and y to be equal and opposite and with vanishing initial velocity. Can you see this? For generic initial conditions the motion can be quite complicated as we discussed for the anisotropic oscillator. 12

### L = 1 2 a(q) q2 V (q).

Physics 3550 Motion near equilibrium - Small Oscillations Relevant Sections in Text: 5.1 5.6, 11.1 11.3 Motion near equilibrium 1 degree of freedom One of the most important situations in physics is motion

### Damped 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,

### Coupled Oscillators Monday, 29 October 2012 normal mode Figure 1:

Coupled Oscillators Monday, 29 October 2012 In which we count degrees of freedom and find the normal modes of a mess o masses and springs, which is a lovely model of a solid. Physics 111 θ 1 l m k θ 2

### Seminar 6: COUPLED HARMONIC OSCILLATORS

Seminar 6: COUPLED HARMONIC OSCILLATORS 1. Lagrangian Equations of Motion Let consider a system consisting of two harmonic oscillators that are coupled together. As a model, we will use two particles attached

### In which we count degrees of freedom and find the normal modes of a mess o masses and springs, which is a lovely model of a solid.

Coupled Oscillators Wednesday, 26 October 2011 In which we count degrees of freedom and find the normal modes of a mess o masses and springs, which is a lovely model of a solid. Physics 111 θ 1 l k θ 2

### Damped Harmonic Oscillator

Damped Harmonic Oscillator Wednesday, 23 October 213 A simple harmonic oscillator subject to linear damping may oscillate with exponential decay, or it may decay biexponentially without oscillating, or

### the EL equation for the x coordinate is easily seen to be (exercise)

Physics 6010, Fall 2016 Relevant Sections in Text: 1.3 1.6 Examples After all this formalism it is a good idea to spend some time developing a number of illustrative examples. These examples represent

### Mechanics 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

### The 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

### 4. 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

### Chapter 6. Second order differential equations

Chapter 6. Second order differential equations A second order differential equation is of the form y = f(t, y, y ) where y = y(t). We shall often think of t as parametrizing time, y position. In this case

### 1 Simple Harmonic Oscillator

Physics 1a Waves Lecture 3 Caltech, 10/09/18 1 Simple Harmonic Oscillator 1.4 General properties of Simple Harmonic Oscillator 1.4.4 Superposition of two independent SHO Suppose we have two SHOs described

### M2A2 Problem Sheet 3 - Hamiltonian Mechanics

MA Problem Sheet 3 - Hamiltonian Mechanics. The particle in a cone. A particle slides under gravity, inside a smooth circular cone with a vertical axis, z = k x + y. Write down its Lagrangian in a) Cartesian,

### Thursday, August 4, 2011

Chapter 16 Thursday, August 4, 2011 16.1 Springs in Motion: Hooke s Law and the Second-Order ODE We have seen alrealdy that differential equations are powerful tools for understanding mechanics and electro-magnetism.

### 01 Harmonic Oscillations

Utah State University DigitalCommons@USU Foundations of Wave Phenomena Library Digital Monographs 8-2014 01 Harmonic Oscillations Charles G. Torre Department of Physics, Utah State University, Charles.Torre@usu.edu

### Problem 1: Lagrangians and Conserved Quantities. Consider the following action for a particle of mass m moving in one dimension

105A Practice Final Solutions March 13, 01 William Kelly Problem 1: Lagrangians and Conserved Quantities Consider the following action for a particle of mass m moving in one dimension S = dtl = mc dt 1

### Lecture 9: Eigenvalues and Eigenvectors in Classical Mechanics (See Section 3.12 in Boas)

Lecture 9: Eigenvalues and Eigenvectors in Classical Mechanics (See Section 3 in Boas) As suggested in Lecture 8 the formalism of eigenvalues/eigenvectors has many applications in physics, especially in

### Assignment 6. Using the result for the Lagrangian for a double pendulum in Problem 1.22, we get

Assignment 6 Goldstein 6.4 Obtain the normal modes of vibration for the double pendulum shown in Figure.4, assuming equal lengths, but not equal masses. Show that when the lower mass is small compared

### Linear Second Order ODEs

Chapter 3 Linear Second Order ODEs In this chapter we study ODEs of the form (3.1) y + p(t)y + q(t)y = f(t), where p, q, and f are given functions. Since there are two derivatives, we might expect that

### P321(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

### Linear 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

### Analytical Mechanics ( AM )

Analytical Mechanics ( AM ) Olaf Scholten KVI, kamer v8; tel nr 6-55; email: scholten@kvinl Web page: http://wwwkvinl/ scholten Book: Classical Dynamics of Particles and Systems, Stephen T Thornton & Jerry

### Consider a particle in 1D at position x(t), subject to a force F (x), so that mẍ = F (x). Define the kinetic energy to be.

Chapter 4 Energy and Stability 4.1 Energy in 1D Consider a particle in 1D at position x(t), subject to a force F (x), so that mẍ = F (x). Define the kinetic energy to be T = 1 2 mẋ2 and the potential energy

### 1 Pushing your Friend on a Swing

Massachusetts Institute of Technology MITES 017 Physics III Lecture 05: Driven Oscillations In these notes, we derive the properties of both an undamped and damped harmonic oscillator under the influence

### ! 4 4! o! +! h 4 o=0! ±= ± p i And back-substituting into the linear equations gave us the ratios of the amplitudes of oscillation:.»» = A p e i! +t»»

Topic 6: Coupled Oscillators and Normal Modes Reading assignment: Hand and Finch Chapter 9 We are going to be considering the general case of a system with N degrees of freedome close to one of its stable

### Physics 6010, Fall Relevant Sections in Text: Introduction

Physics 6010, Fall 2016 Introduction. Configuration space. Equations of Motion. Velocity Phase Space. Relevant Sections in Text: 1.1 1.4 Introduction This course principally deals with the variational

### A 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

### Variation Principle in Mechanics

Section 2 Variation Principle in Mechanics Hamilton s Principle: Every mechanical system is characterized by a Lagrangian, L(q i, q i, t) or L(q, q, t) in brief, and the motion of he system is such that

### Keble College - Hilary 2014 CP3&4: Mathematical methods I&II Tutorial 5 - Waves and normal modes II

Tomi Johnson 1 Keble College - Hilary 2014 CP3&4: Mathematical methods I&II Tutorial 5 - Waves and normal modes II Prepare full solutions to the problems with a self assessment of your progress on a cover

### Introductory Physics. Week 2015/05/29

2015/05/29 Part I Summary of week 6 Summary of week 6 We studied the motion of a projectile under uniform gravity, and constrained rectilinear motion, introducing the concept of constraint force. Then

### Chapter 2: Complex numbers

Chapter 2: Complex numbers Complex numbers are commonplace in physics and engineering. In particular, complex numbers enable us to simplify equations and/or more easily find solutions to equations. We

### Physics 351 Monday, January 22, 2018

Physics 351 Monday, January 22, 2018 Phys 351 Work on this while you wait for your classmates to arrive: Show that the moment of inertia of a uniform solid sphere rotating about a diameter is I = 2 5 MR2.

### Selected Topics in Physics a lecture course for 1st year students by W.B. von Schlippe Spring Semester 2007

Selected Topics in Physics a lecture course for st year students by W.B. von Schlippe Spring Semester 7 Lecture : Oscillations simple harmonic oscillations; coupled oscillations; beats; damped oscillations;

### 221B Lecture Notes on Resonances in Classical Mechanics

1B Lecture Notes on Resonances in Classical Mechanics 1 Harmonic Oscillators Harmonic oscillators appear in many different contexts in classical mechanics. Examples include: spring, pendulum (with a small

### Oscillations and Waves

Oscillations and Waves Somnath Bharadwaj and S. Pratik Khastgir Department of Physics and Meteorology IIT Kharagpur Module : Oscillations Lecture : Oscillations Oscillations are ubiquitous. It would be

### 1. (10 points) Find the general solution to the following second-order differential equation:

Math 307A, Winter 014 Midterm Solutions Page 1 of 8 1. (10 points) Find the general solution to the following second-order differential equation: 4y 1y + 9y = 9t. To find the general solution to this nonhomogeneous

### Oscillations. Simple Harmonic Motion of a Mass on a Spring The equation of motion for a mass m is attached to a spring of constant k is

Dr. Alain Brizard College Physics I (PY 10) Oscillations Textbook Reference: Chapter 14 sections 1-8. Simple Harmonic Motion of a Mass on a Spring The equation of motion for a mass m is attached to a spring

### Wave Phenomena Physics 15c. Lecture 2 Damped Oscillators Driven Oscillators

Wave Phenomena Physics 15c Lecture Damped Oscillators Driven Oscillators What We Did Last Time Analyzed a simple harmonic oscillator The equation of motion: The general solution: Studied the solution m

### In the presence of viscous damping, a more generalized form of the Lagrange s equation of motion can be written as

2 MODELING Once the control target is identified, which includes the state variable to be controlled (ex. speed, position, temperature, flow rate, etc), and once the system drives are identified (ex. force,

### 28. 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.

### Solution 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.......................................

### Resonance 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.

### AA242B: MECHANICAL VIBRATIONS

AA242B: MECHANICAL VIBRATIONS 1 / 50 AA242B: MECHANICAL VIBRATIONS Undamped Vibrations of n-dof Systems These slides are based on the recommended textbook: M. Géradin and D. Rixen, Mechanical Vibrations:

### Math 3313: Differential Equations Second-order ordinary differential equations

Math 3313: Differential Equations Second-order ordinary differential equations Thomas W. Carr Department of Mathematics Southern Methodist University Dallas, TX Outline Mass-spring & Newton s 2nd law Properties

### The two body problem involves a pair of particles with masses m 1 and m 2 described by a Lagrangian of the form:

Physics 3550, Fall 2011 Two Body, Central-Force Problem Relevant Sections in Text: 8.1 8.7 Two Body, Central-Force Problem Introduction. I have already mentioned the two body central force problem several

### 1 Review of simple harmonic oscillator

MATHEMATICS 7302 (Analytical Dynamics YEAR 2017 2018, TERM 2 HANDOUT #8: COUPLED OSCILLATIONS AND NORMAL MODES 1 Review of simple harmonic oscillator In MATH 1301/1302 you studied the simple harmonic oscillator:

### Vibrations and waves: revision. Martin Dove Queen Mary University of London

Vibrations and waves: revision Martin Dove Queen Mary University of London Form of the examination Part A = 50%, 10 short questions, no options Part B = 50%, Answer questions from a choice of 4 Total exam

### Physics 6010, Fall 2016 Constraints and Lagrange Multipliers. Relevant Sections in Text:

Physics 6010, Fall 2016 Constraints and Lagrange Multipliers. Relevant Sections in Text: 1.3 1.6 Constraints Often times we consider dynamical systems which are defined using some kind of restrictions

### Undetermined Coefficents, Resonance, Applications

Undetermined Coefficents, Resonance, Applications An Undetermined Coefficients Illustration Phase-amplitude conversion I Phase-amplitude conversion II Cafe door Pet door Cafe Door Model Pet Door Model

### Classical Mechanics Phys105A, Winter 2007

Classical Mechanics Phys5A, Winter 7 Wim van Dam Room 59, Harold Frank Hall vandam@cs.ucsb.edu http://www.cs.ucsb.edu/~vandam/ Phys5A, Winter 7, Wim van Dam, UCSB Midterm New homework has been announced

### T = 1 2 mv2. For a system of N particles, the kinetic energy of the system is obtained by adding up the individual kinetic energies: T = 2 m ivi 2.

Physics 3550, Fall 2011 Energy. Relevant Sections in Text: Chapter 4 Kinetic Energy. I probably don t need to impress upon you the importance of the notion of energy in analyzing/understanding physical

University of Rhode Island DigitalCommons@URI Classical Dynamics Physics Course Materials 015 14. Oscillations Gerhard Müller University of Rhode Island, gmuller@uri.edu Creative Commons License This wor

### Damped Oscillation Solution

Lecture 19 (Chapter 7): Energy Damping, s 1 OverDamped Oscillation Solution Damped Oscillation Solution The last case has β 2 ω 2 0 > 0. In this case we define another real frequency ω 2 = β 2 ω 2 0. In

### Vibrations and Waves Physics Year 1. Handout 1: Course Details

Vibrations and Waves Jan-Feb 2011 Handout 1: Course Details Office Hours Vibrations and Waves Physics Year 1 Handout 1: Course Details Dr Carl Paterson (Blackett 621, carl.paterson@imperial.ac.uk Office

### Lecture 6: Differential Equations Describing Vibrations

Lecture 6: Differential Equations Describing Vibrations In Chapter 3 of the Benson textbook, we will look at how various types of musical instruments produce sound, focusing on issues like how the construction

### Topic 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

### Chapter 1. Harmonic Oscillator. 1.1 Energy Analysis

Chapter 1 Harmonic Oscillator Figure 1.1 illustrates the prototypical harmonic oscillator, the mass-spring system. A mass is attached to one end of a spring. The other end of the spring is attached to

### NORMAL MODES, WAVE MOTION AND THE WAVE EQUATION. Professor G.G.Ross. Oxford University Hilary Term 2009

NORMAL MODES, WAVE MOTION AND THE WAVE EQUATION Professor G.G.Ross Oxford University Hilary Term 009 This course of twelve lectures covers material for the paper CP4: Differential Equations, Waves and

### PHY 3221 Fall Homework Problems. Instructor: Yoonseok Lee. Submit only HW s. EX s are additional problems that I encourage you to work on.

PHY 3221 Fall 2012 Homework Problems Instructor: Yoonseok Lee Submit only HW s. EX s are additional problems that I encourage you to work on. Week 1: August 22-24, Due August 27 (nothing to submit) EX:

### Mass on a spring: including a driving force, and friction proportional to velocity (eg, from a dashpot).

Physics 06a, Caltech 8 October, 208 Lecture 6: Equilibria and Oscillations We study equilibria and linear oscillations. Motion of dynamical systems can be classified as secular (long-term, non-periodic)

### Honors Differential Equations

MIT OpenCourseWare http://ocw.mit.edu 18.034 Honors Differential Equations Spring 009 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. LECTURE 7. MECHANICAL

### Physics 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

### 2 Canonical quantization

Phys540.nb 7 Canonical quantization.1. Lagrangian mechanics and canonical quantization Q: How do we quantize a general system?.1.1.lagrangian Lagrangian mechanics is a reformulation of classical mechanics.

### Two-Body Problem. Central Potential. 1D Motion

Two-Body Problem. Central Potential. D Motion The simplest non-trivial dynamical problem is the problem of two particles. The equations of motion read. m r = F 2, () We already know that the center of

### Physics 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

### dx n a 1(x) dy

HIGHER ORDER DIFFERENTIAL EQUATIONS Theory of linear equations Initial-value and boundary-value problem nth-order initial value problem is Solve: a n (x) dn y dx n + a n 1(x) dn 1 y dx n 1 +... + a 1(x)

### Physics 5153 Classical Mechanics. Canonical Transformations-1

1 Introduction Physics 5153 Classical Mechanics Canonical Transformations The choice of generalized coordinates used to describe a physical system is completely arbitrary, but the Lagrangian is invariant

### Classical mechanics of particles and fields

Classical mechanics of particles and fields D.V. Skryabin Department of Physics, University of Bath PACS numbers: The concise and transparent exposition of many topics covered in this unit can be found

### Forced Oscillation and Resonance

Chapter Forced Oscillation and Resonance The forced oscillation problem will be crucial to our understanding of wave phenomena Complex exponentials are even more useful for the discussion of damping and

### for changing independent variables. Most simply for a function f(x) the Legendre transformation f(x) B(s) takes the form B(s) = xs f(x) with s = df

Physics 106a, Caltech 1 November, 2018 Lecture 10: Hamiltonian Mechanics I The Hamiltonian In the Hamiltonian formulation of dynamics each second order ODE given by the Euler- Lagrange equation in terms

### Parametric Excitation of a Linear Oscillator

Parametric Excitation of a Linear Oscillator Manual Eugene Butikov Annotation. The manual includes a description of the simulated physical system and a summary of the relevant theoretical material for

### Assignments VIII and IX, PHYS 301 (Classical Mechanics) Spring 2014 Due 3/21/14 at start of class

Assignments VIII and IX, PHYS 301 (Classical Mechanics) Spring 2014 Due 3/21/14 at start of class Homeworks VIII and IX both center on Lagrangian mechanics and involve many of the same skills. Therefore,

### 4.2 Homogeneous Linear Equations

4.2 Homogeneous Linear Equations Homogeneous Linear Equations with Constant Coefficients Consider the first-order linear differential equation with constant coefficients a 0 and b. If f(t) = 0 then this

### 22.2. Applications of Eigenvalues and Eigenvectors. Introduction. Prerequisites. Learning Outcomes

Applications of Eigenvalues and Eigenvectors 22.2 Introduction Many applications of matrices in both engineering and science utilize eigenvalues and, sometimes, eigenvectors. Control theory, vibration

### Identification Methods for Structural Systems. Prof. Dr. Eleni Chatzi Lecture March, 2016

Prof. Dr. Eleni Chatzi Lecture 4-09. March, 2016 Fundamentals Overview Multiple DOF Systems State-space Formulation Eigenvalue Analysis The Mode Superposition Method The effect of Damping on Structural

### Rotational 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

### Non-Linear Response of Test Mass to External Forces and Arbitrary Motion of Suspension Point

LASER INTERFEROMETER GRAVITATIONAL WAVE OBSERVATORY -LIGO- CALIFORNIA INSTITUTE OF TECHNOLOGY MASSACHUSETTS INSTITUTE OF TECHNOLOGY Technical Note LIGO-T980005-01- D 10/28/97 Non-Linear Response of Test

### Solutions 2: Simple Harmonic Oscillator and General Oscillations

Massachusetts Institute of Technology MITES 2017 Physics III Solutions 2: Simple Harmonic Oscillator and General Oscillations Due Wednesday June 21, at 9AM under Rene García s door Preface: This problem

### Springs: Part I Modeling the Action The Mass/Spring System

17 Springs: Part I Second-order differential equations arise in a number of applications We saw one involving a falling object at the beginning of this text (the falling frozen duck example in section

### Oscillations Simple Harmonic Motion

Oscillations Simple Harmonic Motion Lana Sheridan De Anza College Dec 1, 2017 Overview oscillations simple harmonic motion (SHM) spring systems energy in SHM pendula damped oscillations Oscillations and

### Math 312 Lecture Notes Linear Two-dimensional Systems of Differential Equations

Math 2 Lecture Notes Linear Two-dimensional Systems of Differential Equations Warren Weckesser Department of Mathematics Colgate University February 2005 In these notes, we consider the linear system of

### MIT Weakly Nonlinear Things: Oscillators.

18.385 MIT Weakly Nonlinear Things: Oscillators. Department of Mathematics Massachusetts Institute of Technology Cambridge, Massachusetts MA 02139 Abstract When nonlinearities are small there are various

### Work sheet / Things to know. Chapter 3

MATH 251 Work sheet / Things to know 1. Second order linear differential equation Standard form: Chapter 3 What makes it homogeneous? We will, for the most part, work with equations with constant coefficients

### Math 211. Substitute Lecture. November 20, 2000

1 Math 211 Substitute Lecture November 20, 2000 2 Solutions to y + py + qy =0. Look for exponential solutions y(t) =e λt. Characteristic equation: λ 2 + pλ + q =0. Characteristic polynomial: λ 2 + pλ +

### Wave Phenomena Physics 15c

Wave Phenomena Physics 15c Lecture Harmonic Oscillators (H&L Sections 1.4 1.6, Chapter 3) Administravia! Problem Set #1! Due on Thursday next week! Lab schedule has been set! See Course Web " Laboratory

### kg meter ii) Note the dimensions of ρ τ are kg 2 velocity 2 meter = 1 sec 2 We will interpret this velocity in upcoming slides.

II. Generalizing the 1-dimensional wave equation First generalize the notation. i) "q" has meant transverse deflection of the string. Replace q Ψ, where Ψ may indicate other properties of the medium that

### Physics 351, Spring 2017, Homework #3. Due at start of class, Friday, February 3, 2017

Physics 351, Spring 2017, Homework #3. Due at start of class, Friday, February 3, 2017 Course info is at positron.hep.upenn.edu/p351 When you finish this homework, remember to visit the feedback page at

### HB 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

### Wave Phenomena Physics 15c

Wave Phenomena Physics 15c Lecture Harmonic Oscillators (H&L Sections 1.4 1.6, Chapter 3) Administravia! Problem Set #1! Due on Thursday next week! Lab schedule has been set! See Course Web " Laboratory

### Forced Mechanical Vibrations

Forced Mechanical Vibrations Today we use methods for solving nonhomogeneous second order linear differential equations to study the behavior of mechanical systems.. Forcing: Transient and Steady State

### Second-Order Linear ODEs

C0.tex /4/011 16: 3 Page 13 Chap. Second-Order Linear ODEs Chapter presents different types of second-order ODEs and the specific techniques on how to solve them. The methods are systematic, but it requires

### x(t+ δt) - x(t) = slope δt t+δt

Techniques of Physics Worksheet 2 Classical Vibrations and Waves Introduction You will have encountered many different examples of wave phenomena in your courses and should be familiar with most of the

### PHYSICS. 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

### Chapter 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

### Imaginary. Axis. Real. Axis

Name ME6 Final. I certify that I upheld the Stanford Honor code during this exam Monday December 2, 25 3:3-6:3 p.m. ffl Print your name and sign the honor code statement ffl You may use your course notes,

### Classical Mechanics Comprehensive Exam Solution

Classical Mechanics Comprehensive Exam Solution January 31, 011, 1:00 pm 5:pm Solve the following six problems. In the following problems, e x, e y, and e z are unit vectors in the x, y, and z directions,

### Differential Equations

Differential Equations A differential equation (DE) is an equation which involves an unknown function f (x) as well as some of its derivatives. To solve a differential equation means to find the unknown

### Math 221 Topics since the second exam

Laplace Transforms. Math 1 Topics since the second exam There is a whole different set of techniques for solving n-th order linear equations, which are based on the Laplace transform of a function. For