Classical Mechanics, cycloid pendulum

In summary, the cycloid pendulum is a system where a particle is under the effect of a constant gravitational field and moves without friction over a curve described parametrically by x=a(\theta + \sin \theta) and y=a(1-\cos \theta ). The Lagrangian and equations of motion can be derived using the generalized coordinate \theta, and for small angles, the motion equation becomes \ddot \theta +K \theta=0, where K is a constant. The solution to this differential equation is a periodic function with a period of \frac{2 \pi}{\sqrt {\frac{g}{2a+1}}}.
  • #1
fluidistic
Gold Member
3,949
264

Homework Statement


The cycloid pendulum consists of a particle under the effect of a constant gravitational field [itex](\vec g = -g \hat y)[/itex] that moves without friction over a curve described parametrically by [itex]x=a(\theta + \sin \theta)[/itex] and [itex]y=a(1-\cos \theta )[/itex].
1)Write down the Lagrangian and the equations of motion. (Hint, use [itex]\theta[/itex] as a generalized coordinate)
2)Solve the motion equation and verify that the pendulum is rigorously chronicle.

Homework Equations



L=T-V. Lagrange's equation.

The Attempt at a Solution


1)[itex]L=T-V. V=mgh=mgy[/itex].
[itex]v^2=2a^2 \dot \theta ^2 [1+\cos \theta ][/itex].
Thus the Lagrangian [itex]L=ma ^2 \dot \theta ^2 [1+ \cos \theta ] + mga [\cos \theta -1 ][/itex].
Using [itex]\theta[/itex] as the only generalized coordinate, I derive the equation of motion as being [itex]2a [\cos ( \theta )\ddot \theta + \ddot \theta - \sin (\theta )\dot \theta ^2 ]+\sin \theta (a \dot \theta ^2 +g )=0[/itex].
I must solve this... is this really serious? :/
If I consider small oscillations, [itex]\sin \theta \approx \theta[/itex] and [itex]\cos \theta \approx 1[/itex]. The motion equation reduces to [itex]\ddot \theta - \frac{\theta \dot \theta}{2}+ \frac{\theta \dot \theta ^2}{4}+\frac{g \theta }{4a}=0[/itex]. Which still doesn't look trivial to me to solve. :/

Did I go wrong somewhere? If not, how do I tackle the equation of motion? (and which one do I choose?)
Thanking you!
 
Physics news on Phys.org
  • #2
The reduced equation doesn't look really okay to me (check the dimensions of the 2nd and the 3rd terms). Anyway, in small angle approximation, we usually omit all the terms of high degrees of not only [itex]\theta[/itex] but also [itex]d\theta/dt[/itex] and [itex]d^2\theta/dt^2[/itex]. So the 2nd and 3rd terms should be omitted.

That fact is, I believe, for this reason: intuitively, if the distance is short, and you know the time taken is not so small, you should expect that the speed and the acceleration are also small. Usually this argument applies well to oscillation, as when we reduce the oscillation equation to the 1st degree (to the degree of linearity) - or harmonic oscillation, the "time taken" is specified by the period which doesn't depend on how short or long the amplitude (or "distance") is and instead depends on the choice of the inherent characteristics of the system. So normally we cannot assume that the period (or "time taken") is small.
 
  • #3
hikaru1221 said:
The reduced equation doesn't look really okay to me (check the dimensions of the 2nd and the 3rd terms). Anyway, in small angle approximation, we usually omit all the terms of high degrees of not only [itex]\theta[/itex] but also [itex]d\theta/dt[/itex] and [itex]d^2\theta/dt^2[/itex]. So the 2nd and 3rd terms should be omitted.

That fact is, I believe, for this reason: intuitively, if the distance is short, and you know the time taken is not so small, you should expect that the speed and the acceleration are also small. Usually this argument applies well to oscillation, as when we reduce the oscillation equation to the 1st degree (to the degree of linearity) - or harmonic oscillation, the "time taken" is specified by the period which doesn't depend on how short or long the amplitude (or "distance") is and instead depends on the choice of the inherent characteristics of the system. So normally we cannot assume that the period (or "time taken") is small.
Thanks a lot for your reply.
Well the dimensions will never match for each term I think. g is in m/s while the phi's are rad, rad/s, rad/s², etc. So basically my equations is all wrong :/
You said that I should depreciate the 2nd and 3rd term. But the first term is angular acceleration I think, which should be considered small too? I mean, why not depreciate the first term too?
 
  • #4
Then you have to omit the last term too :biggrin: It depends on the degree of approximation. The 2nd and 3rd terms are omitted because they are of higher degree (small quantity x small quantity; while for the 1st and the last ones, it's just one small quantity).
But you'd better check the equation first.
 
  • #5
hikaru1221 said:
Then you have to omit the last term too :biggrin: It depends on the degree of approximation. The 2nd and 3rd terms are omitted because they are of higher degree (small quantity x small quantity; while for the 1st and the last ones, it's just one small quantity).
But you'd better check the equation first.

Ahh thanks a lot, now I think I understand what you meant. Small quantity x small quantitie = small quantity squared = very small quantity= negligible quantity.
Ok I'll redo the algebra...
 
  • #6
I've redone the algebra and I find a simpler expression, so I'm very grateful to you! :D
Equation of motion: [itex]\ddot \theta [2a + \cos \theta ]+g \sin \theta =0[/itex].
For small angles this becomes [itex]\ddot \theta (2a+1) +g \theta =0[/itex]. Or written even simpler, [itex]\ddot \theta +K \theta=0[/itex].
I hope I'm right this time...

Edit: Solving this diff. eq. I get [itex]\theta (t)=c_1 \sin (\sqrt K t)+c_2 \cos (\sqrt K t)[/itex]. It's obvious to me now that it's a periodic function. I'll try to find the period.

Edit 2: Done, period is [itex]\frac{2 \pi}{\sqrt {\frac{g}{2a+1}}}[/itex]. I'd love if someone could verify my answer. (It's not even homework, I'm preparing myself for the final exam on mid/end July. I just want to be sure I'm doing things the right way :) )
 
Last edited:
  • #7
Check the dimension of the equation of motion :smile:
I guess it should be [itex]\frac{d^2\theta}{dt^2}2a(1+cos\theta)+gsin\theta = 0[/itex]
 
  • #8
fluidistic said:
Equation of motion: [itex]\ddot \theta [2a + \cos \theta ]+g \sin \theta =0[/itex].

Are you sure? I got [itex]2a\ddot \theta [1 + \cos \theta ]+g \sin \theta =0[/itex]

Edit: Solving this diff. eq. I get [itex]\theta (t)=c_1 \sin (\sqrt K t)+c_2 \cos (\sqrt K t)[/itex]. It's obvious to me now that it's a periodic function. I'll try to find the period.

Yes, that's right, assuming you use the right value of K. Typically, solutions of differential equations that look like [itex]\ddot x + w^2x =0[/itex] are written as Asin(wx+phi). This is mathematically identical to your solution, but it's easier to understand intuitively.
 
  • #9
fluidistic said:
Using [itex]\theta[/itex] as the only generalized coordinate, I derive the equation of motion as being [itex]2a [\cos ( \theta )\ddot \theta + \ddot \theta - \sin (\theta )\dot \theta ^2 ]+\sin \theta (a \dot \theta ^2 +g )=0[/itex].

It is right, but you made a mistake in the next line, it should be:

fluidistic said:
[itex]\ddot \theta - \frac{\theta \dot \theta^2}{2}+ \frac{\theta \dot \theta ^2}{4}+\frac{g \theta }{4a}=0[/itex].

ehild
 
  • #10
Thanks guys for your help. I've redone the algebra 3 times now, starting from the Lagrangian. My equation of motion is [itex]2a \ddot \theta [1+ \cos \theta]-a\dot \theta ^2 \sin \theta +g \sin \theta =0[/itex].
It seems like I have an extra "[itex]-a \dot \theta ^2 \sin \theta[/itex]" term. I don't know how you could have got rid of it.
[itex]\frac{\partial L}{\partial \theta}=-ma\sin \theta [a \dot \theta ^2 -g][/itex].
While [itex]\frac{\partial L}{\partial \dot \theta}=2ma^2 [1+\cos \theta ]\dot \theta[/itex].
Deriving this with respect to time gave me[itex]\frac{d}{dt} \left ( \frac{\partial L}{\partial \dot \theta} \right ) =2ma^2 \{ \ddot \theta [1+\cos \theta ] - \dot \theta ^2 \sin \theta \}[/itex]. Hence my equation of motion... do you spot any error in my derivation?
And I don't really know what is meant by "check out the dimension of the equation of motion"
Some term has rad/s², another has m/s² because of the g term so the dimension seems like ill defined in my (our?!) equation of motion.
Edit: For small oscillations, our motion equation is the same, namely [itex]\ddot \phi+\frac{g \phi}{4a}=0[/itex]. Solving it gives [itex]\theta (t)=A\sin \left ( \frac{1}{2} \sqrt{\frac{g}{a}} t + \phi \right )[/itex] with period [itex]\pi \sqrt{ \frac{a}{g}}[/itex] if I didn't make any error. So yes, the pendulum has a period, it's chronicle.
 
Last edited:
  • #11
fluidistic said:
It seems like I have an extra "[itex]-a \dot \theta ^2 \sin \theta[/itex]" term. I don't know how you could have got rid of it.

You can assume it's negligible because it's the product of sin(theta), a small quantity, with [itex]-a \dot \theta ^2[/itex], an even smaller quantity.

[itex]\frac{\partial L}{\partial \theta}=-ma\sin \theta [a \dot \theta ^2 -g][/itex].
While [itex]\frac{\partial L}{\partial \dot \theta}=2ma^2 [1+\cos \theta ]\dot \theta[/itex].
Deriving this with respect to time gave me[itex]\frac{d}{dt} \left ( \frac{\partial L}{\partial \dot \theta} \right ) =2ma^2 \{ \ddot \theta [1+\cos \theta ] - \dot \theta ^2 \sin \theta \}[/itex]. Hence my equation of motion... do you spot any error in my derivation?

Nope, it's correct.

And I don't really know what is meant by "check out the dimension of the equation of motion"
Some term has rad/s², another has m/s² because of the g term so the dimension seems like ill defined in my (our?!) equation of motion.

This isn't possible. In any equation, all of the terms must have exactly the same units, or else the equation is guaranteed to be wrong. However, your equation of motion has consistent units, because a^2 has units of m^2.

Edit: For small oscillations, our motion equation is the same, namely [itex]\ddot \phi+\frac{g \phi}{4a}=0[/itex]. Solving it gives [itex]\theta (t)=A\sin \left ( \frac{1}{2} \sqrt{\frac{g}{a}} t + \phi \right )[/itex] with period [itex]\pi \sqrt{ \frac{a}{g}}[/itex] if I didn't make any error. So yes, the pendulum has a period, it's chronicle.

I think the period should have an extra factor of 4 in front of it. Otherwise, your answer is correct.
 
  • #12
Thank you very much for all guys. I didn't know a had units. Now everything is clear.
Problem solved.
 

FAQ: Classical Mechanics, cycloid pendulum

What is classical mechanics?

Classical mechanics is a branch of physics that studies the motion of objects and systems under the influence of forces. It is based on Isaac Newton's laws of motion and includes topics such as kinematics, dynamics, and energy.

What is a cycloid pendulum?

A cycloid pendulum is a type of pendulum where the bob (weight) moves along a cycloid curve instead of a simple arc. This results in a more complex motion and can be used to demonstrate principles of classical mechanics.

What factors affect the period of a cycloid pendulum?

The period of a cycloid pendulum is affected by the length of the pendulum, the mass of the bob, and the amplitude of the cycloid curve. It is also affected by the acceleration due to gravity and any external forces acting on the pendulum.

How is energy conserved in a cycloid pendulum?

In a cycloid pendulum, energy is conserved because the potential energy of the bob is converted into kinetic energy as it moves along the cycloid curve. This kinetic energy is then converted back into potential energy as the bob swings back and forth.

What are some real-world applications of cycloid pendulums?

Cycloid pendulums have been used in various applications, including as timekeeping devices in clocks and as demonstration tools in physics education. They are also used in seismographs to measure earthquake activity and in some amusement park rides.

Back
Top