- +1-530-264-8006
- info@tutorsglobe.com

18,76,764

Questions

Asked

21,311

Experts

9,67,568

Questions

Answered

Start Excelling in your courses, Ask an Expert and get answers for your homework and assignments!!

Submit Assignment2015 © Tutors Globe. All rights reserved.

## Damped Harmonic Motion, Physics tutorial

Introduction:In real systems, masses on springs do not continue to oscillate forever at same amplitude; ultimately oscillations die away and object stops. This is because of fact that springs are not very effective at storing and releasing energy. Much of the energy is dissipated as heat because of friction within spring.

Let us suppose that mass moves horizontally in the viscous medium, say inside lubricated cylinder, as shown in figure given below. As mass moves, it will experience drag, that we denote by F

_{d}. For oscillations of adequately small amplitude, it is quite sensible to model damping force after Stokes' law. That is, we take F_{d}to be proportional to velocity and writeF

_{d}= -γv.....................................Eq.1The negative sign signifies that damping force opposes motion. Constant of proportionality γ is known as damping coefficient. Numerically, it is equal to force per unit velocity and is estimated in N/ms

^{-1}= kgms^{-2}/ms^{-1}= kgs^{-1}Let us take the x -axis to be along length of spring. We state origin of axis (x = 0) as equilibrium position of mass. Imagine that mass (in spring-mass system) is pulled longitudinally and then released. It gets displaced from equilibrium position. At any moment, forces applying on spring-mass system are:

(i) The restoring force: -kx where k is spring factor and

(ii) The damping force: -γv where v = dx/dt is instantaneous velocity of oscillator.

This signifies that for damped harmonic oscillator, equation of motion should comprise restoring force and the damping force. Therefore, in this case Eq. 2 is modified to

md

^{2}x/dt^{2}= -kx - γdx/dt.....................................Eq.2After rearranging terms and dividing throughout by m, equation of motion of the damped oscillator takes form

d

^{2}x/dt^{2}+ 2bdx/dt + ω_{0}^{2}x = 0.....................................Eq.3Where ω

_{0}^{2}= k/mConstant b has dimensions of force/(velocity x mass) = MLT

^{-2}/LT^{-1}= T^{-1}Therefore, its unit is s

^{-1}, which is same as that of ω_{0}.:Solutions of differential equationTo discover how damping influences amplitude of oscillation to solve Eq. when both the restoring force and the damping force are present. The general solution, as discussed above, should involve both exponential and harmonic terms. Let us therefore take a solution of the form

X(t) = a exp(αt) .....................................Eq.4

Where a and α

Differentiating Eq.4 twice with respect to time, we attain

dx/dt = aα exp(αt)

and d

^{2}x/dt^{2}= aα^{2}exp(αt)Substituting the expressions in Eq.3, we get

(α

^{2}+ 2bα + ω_{0}^{2})a exp(αt) = 0.....................................Eq.5For this equation to hold at all times, we should either have

a = 0

Which is trivial, or

α

^{2}+ 2bα + ω_{0}^{2}= 0.....................................Eq.6This equation is quadratic in α. Let us call the two roots α

_{1}and α_{2}α

_{1}= -b + (b^{2}- ω_{0}^{2})^{1/2}.....................................Eq.7(a)α

_{2}= -b - (b^{2}- ω_{0}^{2})^{1/2}.....................................Eq.7(b)Since Eq.3 is linear, principle .of superposition is appropriate. Therefore, general solution is get by the superposition of x

_{1}and x_{2}:x(t) = exp(-bt)[a

_{1}exp{(b^{2}- ω_{0}^{2})^{1/2}}t + a_{2}exp{exp{-(b^{2}- ω_{0}^{2})^{1/2}}].....................................Eq.8Heavy Damping:The system returns to equilibrium position very slowly, without any oscillation. Heavy damping takes place when resistive forces exceed those of critical damping.

When resistance to motion is very strong, system is said to be heavily damped. Springs joining wagons of the train constitute the most significant heavily damped system.

Mathematically, a system is said to be heavily damped if b > ω0. Then the quantity is positive definite. If we put

β = √b

^{2}- ω_{0}^{2}The general solution for damped oscillator given by Eq.9 reduces to

x(t) = exp(-bt)[a

_{1}exp(βt)[a_{1}exp(βt) + a_{2}exp(-βt)].....................................Eq.9This represents non-oscillatory behavior. Such a motion is called dead-beat. The actual displacement will, however, be determined by the initial conditions. Let us suppose that to begin with, the oscillator is at its equilibrium position, i.e x = 0 at t = 0. Then we give it a sudden kick so that it acquires a velocity v0, i.e. v = v0 at t = 0. Then from Eq.9 we have

a

_{1}+ a_{2}= 0 and -b(a_{1}+ a_{2}) +β(a_{1}- a_{2}) = v0These equation may be solved to give

a

_{1}= -a_{2}= v_{0}/2βOn substituting these results in Eq.9, we can write solution in compact form:

x(t) = v

_{0}/2βexp(-bt)[exp(βt)-exp(-βt)]= v

_{0}/βexp(-bt)sinhβt...............Eq.10Where sinhβt = (1/2)[exp(βt) - exp(-βt)] is hyperbolic sine function. From Eq. 10 it is clear that x(t) will be determined by interplay of the increasing hyperbolic function and decaying exponential. These are plotted separately. Eq.10 shows for heavily damped system when it is rapidly disturbed from equilibrium position. You will note that at first displacement increases with time. But soon exponential term becomes significant and displacement starts to decrease slowly.

Critical Damping:You may have observed that on hitting an isolated road bump, a car bounces up and down and the occupants feel uncomfortable. To minimise this discomfort, the bouncing caused by the road bumps must be damped very rapidly and the automobile restored to equilibrium quickly. For this we use critically damped shock absorbers. Critical damping is also useful in recording instruments such as a galvanometer (pointer type as well as suspended coil type) which experience sudden impulses. We require the pointer to move to the correct position in minimum time and stay there without executing oscillations. Similarly, a ballistic galvanometer coil is required to return to zero displacement immediately.

Mathematically, we say that a system is critically damped if b is equal to the natural frequency, ω

_{0}, of the system. This means that b^{2}- ω_{0}^{2}= 0, so that eq.8 reduces to x(t) = (a_{1}+ a_{2})exp(-βt)= aexp(-bt)

Where a = a

_{1}+ a_{2}So, two terms in Eq.8 give same time dependence and reduce to one term. It can be easily verified that in this case general solution of Eq.3

x(t) = (P + qt)exp(-bt)

Where p and q are constants, p has dimensions of length and q those of velocity. These can be find out easily from initial conditions.

Let us suppose that system is disturbed from mean equilibrium position by the sudden impulse. (Coil of a suspended type galvanometer receives some electric charge at t = 0.) That is, at t = 0, x(0) and dx/dt|

_{t=0}= v_{0}. This provides p = 0 and q = v_{0}, so that complete solution isx(t) = v

_{0}texp(-bt)Average energy of a weakly damped oscillator:The presence of damping the amplitude of oscillation decreases with the passage of time. This means that energy is dissipated in overcoming resistance to motion. From Unit I we recall that at any time, the total energy of a harmonic oscillator is made up of kinetic and potential components. We can still use the same definition and write

E(t) = K.E.(t) + U(t)

= 1/2m(dx/dt)

^{2}+ 1/2kx^{2}Where dx/dt denotes instantaneous velocity

For the weakly damped harmonic oscillator, instantaneous displacement is given by Equation:

x(t) = a

_{0}exp(-bt)cos(w_{d}t + Φ). By differentiating it with respect to time, instantaneous velocity is:dx(t)/dt = v = -a

_{0}exp(-bt)[bcos(ω_{d}t + Φ) + ω_{d}sin(ω_{d}t + Φ)]^{2}Therefore, kinetic energy of oscillator is

K.E. = 1/2ma

^{2}_{0}exp(-2bt)[b^{2}cos^{2}(ω^{d}t + Φ) + ω_{d}^{2}sin^{2}(ω_{d}t + Φ) + bω_{d}sin2(ωdt + Φ)]Potential energy of the oscillator is:

U = 1/2kx2 = 1/2mω02x2

Therefore, total energy of oscillator at any time t is provided by

E(t) = 1/2ma

^{2}_{0}exp(-2bt)[(b^{2}+ ω_{0}^{2})cos^{2}(ω_{d}t + Φ) + ω_{d}^{2}sin^{2}(ωdt + Φ) + bωdsin2(ωdt + Φ)]:Methods of characterizing damped systemsIn the viscous damping model, damped oscillator is characterized by γ and ω

_{0}. We also know that model applies to greatly different physical systems. In certain cases it is more suitable to employ other parameters to characterize damped motion.Logarithmic Decrement:The most suitable method to find out amount of damping present in the system is to estimate rate at which amplitude of oscillation dies away. Let us assume damped vibration shown graphically in Figure given below. Let a

_{0}and a_{1}be first two successive amplitudes of oscillation separated by one period.You will note that the amplitudes lie in same direction/quadrant. If T is period of oscillation, then using Equation for weakly damped oscillator, we can write

a1 = a

_{0}exp(-bt)So that a

_{0}/a_{1}= exp(bt) = exp(γT/2m)The ratio a

_{0}/a_{1}, larger amplitude is in numerator. That is why this ratio is known as decrement. It is denoted by the symbol d.Relaxation Time:In physics we estimate decay of the quantity in terms of fraction e

^{-1}of initial value. This provides us another way of defining damping effect by means of time taken by amplitude to decay to e^{-1}= 0.368 of original value. This time is known as relaxation time. To understand this, we recall that amplitude of damped oscillation is given bya(t) = a

_{0}exp(-bt)If we signify amplitude of oscillation after the interval of time r by a (t + T), we can write

a(t + τ) = a

_{0}exp(-b(t + T))By taking the ratio a(t + T)/a(t), we get

1/e for bT = 1

This illustrates that for b = T

^{-1}amplitude drops to 1/ e = 0.368 of initial value. Using result in Eq., we obtain:<P> = 2<E>/T

Relaxation time, T, is thus a measure of rapidity with which motion is damped.

Quality Factor:Yet another way of stating damping effect is by means of rate of decay of energy. Average energy of weakly damped oscillator decays to E

_{0}e^{-1}in time t = 1/2b = m/γ seconds. If ω_{d}is angular frequency, then in this time oscillator will vibrate through ω_{d}m/γ radians. Number of radians by which weakly damped system oscillates as its average energy decays to E_{0}e^{-1}is the measure of quality factor, Q:Q = ω

_{d}m/γ = ω_{d}/2b = ω_{d}T/2For weakly damped mechanical oscillator, quality factor can be stated in terms of spring factor and damping constant. For weak damping,

ω

_{d}≈ ω_{0}= √k/mTherefore Q = √km/γ

That is, quality factor of the weakly damped oscillator is directly proportional to square root of k and inversely proportional to γ

Examples of damped systems:All harmonic oscillators in nature have some damping, which in general, is quite small. To understand effect of damping, consider two specific cases: (a) Oscillations of charge in the LCR circuit, and (b) motion of coil in the suspension type galvanometer. These are of particular interest to us as former has wide applications in radio engineering and latter is utilized in physics laboratory.

An LCR Circuit:In the ideal LC circuit, charge excites SHM. If a current I flows through circuit because of discharging/charging of capacitor, voltage drop across resistor is RI. Therefore Equation modifies to -

q/c = -LdI/dt - dq/dt

Equation can be written as:

Ld

^{2}q/dt^{2}+ Rdq/dt + q/c = 0L, R and 1/C are respectively analogous to m, γ and k.

This signifies that resistor in electric circuit has exactly analogous effect as that of viscous force in the mechanical system. To continue further, divide Equation given above throughout by L attaining

d

^{2}q/dt^{2}+ (R/L)dq/dt + 1q/LC = 0ω

_{0}^{2}= 1/LC and b = R/2LWe know that b has dimensions of time inverse. This means that R/L has unit of s

^{-1}, same as that of ω_{0}. That is why ω_{0}L is measured in ohm.For the weakly damped circuit, charge on capacitor plates at time t is

q(t) = q

_{0}exp(-Rt/2L)cos(ω_{d}t + Φ) with angular frequencyωd = √(1/LC - R

^{2}/4L^{2})Equation given above illustrates that charge amplitude q

_{0}exp(-R/2Lt) will decay at rate that depends on resistance. Therefore in LCR circuit, resistance is only dissipative element; an increase in R increases rate of decay of charge and decreases frequency of oscillations.Since ω

_{0}L is measured in ohms, 1/ω_{0}C is also measured in ohms. These are respectively referred to as inductive reactance and capacitive reactance.Q value of weakly damped LCR circuit is

Q = ω

_{d}/2b ≈ ω_{0}L/R = 1/R√L/CThis equation shows that for a purely inductive circuit (R = 0), quality factor will be infinite.

Tutorsglobe: A way to secure high grade in your curriculum (Online Tutoring)Expand your confidence, grow study skills and improve your grades.

Since 2009, Tutorsglobe has proactively helped millions of students to get better grades in school, college or university and score well in competitive tests with live, one-on-one online tutoring.

Using an advanced developed tutoring system providing little or no wait time, the students are connected on-demand with a tutor at www.tutorsglobe.com. Students work one-on-one, in real-time with a tutor, communicating and studying using a virtual whiteboard technology. Scientific and mathematical notation, symbols, geometric figures, graphing and freehand drawing can be rendered quickly and easily in the advanced whiteboard.

Free to know our price and packages for online physics tutoring. Chat with us or submit request at info@tutorsglobe.com