Equations of motion of a particle over a cone+conserved quantities

Click For Summary

Homework Help Overview

The discussion revolves around the equations of motion for a particle constrained to move on the surface of a cone, utilizing Lagrangian mechanics. The original poster draws parallels to a previous problem involving a cylinder and explores the implications of conservation laws in this new context.

Discussion Character

  • Exploratory, Conceptual clarification, Mathematical reasoning, Assumption checking

Approaches and Questions Raised

  • The original poster attempts to derive the equations of motion using cylindrical coordinates and expresses concerns about the conservation of momentum in relation to the constraints imposed by the geometry of the cone. Participants question the independence of variables in the Lagrangian and discuss methods for incorporating constraints.

Discussion Status

Participants are actively engaging with the mathematical formulation of the problem, offering corrections and alternative approaches. There is recognition of the need for additional equations to account for the dependencies between variables, and some guidance is provided on how to express relationships between the coordinates.

Contextual Notes

There is an ongoing discussion about the implications of substituting one variable for another in the Lagrangian and how this affects the conservation of momentum. Participants also reflect on the appropriateness of using certain coordinates in the context of Lagrangian mechanics.

fluidistic
Gold Member
Messages
3,934
Reaction score
286

Homework Statement


Exactly the same problem as https://www.physicsforums.com/showthread.php?p=3335113#post3335113 but instead of a cylinder, the surface is a cone.

Homework Equations


Same as previous thread.

The Attempt at a Solution


I used cylindrical coordinates [itex](r, \phi , z)[/itex].
By intuition I know the angular momentum with respect to the z-axis (axis of symmetry of the cone) must be conserved, thus the Lagrangian must not contain [itex]\phi[/itex].
I reached that the Lagrangian [itex]L=\frac{m}{2}(\dot r ^2 + r^2 \dot \phi ^2 + \dot z ^2)[/itex].
Now, I believe I must express r in function of z. I notice that if [itex]\theta[/itex] is the angle worth half the angle of the vertex of the cone, then [itex]\tan \theta = constant = \frac{r}{z} \Rightarrow r=k_1z[/itex].
So now I have 2 cyclic coordinates and it means that the momentum [itex]P_r[/itex] is a constant. Since r is directly related to z, this also means that the z-component of the linear momentum is constant... Well I think so.

Calculating Lagrange's equations, I reach as equation of motion: [itex]\ddot z (K+1) - k_2 z \dot \phi ^2=0[/itex]. This can be written as [itex]K_3 \ddot z + K_2 z \dot \phi =0[/itex]. Where [itex]K_3 >1[/itex].
Am I right?
If so, how can I solve it?! What method would do the job?
 
Physics news on Phys.org
So now I have 2 cyclic coordinates and it means that the momentum [itex]P_r[/itex] is a constant. Since r is directly related to z, this also means that the z-component of the linear momentum is constant... Well I think so.

Your Lagrangian and your equation for r are both right, but this conclusion isn't. The Lagrangian's variables are not independent, so you need to somehow include the restraint r=k_1*z. There are two ways to do this: by directly substituting r=k_1*z into the Lagrangian, or by using Lagrange multipliers. If you use the first method, you get rid of r, so it becomes meaningless to talk about P_r. You also introduce an explicit dependence on z, so P_z isn't constant. If you use the second method, the Lagrangian would have an explicit dependence on both r and z, so neither P_r nor P_z are constant. Either way, P_r and P_z are not constants of motion.

Calculating Lagrange's equations, I reach as equation of motion: [itex]\ddot z (K+1) - k_2 z \dot \phi ^2=0[/itex]. This can be written as [itex]K_3 \ddot z + K_2 z \dot \phi =0[/itex]. Where [itex]K_3 >1[/itex].
Am I right?
If so, how can I solve it?! What method would do the job?

That's right, but you need another equation to deal with the phi dependence--the equation you get by differentiating the Lagrangian with respect to phi and its time derivative.
 
ideasrule said:
Your Lagrangian and your equation for r are both right, but this conclusion isn't. The Lagrangian's variables are not independent, so you need to somehow include the restraint r=k_1*z. There are two ways to do this: by directly substituting r=k_1*z into the Lagrangian, or by using Lagrange multipliers. If you use the first method, you get rid of r, so it becomes meaningless to talk about P_r. You also introduce an explicit dependence on z, so P_z isn't constant. If you use the second method, the Lagrangian would have an explicit dependence on both r and z, so neither P_r nor P_z are constant. Either way, P_r and P_z are not constants of motion.
Thank you very much. I hadn't been clear enough, I had replaced [itex]r[/itex] by [itex]k_1z[/itex] into my Lagrangian and I used it in order to calculate the Lagrange's equations.


That's right, but you need another equation to deal with the phi dependence--the equation you get by differentiating the Lagrangian with respect to phi and its time derivative.
Ah I see. This gives me [itex]2 \dot r \dot \phi + r \ddot \phi =0[/itex].
If I do the same for the [itex]r[/itex] coordinate, I get [itex]\ddot z =0[/itex]. I am not sure I can use this coordinate for the Lagrange's equations however. (In the link I gave in the first post of this thread, I've learned that if the coordinate is constant then I must NOT use it for the Lagrange's equations. Now, if a coordinate is dependent from another, I do not know, hence my doubt).
 
Ah I see. This gives me [itex]2 \dot r \dot \phi + r \ddot \phi =0[/itex].
If I do the same for the [itex]r[/itex] coordinate, I get [itex]\ddot z =0[/itex]. I am not sure I can use this coordinate for the Lagrange's equations however. (In the link I gave in the first post of this thread, I've learned that if the coordinate is constant then I must NOT use it for the Lagrange's equations. Now, if a coordinate is dependent from another, I do not know, hence my doubt).

Why is r still in the Lagrangian? I thought you replaced r by k_1*z.
 
ideasrule said:
Why is r still in the Lagrangian? I thought you replaced r by k_1*z.

My bad... sorry about this error. I did replace it but when I derived the Lagrangian with respect to phi I took my old Lagrangian by error.
I'll fix this as soon as possible. :redface: Thanks for pointing this out.
 
So I get [itex]2 \dot z \dot \phi +z \ddot \phi =0[/itex].
For z, I get [itex]\ddot z K_3+K_2 z \dot \phi ^2 =0[/itex].
 
fluidistic said:
So I get [itex]2 \dot z \dot \phi +z \ddot \phi =0[/itex].

I think you made an algebra error somewhere. The Lagrangian has no explicit dependence on phi, so the partial of L wrt the time derivative of phi should be constant.
 
ideasrule said:
I think you made an algebra error somewhere. The Lagrangian has no explicit dependence on phi, so the partial of L wrt the time derivative of phi should be constant.

Here is my Lagrangian: [itex]\frac{m}{2}(K \dot z ^2 + k_2 z^2 \dot \phi ^2 + \dot z ^2)[/itex].
[itex]\frac{\partial L}{\partial \dot \phi}= mk_2 z^2 \dot \phi[/itex].
Thus [itex]\frac{d}{dt} \left ( \frac{\partial L}{\partial \dot \phi} \right )=mk_2(2\dot z z \dot \phi+ z^2 \ddot \phi )[/itex]. Since [itex]\frac{\partial L}{\partial \phi}=0[/itex], I have that [itex]mk_2(2\dot z z \dot \phi+ z^2 \ddot \phi )=0[/itex] and since [itex]m[/itex], [itex]z[/itex] and [itex]k_2 \neq 0[/itex], I reach the expression I wrote in my previous post.
 
fluidistic said:
Here is my Lagrangian: [itex]\frac{m}{2}(K \dot z ^2 + k_2 z^2 \dot \phi ^2 + \dot z ^2)[/itex].
[itex]\frac{\partial L}{\partial \dot \phi}= mk_2 z^2 \dot \phi[/itex].
Thus [itex]\frac{d}{dt} \left ( \frac{\partial L}{\partial \dot \phi} \right )=mk_2(2\dot z z \dot \phi+ z^2 \ddot \phi )[/itex]. Since [itex]\frac{\partial L}{\partial \phi}=0[/itex], I have that [itex]mk_2(2\dot z z \dot \phi+ z^2 \ddot \phi )=0[/itex] and since [itex]m[/itex], [itex]z[/itex] and [itex]k_2 \neq 0[/itex], I reach the expression I wrote in my previous post.

My bad, that equation is actually correct. However, I don't think it's as useful as writing [itex]mk_2 z^2 \dot \phi=C[/itex], since the time derivative of [itex]\frac{\partial L}{\partial \dot \phi}[/itex] is 0. That way, you'd be able to express the time derivative of phi in terms of z, and substitute the result into the equation of motion for z.
 
  • #10
ideasrule said:
My bad, that equation is actually correct. However, I don't think it's as useful as writing [itex]mk_2 z^2 \dot \phi=C[/itex], since the time derivative of [itex]\frac{\partial L}{\partial \dot \phi}[/itex] is 0. That way, you'd be able to express the time derivative of phi in terms of z, and substitute the result into the equation of motion for z.
Oh this is bright.
I reach [itex]\ddot z z K_3+C_2=0[/itex]. Is this a Cauchy-Euler equation?
So I have only 1 equation of motion and I have to solve it in order to get the trajectory (since they ask for it)?
 

Similar threads

  • · Replies 4 ·
Replies
4
Views
3K
Replies
6
Views
2K
Replies
2
Views
1K
  • · Replies 3 ·
Replies
3
Views
4K
  • · Replies 3 ·
Replies
3
Views
2K
Replies
1
Views
3K
Replies
0
Views
1K
  • · Replies 9 ·
Replies
9
Views
4K
  • · Replies 11 ·
Replies
11
Views
3K
Replies
10
Views
2K