2 Damped and Driven Harmonic Oscillator
Total Page:16
File Type:pdf, Size:1020Kb
Physics 12a Waves Lecture 4 Caltech, 10/11/19 2 Damped and Driven Harmonic Oscillator 2.2 Driven Harmonic Oscillation In order to have a persistent oscillatory motion in a system with dissipation, we can drive it by pushing on the system periodically. Consider again the classic example of a mass on a spring, but now with a horizontal force FD(t) as the drive. Newton's law now reads d2 d m x(t) = −kx(t) − mΓ x(t) + F (t) (1) dt2 dt D Reorganizing the terms, we get d2 d x(t) + Γ x(t) + !2x(t) = F (t)=m (2) dt2 dt 0 D 2 where we have used !0 = k=m. Let's consider first the case where FD(t) = F0 cos(!Dt) is a periodic force, but with a frequency !D that can be different from the intrinsic frequency !0. To solve the equation, we solve instead its complex version 2 d d F0 z(t) + Γ z(t) + !2z(t) = ei!Dt (3) dt2 dt 0 m i! t Let's try a solution of the form zD(t) = Ae D which describes an oscillatory motion at the driving frequency. Plugging this solution into the equation, we get F0 −!2 Aei!Dt + iΓA! ei!Dt + !2Aei!Dt = ei!Dt (4) D D 0 m From which we find F F 0 0 i!Dt A = 2 2 ; zD(t) = 2 2 e (5) m(!0 + iΓ!D − !D) m(!0 + iΓ!D − !D) and xD(t) = Re(A) cos(!Dt) − Im(A) sin(!Dt) (6) That is, xD(t) has two oscillatory parts, one is in phase with FD(t) and the other is 90 degree out of phase. The amplitude of the in phase part is 2 2 F0 j!0 − !Dj jRe(A)j = 2 2 2 2 2 (7) m (!0 − !D) + Γ !D 1 and the amplitude of the out of phase part is F0 Γ!D jIm(A)j = 2 2 2 2 2 (8) m (!0 − !D) + Γ !D Let's try to understand what this solution means: 1. The amplitude of both parts are proportional to F0, which makes sense. The larger the driving force, the larger the driven motion. F0 2. When !D << !0, Im(A) ≈ 0 while Re(A) ≈ 2 . The driven motion is given by m!0 F0 x(t) ≈ 2 cos(!Dt) (9) m!0 which is completely in phase with the driving force. 3. When !D >> !0 and !D >> Γ, both Re(A) and Im(A) approaches zero, but Im(A) goes to zero faster. So the driven motion is approximately given by F0 F0 x(t) ≈ − 2 cos(!Dt) = 2 cos(!Dt + π) (10) m!D m!D which is completely out of phase with the driving force. 4. Another important point in the parameter range of !D is when !0 = !D. At this point, Re(A) = 0, Im(A) = − F0 . The driven motion is given by mΓ!D F0 F0 π x(t) = sin(!Dt) = cos(!Dt − ) (11) mΓ!D mΓ!D 2 That is, the driven motion is 90 degree out of phase with the driving force. Across the full parameter range of 0 < !D < 1, the phase shift between the driving force and driven motion goes as While the displacement of the driven motion is 90 degree out of phase with the driving force when !D = !0, the velocity of the driven motion is completely in phase with the driving force. This can be seen from dx F v(t) = = 0 cos(! t) (12) dt mΓ D which is directly proportional to the driving force. As the driving force and the velocity are always in the same direction, the driving force is always doing positive work, the system absorbs the most energy from the driving force and goes into "resonance" and oscillates with (approximately) the largest amplitude. 2 Let's calculate how much energy is absorbed by the system in the driven motion for a generic value of !D. dx P = F · v = F · dt = F0 cos(!Dt)(−Re(A)!D sin(!Dt) − Im(A)!D cos(!Dt)) (13) Re(A) 2 = F0!D − 2 sin(2!Dt) − Im(A) cos (!Dt) Over time, the average power that is being absorbed is 2 2 1 F0 Γ !D < P >= − F0!DIm(A) = 4 2 2 2 4 (14) 2 2m !D + (Γ − 2!0)!D + !0 For what value of !D is < P > maximum? Notice that 2 F0 Γ 1 < P >= 4 (15) 2m 2 2 2 !0 !D + (Γ − 2!0) + 2 !D 4 2 !0 When !D = !0, !D + 2 reaches its minimum, therefore, < P > is at its max. !D If we plot < P > wrt to !D, we can see the peak of resonance, which has its maximum at !D = !0, 2 F0 with the maximum value being < P >max= 2mΓ . How sharp is the peak? We can characterize the width of the peak using the quantity called FWHM: 2 F0 Full Width at Half Maximum. As the maximum is < P >max= 2mΓ , the height at half maximum 2 2 F0 F0 is < P >half-max= 4mΓ . Correspondingly, we can find the two values of !D at which < P >= 4mΓ to be 1 q ! = ±Γ + Γ2 + 4!2 (16) 1=2 2 0 Therefore, the FWHM is simply given by Γ. 2 F0 To summarize, resonance occurs at !D = !0, with average absorbed power < P >= 2mΓ which is inversely proportional to Γ. The average absorbed power reaches its maximum because the displacement of oscillation is 90 degree out of phase with the driving force, therefore velocity of the oscillation is completely in phase with the driving force and the driving force is always doing positive work. Away from !D = !0, the absorbed power decays with FWHM Γ. An important quantity characterizing a real oscillator is its Q factor (quality factor) which is defined as ! Q = 0 (17) Γ which is the ratio between the intrinsic frequency and the friction constant. Note that as !0 and Γ have the same dimension, Q is a dimensionless number. With higher Q, the system is going to oscillate more times before the oscillation decays if it is not driven by external forces. At the same 3 time, it is going to have a narrower resonance peak, which means it is harder to tune to resonance and being more selective on frequency. Note that if we define resonance to be when the driven oscillation has the maximum amplitude, q 2 Γ2 then it happens at a slightly different location at !D = !0 − 2 . In practice, we can often ignore this difference as long as Γ is small. We are going to look in more detail about this in the homework. 4.