Harmonic oscillator with slight non-linearity

  • Thread starter groinsmash
  • Start date
  • #1
I have an interesting problem I have come across in my research. It results in the differential equation as follows:

[itex]x''+2γ(x')^\nu+\omega_{o}^2x=g(t)[/itex]

Primes indicate the derivative with respect to [itex]t[/itex]. [itex]\gamma[/itex] and [itex]\omega[/itex] are constants. The non-linearity comes from the first derivative [itex]x'[/itex] which is raised to the power of [itex]\nu[/itex]. [itex]\nu[/itex] is known to be 0.12 but can be between 0 and 1. The cases where [itex]\nu=0[/itex] or [itex]\nu=1[/itex] are easy enough. But how to go about tackling an arbitrary [itex]\nu[/itex]?

The problem may be made easier by noting that [itex]g(t)=1[/itex] for [itex]t\geq0[/itex] and [itex]0[/itex] for [itex]t<0[/itex].

Any ideas on how to go about solving it? Numerically or analytically (which would be amazing).

Thanks!
 

Answers and Replies

  • #2
1,800
53
You gotta' get Mathematica. Here's some Mathematica code:

Code:
\[Gamma] = 0.25
\[Nu] = 0.5; 
\[Omega] = 0.1
mysol = NDSolve[{Derivative[2][x][t] + 
      2*\[Gamma]*Derivative[1][x][t]^\[Nu] + \[Omega]*x[t] == 
     1, x[0] == 1, Derivative[1][x][0] == 1}, x, 
   {t, 0, 5}]
Plot[x[t] /. mysol, {t, 0, 5}]
 

Related Threads on Harmonic oscillator with slight non-linearity

Replies
1
Views
3K
Replies
3
Views
3K
  • Last Post
Replies
4
Views
9K
Replies
0
Views
1K
  • Last Post
Replies
1
Views
1K
Replies
11
Views
830
Replies
2
Views
3K
  • Last Post
Replies
2
Views
1K
Replies
3
Views
2K
  • Last Post
Replies
4
Views
1K
Top