PHY5210 W15

Chapter 11: Coupled Oscillators and Normal Modes

Reading:

In Taylor, read sections 11.3 to 11.4 for today, and 11.5 for Wednesday.

Recall

For two masses and three springs we found

M xddot = -K x

where M is the diagonal matrix [(m1, 0), (0, m2)] and K is the matrix [((k1 + k2), -k2), (-k2, (k2 + k3))].

Two Weakly Coupled Oscillators: One Spring Different from the Other Two

Another situation of particular interest occurs when a set of (identical) harmonic oscillators are coupled, but only weakly. This is a technique commonly applied to solving problem: if you have a difficult problem that involves two simpler situations that you do know how to solve, coupled together, first consider the problem with weak coupling. In our simple model, this situation would correspond to the case of two equal masses, and identical springs 1 and 3, k1 = k3 = k, and a much weaker spring coupling them, k2 << k. This implies that at some point we will make an approximation and drop terms that are small.

Before beginning, we expect to find solutions that are close to two uncoupled oscillators. There will be some small effect due to the weak coupling.

The matrix M is unchanged from before, and we can now write K as K = [((k + k2), -k2), (-k2, (k + k2))]. The determinant equation becomes

(k + k2 - m ω²)² - k2² = (k - m ω²)² + 2k2 (k - m ω²) = (k - m ω²)(k + 2k2 - m ω²) = 0

This equation yields two normal mode frequencies

ω1 = √(k/m)

and

ω2 = √((k + 2k2)/m).

Thus far we've made no approximations. In the weak coupling limit, k2 << k, we find two normal frequencies that are nearly equal, as suggested at the beginning. The first frequency is the same as in the first example, and the same frequency we would expect for two independent, but identical, oscillators. The effect of the weak coupling is expressed in the second frequency ω2 being slightly higher than the uncoupled frequency. It is convenient to express this result in a different way that emphasizes the closeness of ω1 and ω2. Start by taking the average of ω1 and ω2

ω0 = (ω1 + ω2)/2

Now we can write the two normal frequencies as being a little below and a little above the average frequency ω0

ω1 = ω0 - ε  and  ω2 = ω0 + ε

where some calculation will reveal that ε ≈ k2/2m.

x = C1 [1, 1] ei0 - ε)t  and  x = C2 [1, -1] ei0 + ε)t.
x(t) = {C1 [1, 1] e-iεt + C2 [1, -1] eiεt} eiω0t.
Let C1 = C2 = A/2, then the position as a function of time becomes
x(t) = A [cos εt, -i sin εt] eiω0t.

The actual motion is the real part of x (or equally valid, the imaginary part). Taking the real part gives

x1(t) = A cos εt cos ωt
x2(t) = A sin εt sin ωt

Let's look at this result in detail. Since ε << ω we can consider the factors A cos εt and A sin εt as a slowly changing amplitude for the cosωt and sinωt oscillations. That is, we can approximate the overall motion as (rapid) oscillations at the uncoupled frequency occurring with an envelope of slowly varying amplitude. The amplitude envelopes vary out of phase, so that first mass 1 swings at full amplitude while mass 2 is at rest; then the energy in mass 1 is transferred to mass 2, with a corresponding increase of the amplitude of mass 2 and reduction of amplitude of mass 1; until finally mass 2 swings at full amplitude and mass 1 comes to rest; then the whole process repeats but with energy being transferred from mass 2 to mass 1.

This phenomenon is known as beats. It arises when two modes exist that are close in frequency. It occurs with sound waves, electrical waves, and quantum mechanical probability waves.

Lagrangian Approach

Lagrangian Approach for Two Carts and Three Springs

Already done in lecture.

The Double Pendulum

A double pendulum is a pendulum suspended from a pendulum. Call the upper mass m1, the lower mass m2, the length of the upper pendulum is L1 and the length of the lower pendulum is L2. The upper pendulum hangs at an angle φ1 to the vertical and the lower makes an angle φ2 to the vertical.

As is the case for the single pendulum, it is possible to analytically solve for the motion only in the case of small angle oscillations, so we will make the small angle approximation.

The potential energy of the upper pendulum is

U1 = m1 g L1 (1 - cosφ1)

The potential energy of the lower pendulum is

U2 = m2 g [L1 (1 - cosφ1) + L2 (1 - cosφ2)]

In the small angle approximation, the total potential energy is

U1, φ2) = ½ (m1 + m2) gL1φ1² + ½m2 gL2φ2².

For the first mass, the position is (x1, y1) = L1 (sinφ1, 1-cosφ1), and the kinetic energy is T1 = ½ m1 L1² φ1dot².

The position of the second mass is (x2, y2) = (x1, y1) + L2 (sinφ2, 1-cosφ2). The square of the velocity is v2² = x2dot² + y2dot² = L1² φ1dot² + L2² φ2dot² + L1 L2 φ1dot φ2dot (cosφ1 cosφ2 + sinφ1 sinφ2). The final trigonometric terms can be written as cos(φ1 - φ2) ≈ 1 for small angles.

Putting these together, we can write the kinetic energy as

T = ½(m1 + m2) L1² φ1dot² + ½m2 L2² φ2dot² + m2 L1 L2 φ1dot φ2dot.

And the Lagrangian is T - U. Lagrange's equations yield, for φ1

(m1 + m2) L1² φ1ddot + m2 L1 L2 φ2ddot = -(m1 + m2) g L1 φ1

and for φ2

m2 L1 L2 φ1ddot + m2 L2² φ2ddot = -m2 g L1 φ2.

These can be written in our standard matrix form

M φddot = -K φ

where

φ = [φ1, φ2]

and

M = [((m1 + m2) L1², m2 L1 L2), (m2 L1 L2, m2 L2²)]

and

K = [((m1 + m2) g L1, 0), (0, m2 g L1)].

Try the same trick again to solve this equation: try the solution φ = aeiωt, and get the characteristic equation (K - ω² M) a = 0 which must be solved for ω and a.

Equal Lengths and Masses

As a simple example, consider the case when the masses and lengths are equal. This simplifies K and M to

M = m L² [(2, 1), (1, 1)]

and

K = m L² ω0² [(2, 0), (0, 1)]

where ω0 = √(g/L) is the natural frequency of a pendulum of length L. The characteristic equation (determinant) is

m² L4 [2(ω0² - ω²)² - ω4] = m² L44 - 4ω0²ω² + 2ω0²] = 0

This has solutions

ω1² = (2 - √2)ω0²  and  ω2² = (2 + √2)ω0²

The normal modes are found by substituting back into the characteristic equation. For ω1 we get

(K - ω² M) a = m L² ω0² (√2 - 1) [(2, -√2), (-√2, 1)] [a1, a2] = 0,

implying that a2 = √2 a1. For ω2 we get

(K - ω² M) a = m L² ω0² (√2 + 1) [(2, √2), (√2, 1)] [a1, a2] = 0,

implying that a2 = -√2 a1.

In the first mode, the masses oscillate together, but the amplitude of the lower mass is √2 times larger than the amplitude of the upper mass. In the second mode, the masses oscillate 180° out of phase, again with the lower mass swinging with an amplitude √2 times greater than the upper mass.

© 2015 Robert Harr