Sunday, January 24, 2016

homework and exercises - Nonlinear spring $F=-kx^3$


A nonlinear spring whose restoring force is given by $F=-kx^3$ where $x$ is the displacement from equilibrium , is stretched a distance $A$. Attached to its end is a mass $m$. Calculate....(I can do that) ..suppose the amplitude of oscillation is increased, what happens to the period?


Here's what I think: If the amplitude is increased the spring posses more total energy, at equilibrium the spring is traveling faster than before because it posses more kinetic energy. I think in the spring travels faster when it's at a similar displacement from equilibrium, but it has to travel more distance, so I can't conclude anything.


I was think about solving,


$$mx''=-kx^3$$


But realized this is a very hard job.


Any ideas? Thanks.



Answer



The potential energy is $U\left(x\right) = kx^4/4$ since $-d/dx\left(kx^4/4\right) = -kx^3 = F$, and the energy $$ E = \frac{1}{2}m\left(\frac{dx}{dt}\right)^2 + \frac{1}{4}kx^4 $$ is conserved.



From the above you can show that $$ \begin{eqnarray} dt &=& \pm \ dx \sqrt{\frac{m}{2E}}\left(1-\frac{k}{4E}x^4\right)^{-1/2} \\ &=& \pm \ dx \sqrt{\frac{2m}{k}} \ A^{-2} \left[1-\left(\frac{x}{A}\right)^4\right]^{-1/2} \end{eqnarray} $$ where the amplitude $A = \left(4E / k\right)^{1/4}$ can be found from setting $dx/dt = 0$ in the expression for the energy and solving for $x$.


The period is then $$ \begin{eqnarray} T &=& 4 \sqrt{\frac{2m}{k}} \ A^{-2} \int_0^A dx \left[1-\left(\frac{x}{A}\right)^4\right]^{-1/2} \\ &=& 4 \sqrt{\frac{2m}{k}} \ A^{-1} \int_0^1 du \left(1-u^4\right)^{-1/2} \\ &=& \left(4 \sqrt{\frac{2m}{k}} I\right) A^{-1} \\ &\propto& A^{-1} \end{eqnarray} $$ where $u = x/A$ and $I = \int_0^1 du \left(1-u^4\right)^{-1/2} \approx 1.31$ (see this).


You can repeat the above for a more general potential energy $U\left(x\right) = \alpha \left|x\right|^n$, where you should find that


$$ dt = \pm \ dx \sqrt{\frac{m}{2\alpha}} \ A^{-n/2} \left[1-\left(\frac{\left|x\right|}{A}\right)^n\right]^{-1/2} $$


and


$$ \begin{eqnarray} T_n &=& \left(4 \sqrt{\frac{m}{2\alpha}} I_n\right) A^{1-n/2} \\ &\propto& A^{1-n/2} \end{eqnarray} $$


where


$$ I_n = \int_0^1 du \left(1-u^n\right)^{-1/2} $$


can be evaluated in terms of gamma functions (see this).


This is in agreement with the above for $\alpha = k/4$ and $n=4$, and with Landau and Lifshitz's Mechanics problem 2a of section 12 (page 27), where they find that $T_n \propto E^{1/n-1/2} \propto A^{1-n/2}$.



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 \...