Finding the period of a limit cycle

Click For Summary
SUMMARY

The discussion focuses on determining the period of a limit cycle for the Duffing oscillator described by the equation \(\ddot x + x + \epsilon x^3 = 0\). The conserved quantity is identified as \(C = \frac{1}{2}x^2 + \frac{\epsilon}{4}x^4 + \frac{1}{2}\dot{x}^2\). The period \(T\) is expressed as an integral \(T = 2\int_0^{X} \frac{dx}{\sqrt{C - x^2 - \frac{\epsilon}{2}x^4}}\), where \(X\) is the maximum displacement. The main challenge discussed is selecting the appropriate bounds for the integral.

PREREQUISITES
  • Understanding of the Duffing oscillator dynamics
  • Familiarity with conserved quantities in mechanical systems
  • Knowledge of integral calculus, specifically definite integrals
  • Experience with power series expansions and their applications
NEXT STEPS
  • Study the derivation of conserved quantities in nonlinear oscillators
  • Learn about the method of parametrization for closed line integrals
  • Explore power series expansions in the context of differential equations
  • Investigate numerical methods for evaluating integrals with variable bounds
USEFUL FOR

Students and researchers in physics, particularly those studying nonlinear dynamics, mechanical systems, and oscillatory motion, will benefit from this discussion.

alexfloo
Messages
192
Reaction score
0

Homework Statement


We're considering the following special case of the Duffing oscillator, in the limit of small nonlinearity:

[itex]\ddot x + x + \epsilon x^3 = 0[/itex].

The problem asks us to find a conserved quantity, and use it to write the period as an integral. Then, to write it as a power series in [itex]\epsilon[/itex] and evaluate it, then check the first two terms agains a power series that was given in our text. Mine didn't match up. I'm pretty sure I got the integrand right, but I really just don't know how select the bounds.


Homework Equations





The Attempt at a Solution


First of all, this system resembles a force on a unit mass which doesn't depend on velocity. Therefore, there is a well-defined potential energy as a function of only x, which in this case is

[itex]\frac{1}{2}x^2+\frac{\epsilon}{4}x^4[/itex]

Therefore,

[itex]C=\frac{1}{2}x^2+\frac{\epsilon}{4}x^4+\frac{1}{2}\dot{x}^2[/itex]

is conserved. (We can check that will time derivatives but I won't.)

Therefore, along a particular trajectory, there is some constant C such that

[itex]\frac {dx}{dt}=\dot x=\sqrt{C-x^2-\frac{\epsilon}{2}x^4}[/itex].

Therefore

[itex]dt=\frac{dx}{\sqrt{C-x^2-\frac{\epsilon}{2}x^4}}[/itex].

So the period will be

[itex]T=\int \frac{dx}{\sqrt{C-x^2-\frac{\epsilon}{2}x^4}}[/itex].

The problem is I have no idea what to set the bounds of the integral as. The x-value at the start and end of the period are the same, so the integral will equal zero.

I recognize also that it'd make more sense to cast it as a closed line integral on the set

[itex]y=\sqrt{C-x^2-\frac{\epsilon}{2}x^4}[/itex],

but I'm not sure how to effectively parametrize that set, or else I should be writing it in a very different form.
 
Physics news on Phys.org
alexfloo said:

Homework Statement


We're considering the following special case of the Duffing oscillator, in the limit of small nonlinearity:

[itex]\ddot x + x + \epsilon x^3 = 0[/itex].

The problem asks us to find a conserved quantity, and use it to write the period as an integral. Then, to write it as a power series in [itex]\epsilon[/itex] and evaluate it, then check the first two terms agains a power series that was given in our text. Mine didn't match up. I'm pretty sure I got the integrand right, but I really just don't know how select the bounds.


Homework Equations





The Attempt at a Solution


First of all, this system resembles a force on a unit mass which doesn't depend on velocity. Therefore, there is a well-defined potential energy as a function of only x, which in this case is

[itex]\frac{1}{2}x^2+\frac{\epsilon}{4}x^4[/itex]

Therefore,

[itex]C=\frac{1}{2}x^2+\frac{\epsilon}{4}x^4+\frac{1}{2}\dot{x}^2[/itex]

is conserved. (We can check that will time derivatives but I won't.)

Therefore, along a particular trajectory, there is some constant C such that

[itex]\frac {dx}{dt}=\dot x=\sqrt{C-x^2-\frac{\epsilon}{2}x^4}[/itex].

Therefore

[itex]dt=\frac{dx}{\sqrt{C-x^2-\frac{\epsilon}{2}x^4}}[/itex].

So the period will be

[itex]T=\int \frac{dx}{\sqrt{C-x^2-\frac{\epsilon}{2}x^4}}[/itex].

The problem is I have no idea what to set the bounds of the integral as. The x-value at the start and end of the period are the same, so the integral will equal zero.

I recognize also that it'd make more sense to cast it as a closed line integral on the set

[itex]y=\sqrt{C-x^2-\frac{\epsilon}{2}x^4}[/itex],

but I'm not sure how to effectively parametrize that set, or else I should be writing it in a very different form.

Remember that when you take the square root there are two solutions to consider. The integral isn't zero because you have to choose the plus sign along one half cycle and the minus sign on the other half cycle. The period is thus more like

[tex]T = 2\int_0^{X} \frac{dx}{\sqrt{C - x^2 - \epsilon x^4/2}},[/tex]

where X is the maximum value x attains.
 

Similar threads

Replies
4
Views
2K
  • · Replies 31 ·
2
Replies
31
Views
3K
  • · Replies 3 ·
Replies
3
Views
2K
Replies
5
Views
3K
  • · Replies 13 ·
Replies
13
Views
2K
Replies
2
Views
2K
  • · Replies 6 ·
Replies
6
Views
2K
  • · Replies 10 ·
Replies
10
Views
3K
  • · Replies 9 ·
Replies
9
Views
4K
Replies
6
Views
3K