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

In summary, the conversation discusses a problem involving a cone surface and the conservation of angular momentum. The Lagrangian and equation of motion are derived using cylindrical coordinates and the constraint r=k_1*z. However, there is an error in the derivation and the Lagrangian should not have an explicit dependence on phi. The partial of the Lagrangian with respect to the time derivative of phi should be constant.
  • #1
fluidistic
Gold Member
3,923
261

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
  • #2
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.
 
  • #3
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).
 
  • #4
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.
 
  • #5
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.
 
  • #6
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].
 
  • #7
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.
 
  • #8
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.
 
  • #9
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)?
 

1. What are the equations of motion for a particle moving over a cone?

The equations of motion for a particle moving over a cone can be derived using Newton's laws of motion and the geometry of the cone. The equations typically involve the variables of position, velocity, acceleration, and time.

2. How are the equations of motion affected by the shape of the cone?

The shape of the cone can affect the equations of motion in terms of the acceleration and velocity components in different directions. For example, a steeper cone may result in a greater acceleration in the downward direction compared to a shallower cone.

3. What are conserved quantities in this system?

Conserved quantities are physical quantities that remain constant throughout the motion of the particle over the cone. These quantities include energy, momentum, and angular momentum.

4. How do conserved quantities affect the equations of motion?

Conserved quantities play an important role in the equations of motion as they provide insights into the behavior of the system. For example, the conservation of energy can help determine the maximum height that the particle can reach on the cone.

5. Are there any real-life applications of this concept?

Yes, the equations of motion of a particle over a cone and the concept of conserved quantities have various real-life applications. These include studying the motion of objects rolling down a sloped surface, analyzing the flight of projectiles, and understanding the orbit of planets around the sun.

Similar threads

  • Advanced Physics Homework Help
Replies
4
Views
2K
  • Advanced Physics Homework Help
Replies
2
Views
687
Replies
6
Views
959
  • Advanced Physics Homework Help
Replies
0
Views
430
  • Advanced Physics Homework Help
Replies
11
Views
1K
  • Advanced Physics Homework Help
Replies
10
Views
1K
  • Advanced Physics Homework Help
Replies
9
Views
2K
  • Advanced Physics Homework Help
Replies
4
Views
914
  • Advanced Physics Homework Help
Replies
3
Views
1K
  • Advanced Physics Homework Help
Replies
16
Views
903
Back
Top