Small oscillations on a constraint curve

tmode
Messages
3
Reaction score
0

Homework Statement



From Goldstein Classical Mechanics, 6.16:

A mass particle moves in a constant vertical gravitational field along the curve defined by y=ax4 , where y is the vertical direction. Find the equation of motion for small oscillations about the position of equilibrium.

The Attempt at a Solution



When taking a quadratic approximation as in the standard way to solve small oscillations, the eigenfrequency vanishes, giving a result without oscillation, which is clearly not what I'm looking for. I'm honestly not sure how to proceed with this.
 
Physics news on Phys.org
Hi tmode, welcome to PF. Show what you tried in detail, please. How did you take the constraint into account?

ehild
 
Hello, thank you.

I took the constraint into account by substituting y=ax^4 and \dot{y}=4ax^3 into the Lagrangian: L = T-V = \frac{1}{2}m(\dot{x}^2 + \dot{y}^2) - mgy
To yield, without any small oscillation approximation: L = \frac{1}{2}m\dot{x}^2 (1+16 a^2 x^6) -mgax^4
However, when you try to Taylor expand the m(1+16a^2 x^6) and mgax^4 about the equilibrium point x=0 (ie. small x), at quadratic order you get back m in the former and 0 in the latter, resulting in ω=0 for this treatment as a small oscillation approximation, and therefore no oscillation. For the potential, the first non-zero Taylor expansion approximation is the original value mgax^4, which would be very small in any case, on the order of x^4.

Taking the Euler-Lagrange equation to get equation of motion, we end up with 0=m\ddot{x}(1+16a^2 x^6) + 48m\dot{x}^2 a^2 x^5 + 4mga x^3
I don't think this is all the question wants answered, since it's not taking any sort of small oscillation approximation, nor does it give an explicit equation of motion.
 
Last edited:
The equation you derived is correct. I do not see anything you can do more. Even ignoring the terms with x^5 and x^6 it is not an equation for SHM.

ehild
 
The motion is not harmonic, why do you expect to find a equation for harmonic motion?

In first non-vanishing order, we get ##\ddot{x} = -4gax^3##, which is just the derivative of the potential => okay.
 
Alright, thanks both of you. In the context of the level of difficulty of the rest of the problems from that chapter, it seemed far too simple leaving it there, seeming like an incomplete solution. Thanks again.
 
To solve this, I first used the units to work out that a= m* a/m, i.e. t=z/λ. This would allow you to determine the time duration within an interval section by section and then add this to the previous ones to obtain the age of the respective layer. However, this would require a constant thickness per year for each interval. However, since this is most likely not the case, my next consideration was that the age must be the integral of a 1/λ(z) function, which I cannot model.
Back
Top