Period of nonlinear spring-mass system


by jinteni
Tags: nonlinear, period, springmass
jinteni
jinteni is offline
#1
Nov28-07, 03:47 PM
P: 2
I wasn't sure where to post this problem, as it's physics related, but rather advanced in its math content (and it's a problem for my applied math course).

1. The problem statement, all variables and given/known data
Considering a spring-mass system (like here), given that the nonlinear spring has force = qx^3, where q is the spring stiffness, what is the period of the oscillation when the mass is released from rest at x_0?

3. The attempt at a solution

The equation of motion of the system is F=ma=-qx^3, so m*((d^2)x/dt^2)+q*x^3=0. Integrating, I get the energy of the system as (1/2)*m*((dx/dt)^2)+(1/4)*q*x^4.
When released from rest at x_0, the system then has E=(1/4)*q*x_0^4. I haven't found examples of this type of problem anywhere, so any help would be greatly appreciated!
Phys.Org News Partner Science news on Phys.org
SensaBubble: It's a bubble, but not as we know it (w/ video)
The hemihelix: Scientists discover a new shape using rubber bands (w/ video)
Microbes provide insights into evolution of human language
rock.freak667
rock.freak667 is offline
#2
Nov28-07, 04:34 PM
HW Helper
P: 6,214
[tex]m\frac{d^2x}{dt^2}-qx^3=0[/tex]

you simply can't integrate both sides

you would have to form the auxiliary eq'n and use the general solution for complex roots.
jinteni
jinteni is offline
#3
Nov28-07, 05:25 PM
P: 2
I've been shown an example of a nonlinear pendulum model, which gives an answer of something involving the EllipticK integral in Maple, but no information on how to derive that answer. For the linear models, the pendulum and spring-mass are directly related, so I don't know if this helps me in any way.

HallsofIvy
HallsofIvy is offline
#4
Nov28-07, 05:31 PM
Math
Emeritus
Sci Advisor
Thanks
PF Gold
P: 38,904

Period of nonlinear spring-mass system


rock.freak667's point is that you integrated [itex]m d^2x/dt^2[/itex] with respect to t and [itex]-qx^3[/itex] with respect to x! That's what you can't do!

since the independent variable, t, does not appear explicitely, you could use "quadrature". Let v= dx/dt. Then [itex]d^2x/dt^2= dv/dt= (dv/dx)(dx/dt)= v dv/dx[/itex]. The equation reduces to mvdv/dx- qx^3= 0. Now you can integrate with respect to x: [itex](1/2)mv^2- (q/4)x^4[/itex]= Constant. (If you look closely at that, you will see "conservation of energy".) Since v= dx/dt, that is the same as
[tex]dx/dt= \sqrt{\frac{g}{m}}\frac{x^2}{2}[/tex]
That should be easy to integrate. However, I'm not sure that that will answer your question about the period- in general non-linear spring motion is NOT periodic!
14c4b06b82
14c4b06b82 is offline
#5
Jul1-11, 12:35 AM
P: 1
x=displacement
a=amplitude
q=spring constant
m=mass
p=period

The energy stored by the spring for a given displacement is given by the integral of the force over the displacement.

[itex]E_{spring}=\int_{0}^{x} q*x^3 dx [/itex]

which evaluates to

[itex]E_{spring}=q*x^4/4[/itex]

when the mass is at it's maximum displacement and the mass is not moving all the energy in the system is stored in the spring. Thus the total energy in the system is the energy held by the spring when the displacement equals the amplitude.

[itex]E_{system}=q*x^4/4 , x = a[/itex]
[itex]E_{system}=q*(a)^4/4[/itex]
[itex]E_{system}=q/4*a^4[/itex]

applying conservation of energy, any energy lost by the spring must go to the mass as kinetic energy.

The kinetic energy of the mass is given by

[itex]E_{kinetic of mass}=\frac{1}{2}mv^2[/itex]

the energy held by the mass can then be expressed as

[itex]E_{kinetic of mass}=E_{system}-E_{spring}[/itex]

after substituting this gives

[itex]E_{kinetic of mass}=(q/4*a^4)-(q/4*x^4)[/itex]
[itex]E_{kinetic of mass}=\frac{q}{4}*(a^4-x^4)[/itex]

substituting this back into the expression for the kinetic energy of the mass

[itex]\frac{q}{4}*(a^4-x^4)=\frac{1}{2}mv^2[/itex]

solving for velocity gives

[itex]v_{mass}=\sqrt{\frac{q}{2m}*(a^4-x^4))}[/itex]

The mass could of course be traveling in either direction.

recall that

[itex]v=\frac{dx}{dt}[/itex]

thus

[itex]\frac{1}{v}=\frac{dt}{dx}[/itex]

in other words the reciprocal of velocity is the time passed per unit of distance.

applying this to the expression for the velocity of the mass gives

[itex]\frac{1}{v}=\sqrt{\frac{2m}{q*(a^4-x^4)}}[/itex]

or equivalently

[itex]\frac{dt}{dx}=\sqrt{\frac{2m}{q*(a^4-x^4)}}[/itex]

integrating both sides dx gives the time taken to travel between the start and end points of that integration. Integrating between [itex]0[/itex] and [itex]a[/itex] provides the time for a quarter period.

[itex]\int_{0}^{a}\frac{dt}{dx}dx=\int_{0}^{a}\sqrt{\frac{2m}{q*(a^4-x^4)}}dx[/itex]

(for some reason MathJax isn't rendering the above expression right. I hope it works for everyone else)

[itex]\int_{0}^{a}\sqrt{\frac{2m}{q*(a^4-x^4)}}dx[/itex]

[itex]=\sqrt{\frac{2m}{q}}*\int_{0}^{a}\sqrt{\frac{1}{a^4-x^4}}dx[/itex]

I'm not sure exactly how to justify this step but nvertheless it should work

[itex]=\sqrt{\frac{2m}{q}}*\frac{1}{a}*\int_{0}^{1}\sqrt{\frac{1}{1-x^4}}dx[/itex]

I have no idea how to approach that definite integral analytically but Wolfram can (definite integral~=1.31103)

working from there

[itex]p/4=\sqrt{\frac{2m}{q}}*\frac{1}{a}*1.33103[/itex]

[itex]p=\frac{7.41630}{a}*\sqrt{\frac{m}{q}}[/itex]
rock.freak667
rock.freak667 is offline
#6
Jul1-11, 01:27 AM
HW Helper
P: 6,214
You sort of bumped a 4 year old thread there .


Register to reply

Related Discussions
SHM, Determing Spring Constant from Period^2 vs. Mass Introductory Physics Homework 1
period of mass M hanging vertically from a spring Introductory Physics Homework 2
Conceptual - Mass, Period, Spring Constant Introductory Physics Homework 1
Period of a Mass on a Spring in Simple Harmonic Motion Introductory Physics Homework 7
Period of a mass on a spring Introductory Physics Homework 5