Phys42200 Lecture23
Phys42200 Lecture23
sin ≈
≈
ℓ
= ?
Damped Harmonic Motion
• Damping forces remove energy from the system
• We will only consider cases where the force is
proportional to the velocity: = − !
• You should be able to construct a free-body diagram
and write the resulting equation of motion:
+ " + = 0
– You should be able to write it in the standard form:
+ #" + = 0
• You must be able to solve this differential equation!
Damped Harmonic Motion
+ #" + = 0
Let = $ %&
• Characteristic polynomial:
' + #' + =0
• Roots (use the quadratic formula):
# #
'=− ± −
2 4
• Classification of solutions:
– Over-damped: # ⁄4 − >0 (distinct real roots)
– Critically damped: # ⁄4 = (one root)
– Under-damped: # ⁄4 − <0 (complex roots)
Damped Harmonic Motion
• Over-damped motion: # ⁄4 − >0
/ /0 0 / /0
.0& & . 23 .0& .& . 23 0
= $ $ 1+ $ $ 1
A
B(C) D
E
Sum of potential differences:
&
:; 1
−9 − ; < − ? + @ ; : = 0
: >
Initial charge, ? , defines the initial conditions.
Example
&
:; 1
9 +; <+ ? + @ ; : = 0
: >
Differentiate once with respect to time:
: ; :; 1
9 +< + ; =0
: : >
: ; :;
+# + ; =0
: :
Remember, the solution is B(C) but the initial
conditions might be in terms of F C = FG + H B C IC
(See examples from the lecture notes…)
Forced Harmonic Motion
• Now the differential equation is
+ " + = = cos
• Driving function is not always given in terms of a real
force… (think about non-inertial reference frames):
: K
J + #J" + J = − = > cos
:
• General properties:
– Steady state properties: ≫ 1/#
– Solution is J = cos −
– Amplitude, , and phase, , depend on
Forced Harmonic Motion
“Q” quantifies the amount of damping:
?=
#
(large Q means small damping force)
/
=
1
M/
− +
?
1/?
= tan.M
−
#
PQRR = −
( ⁄) 4
Q=5
Q=4
Q=3
Q=2
Q=1
/
Average Power
• The rate at which the oscillator absorbs energy is:
ST( ) =
1
2? 1
− +
?
ST( )
/2
? = 10
Full-Width-at-Half-Max:
?=5
[G
?=3 WXYZ = =\
F
?=1
/
Resonance
• Qualitative features: phase shift
1/?
= tan.M
−
→ 0 at low frequencies
→ ^ at high frequencies
`
_= when [ = [G
a
Coupled Oscillators
• Restoring force on
pendulum A:
ℓ ℓ b = −(b − c )
• Restoring force on
pendulum B:
c = (b − c )
c b
b + b + b − c = 0
ℓ
c + c − b − c = 0
ℓ
Coupled Oscillators
• You must be able to draw the free-body diagram
and set up the system of equations.
b + b + b − c = 0
ℓ
c + c − b − c = 0
ℓ
• You must be able to write this system as a matrix
equation.
b + d − d b ()
+ =0
c − d + d c ()
Coupled Oscillators
• Assume solutions are of the form
b () b
= cos −
c () c
• Then,
+ d − − d b
− d + d − c = 0
• You must be able to calculate the eigenvalues of a 2x2
or 3x3 matrix.
– Calculate the determinant
– Calculate the roots by factoring the determinant or using the
quadratic formula.
• These are the frequencies of the normal modes of
oscillation.
Coupled Oscillators
• You must be able to calculate the eigenvectors of a
2x2 or 3x3 matrix
• General solution:
e = feM cos M − g + he cos −i +⋯
• You must be able to solve for the constants of
integration using the initial conditions.
Coupled Discrete Systems
• The general method of calculating eigenvalues will always
work, but for simple systems you should be able to decouple
the equations by a change of variables.
b + b + b − c = 0
ℓ
1
ℓ ℓ c + c − b − c = 0
ℓ
b + + d b − d c = 0
2 c + + d c − d b = 0
= /ℓ, d = /
kM = b + c
3
k = b − c
c b k M + kM = 0
k + ′ k = 0
4
Forced Oscillations
• We mainly considered the qualitative aspects
– We did not analyze the behavior when damping forces
were significant
• Main features:
– Resonance occurs at each normal mode frequency
– Phase difference is = ^⁄2 at resonance
• Example: b driven by the force = mno
– Calculate force term applied to normal coordinates
M = = cos
– Reduced to two one-dimensional forced oscillators:
k M + kM = / cos
k + ′ k = / cos
Uniformly Distributed Discrete Systems
:;q 1 1
−9 − @ ;q − ;qpM : − @ ;q − ;q.M : = 0
: > >
: ;q
+2 ;q − (;q.M + ;qpM ) = 0
:
29
q =
t
t^!
q =
9
t!
q =
29
Boundary Conditions
• Examples:
– String fixed at both ends: J 0 = J 9 = 0
– Organ pipe open at one end: J" 0 = J" 9 = 0
• Driving end has maximal pressure amplitude
– Organ pipe closed at one end: J" 0 = 0, J 9 = 0
– Transmission line open at one end: ; 9 = 0
5
– Transmission line shorted at one end: ! 9 ∝ =0
&
Fourier Analysis
• Normal modes satisfying J 0 = J 9 = 0:
t^
Jq , = vq sin cos q − 'q
9
• General solution:
t^
J , = u vq sin cos q − 'q
9
qxM
• Initial conditions:
t^ t^
J , 0 = u vq sin cos 'q = u v′q sin
9 9
qxM qxM
t^ t^
J" , 0 = − u vq q sin sin 'q = u ′q sin
9 9
qxM qxM
Fourier Analysis
• Fourier sine transform:
t^
= u v′q sin
9
qxM
2 t^
v′q = @ () sin :
9 9
• Fourier cosine transform:
2 t^
′q = @ !() cos :
9 9
Fourier Analysis
v′q = vq cos 'q
′q = vq q sin 'q
Solve for amplitudes:
′q
vq = v′q +
q