Solving the equations of motions of a particles constrained to a sphere.

  • Thread starter Thread starter VoxCaelum
  • Start date Start date
  • Tags Tags
    Particles Sphere
VoxCaelum
Messages
14
Reaction score
0

Homework Statement


Consider a particle moving on a two dimensional sphere of radius r, whose Lagrangian is given by:
L(q,q')=\frac{1}{2}m\sum(qi')2

(A) Transform the Lagrangian into spherical coordinates and write down the resulting Euler-Lagrange equations.

(B) Solve the Euler-Lagrange equations corresping to the Lagrangian and write down the solution in the original coordinates (qi). Express your answer to include the integration constants and denote them as the initial position and velocity.

Homework Equations


\frac{∂L}{∂q}=\frac{d}{dt}\frac{∂L}{∂q'}


The Attempt at a Solution


I think I have part (A) solved. The problem is part B. I really don't know where to start. But this is what I have for part (A)

mr2\phi''Sin2(θ)+\phi'θ'2Sin(θ)Cos(θ)=0
r2\phi'Sin(θ)Cos(θ)=mr2θ''

I would not know where to start solving these equations.
 
Physics news on Phys.org
First of all, you have some typos in your equations, since 2nd equation has 1 time derivative in left and 2 in right side. Also m should drop out everywhere.

Secondly, you forgot to calculate the last equation (wrt r), which will probably prove to be very useful in solving the set.
 
In the second equation I forgot a power of \phi'.
Also I did not forget to do one of the E.L. equations. Because we are constrained to the surface of a sphere I simply set r'=0 and treated r as a constant parameter instead of a dynamical variable. This reduces the Lagrangian (in spherical coordinates) to a Lagrangian involving two functions instead of three.
L=\frac{1}{2}m(r2θ'2+r2\phi'2Sin2(θ))

I also failed to notice that since the derivative of the lagrangian with respect to \phi is zero one of the E.L. equations will yield an integral of motion which should help out a bit.
 
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