Tuesday, April 7, 2015

harmonic oscillator - How to determine phase angle for a sinusoidal motion?


If I have an over-damped mechanical system that is excited with a sinusoidal motion. That sinusoidal motion starts with a determined frequency then increases frequency over time. Of course, it is known that there will be a phase shift between the driving force and the motion of the hanging mass.


My question is, how to figure out phase lag of mass motion in relation to driving force?



Answer



A damped harmonic oscillator with a sinusoidal driving force is represented by the equation


$$\ddot{x} + \gamma\dot{x} + \omega_0^2x = \frac{F_D \sin(\omega_D t)}{m}$$


where $\gamma = b/m$ ($b$ is the damping coefficient, $b=F/v$) and $\omega_0^2 = k/m$ is the resonant frequency of the oscillator. The particular solution to this equation can be determined by taking the imaginary part of the solution to


$$\ddot{x} + \gamma\dot{x} + \omega_0^2x = \frac{F_D}{m}e^{i\omega_D t}$$



If you assume* the solution takes the form


$$x(t) = A e^{i(\omega_D t + \phi)}$$


and plug that in, you get


$$-A \omega_D^2 + \omega_0^2 A = \frac{F_D}{m}\cos(\phi)$$


and


$$\gamma\omega_D A = \frac{F_D}{m}\sin(\phi)$$


Solving for the phase difference gives


$$\tan\phi = \frac{\gamma\omega_D}{\omega_0^2 - \omega_D^2}$$


This depends on the frequency of the driving force and the resonant frequency of the oscillator, but not on the amplitude of the driving force.


You can express this in terms of the dimensionless variable $x = \omega_D / \omega_0$ as



$$\tan\phi = \frac{\gamma}{\omega_0}\frac{x}{1 - x^2}$$


and if you graph it,


plot of phi versus x


(graph generated by Wolfram Alpha) you'll see how the response of the oscillator jumps from leading to lagging when $\omega_D = \omega_0$ (at $x=1$), that is, when the driving and resonant frequencies are equal.




*The same solution can be obtained from Fourier decomposition without making this assumption.


No comments:

Post a Comment

classical mechanics - Moment of a force about a given axis (Torque) - Scalar or vectorial?

I am studying Statics and saw that: The moment of a force about a given axis (or Torque) is defined by the equation: $M_X = (\vec r \times \...