Can the Inverse Tan Function Justify Phase Lag in a Driven Harmonic Oscillator?

In summary, the equation for delta phase lag of a driven oscillator is W = Wo, with a phase lag of pi radians when W is much larger than Wo.
  • #1
Glenn G
113
12
TL;DR Summary
Confused by the Pi rads phase lag
1666170306103.png
This is an equation I found for the delta phase lag of a driven oscillator. W is the driving angular frequency and Wo is the natural angular frequency of the driven system. Of course this is just a small part of the solution to the differential equation.
Now ... 1) when W is much smaller than Wo the delta phase lag is zero (the driven moves in synch with the driver
2) when W = Wo we have a resonance condition and the phase lag is pi/2
3) when W is much larger than Wo there is a phase lag of pi radians
This can be observed nicely with Barton's pendulae.

The shape of the graph of delta phase lag against w that you find looks like ...
1666170733572.png

But if you plot a graph of a made up function such as ...
1666170813731.png

It seems that if you move the negative portion up you get the correct kind of shape but my issue is that inverse tan only goes from -pi/2 to +pi/2 so the inverse tan of anything, whatever it is never approaches pi.

How can this approach be justified? because it almost seems to me that from the maths that if W is very much bigger than Wo then the arctan is negative so rather that a phase lag, if you are changing sign doesn't this suggest a phase lead? Clearly this is not what happens though when you watch the system oscillating and also it can't be anyway, doesn't make sense.

I would love your input on this ...
 
Physics news on Phys.org
  • #2
This is the usual trouble with this formula to get phase shifts. You shouldn't use the naive arctan formula. It's the same as finding the angle of a point ##(x,y)## (Cartesian coordinates) in polar coordinates. The angle is uniquely defined by the two equations ##\cos \varphi=x/R## and ##\sin \varphi=y/R## with ##R=\sqrt{x^2+y^2}##. Then it's easy to see that the uniquely defined angle is given by
$$\varphi=\begin{cases} \text{sign} y \arccos\left (\frac{x}{R} \right ) & \text{for} \quad y \neq 0, \\ 0 & \text{for} \quad y=0, \quad x>0, \\
\pi &\text{for} \quad y=0, \quad x<0. \end{cases}$$
This angle is in the range ##\varphi \in (-\pi,\pi]##.
 
  • Like
Likes sophiecentaur, DaveE, topsquark and 1 other person
  • #3
Thanks Vanhees71 …. So is there a better ‘match’ of formula for the actual shape of the phase characteristics that does go from 0 to pi as W goes from 0 to infinite as in my sketch if the normal arctan isn’t appropriate ?
 
  • #4
My convention for the equation of motion is
$$\ddot{x}+2 \gamma \dot{x} + \omega_0^2 x=A \cos(\Omega t)$$
Then amplitude and phase shift for the stationary state are
$$\hat{x}=\frac{A}{\sqrt{(\omega_0^2-\Omega^2)^2+4 \gamma^2 \Omega^2}},$$
$$\varphi=+\arccos \left (\frac{\omega_0^2-\Omega^2}{\sqrt{(\omega_0^2-\Omega^2)^2+4 \gamma^2 \Omega^2}} \right ).$$
This is of course only valid for the case of weak damping, i.e., ##\omega_0>\gamma##.
[EDIT] This is, of course, nonsense. It's valid for any case of damping!
 
Last edited:
  • Like
Likes topsquark and Glenn G
  • #5
Yay, that looks great for giving the 0, pi/2 and pi relationship for the correct ‘zones’ - was Barton’s pendulae that got me thinking about this. Thanks.
 
  • Like
Likes vanhees71
  • #6
The relation ## \cot x = 1 / \tan x ## translates to ## \text{arccot}~x = \arctan \frac 1 x ##, and ## \cot x = \tan(\frac \pi 2 - x) ## to ## \text{arccot}~x = \frac \pi 2 - \arctan x ##. So you can use ## \arctan x = \frac \pi 2 - \arctan \frac 1 x ## to rewrite the original puzzling equation as $$
\delta = \frac \pi 2 + \arctan {\omega^2 - \omega_0^2 \over 2 \gamma \omega} \ .
$$ With multi-valued functions one should be careful not to cross branch cuts inadvertently.
 
  • #7
arctan is never a good idea when it comes to find the polar angle in polar coordinates. For that reason most computer languages have a function like atan2 (in Fortran it's ##\phi=\text{atan2}(y,x)## leading to values between ##-\pi## and ##\pi## and with specific rules for ##y=0##).
 
  • #8
vanhees71 said:
arctan is never a good idea
Never? The rewritten formula works just fine. :-)
 
  • Like
Likes vanhees71
  • #9
It's right, because in this case the imaginary part is always positive ;-).
 

Related to Can the Inverse Tan Function Justify Phase Lag in a Driven Harmonic Oscillator?

1. What is a driven harmonic oscillator?

A driven harmonic oscillator is a physical system that exhibits oscillatory motion when subjected to an external force or driving force. It can be described mathematically using the equation of motion for a simple harmonic oscillator, but with an additional term representing the driving force.

2. What are some examples of driven harmonic oscillators?

Some examples of driven harmonic oscillators include a mass-spring system with an external force applied, a pendulum with an external driving force, and an electric circuit with an alternating current source.

3. How does the driving frequency affect a driven harmonic oscillator?

The driving frequency of a driven harmonic oscillator determines the amplitude and phase of its oscillations. When the driving frequency matches the natural frequency of the system, resonance occurs and the amplitude of the oscillations becomes very large. When the driving frequency is different from the natural frequency, the amplitude of the oscillations is smaller.

4. What is the difference between a damped and undamped driven harmonic oscillator?

A damped driven harmonic oscillator experiences a damping force that decreases its amplitude over time, while an undamped driven harmonic oscillator does not experience any damping and its amplitude remains constant. Damping can be caused by friction, air resistance, or other dissipative forces.

5. How is a driven harmonic oscillator used in real-world applications?

Driven harmonic oscillators have many practical applications, such as in electrical circuits, mechanical systems, and musical instruments. They are also used in scientific research to study the effects of external forces on oscillatory systems and to model complex physical phenomena.

Similar threads

Replies
7
Views
710
  • Introductory Physics Homework Help
Replies
7
Views
2K
  • Classical Physics
Replies
17
Views
1K
Replies
3
Views
7K
  • Classical Physics
Replies
12
Views
7K
Replies
17
Views
5K
  • Mechanical Engineering
Replies
2
Views
157
Replies
3
Views
1K
  • Introductory Physics Homework Help
Replies
8
Views
211
Back
Top