1. Limited time only! Sign up for a free 30min personal tutor trial with Chegg Tutors
    Dismiss Notice
Dismiss Notice
Join Physics Forums Today!
The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

Homework Help: Help with example from goldstein (lagrangian)

  1. Apr 14, 2010 #1
    1. The problem statement, all variables and given/known data
    From pages 124-125 in edition 3.

    This is about the restricted three body problem (m3 << m1,m2)

    [PLAIN]http://img718.imageshack.us/img718/7012/3bdy.jpg [Broken]


    2. Relevant equations
    L = T-V

    Euler-Lagrange equations

    3. The attempt at a solution

    I'm interested in m3, the negligible mass object.
    Latex isn't showing up in the preview so note that by xdot I mean d/dt(x), abs = absolute value, and r1,r2,r3 are vectors.

    L = 1/2 m3 (xdot^2+ydot^2) - V

    my first question is: is the potential simply this?
    [PLAIN]http://img408.imageshack.us/img408/9834/potential.jpg [Broken]

    and I'm calling this V(r,theta,t)

    changing to polar coordinates (x=rcos(theta), y=rsin(theta)) and simplifying gives

    [PLAIN]http://img11.imageshack.us/img11/6784/14361215.jpg [Broken]

    This is where the example starts, except the potential is not stated, it is just given as V(r,theta,t)
    the vectors in my potential depend on r, theta and t so I guess that is right?

    Next in the example, the coordinate system is made to rotate at the same frequency as m3 so that m1 and m2 appear stationary.
    theta' = theta + wt (I'm using w for omega)
    theta = theta' - wt
    thetadot = theta'dot - w
    so I replace thetadot above in the equation for L and I eventually get
    [PLAIN]http://img198.imageshack.us/img198/2047/lchanged1.jpg [Broken]

    now my main confusion comes. In the example, as well as making the coordinate system rotate, they say switch to cylindrical coordinates and state that this becomes

    [PLAIN]http://img685.imageshack.us/img685/421/lchanged2.jpg [Broken]

    questions:
    The t dependency in the potential changes into a z dependency?
    z seems a strange choice as it is a two-dimensional problem?

    The only reason I can see for this is to get the lagrangian in the form L = (rho,theta,z,rhodot,thetadot,zdot)

    The book does not say what the transformations were except for theta. Just that it was a cylindrical transform using rho, theta and z

    comparing the two equations (the one I derived and the one given) the transformation must have been

    rho = r (since x and y are already in polar coordinates x=rcos(theta)=rhocos(theta))
    theta = theta' - wt
    z = t
    and then they added the zdot^2 to the kinetic energy term which is fine since it is 0 anyway.

    but can you simply make a change like this? z = t doesn't feel right.
    So I'm stuck up to this point.

    The book then says to find the euler-lagrange equations and look for solutions where rhodot=zdot=thetadot=0 and these will be the five lagrange points.

    So the first thing I need is a complete lagrangian with a defined and transformed potential in order to get the derivatives for the euler-lagrange equations.

    Hopefully someone can help me out, thanks
     
    Last edited by a moderator: May 4, 2017
  2. jcsd
  3. Apr 14, 2010 #2

    gabbagabbahey

    User Avatar
    Homework Helper
    Gold Member

    Shouldn't there be another term in your potential...namely the interaction between [itex]m_1[/itex] and [itex]m_2[/itex]?

    Also, are you treating [itex]x_1[/itex], [itex]x_2[/itex], [itex]r_1[/itex], [itex]r_2[/itex] and [itex]r_3[/itex] as distances or vectors? If you are treating them as distances, then [itex]x_1\neq r_2-r_3[/itex] and [itex]x_2\neq r_1-r_3[/itex]. If you are treating them as vectors, then your potential (which is always a scalar!) involves only the magnitudes of [itex]\textbf{x}_1[/itex], [itex]\textbf{x}_2[/itex] and [itex]\textbf{r}_2-\textbf{r}_1[/itex] (from your missing 3rd term).
     
    Last edited by a moderator: May 4, 2017
  4. Apr 14, 2010 #3
    Thanks for that, those are some pretty big mistakes. They are vectors as that's how they are in the book.

    About the 3rd term in the potential, I agree it would be needed if it was the potential was for the system but as far as I understand I'm only interested in the potentials that act upon m3? Just as I have excluded the kinetic energies of m2 and m1, I exclude the potential that acts between them? Let me know if I'm wrong.

    I'll add the third term for now though as I can always change it later

    [PLAIN]http://img693.imageshack.us/img693/9287/newvandl.jpg [Broken]

    So I'm assuming this is right apart from whether the 3rd term in the potential should be included or not.

    Then the following transformation is applied to make to coordinate system rotate at the same frequency.

    [PLAIN]http://img534.imageshack.us/img534/1777/rotatep.jpg [Broken]

    Then thetadot in the equation for L is changed and expanded giving
    [PLAIN]http://img198.imageshack.us/img198/2047/lchanged1.jpg [Broken]

    Then the book says switch to cylindrical coordinates I'll quote here:
    Which is exactly the same, just changing r to rho and t to z?
    In the standard cylindrical coordinates, z = z. So why does the t vanish and become a z in the book? z is not just another variable for the time either as they add in the zdot^2 term in the kinetic energy implying they mean the kinetic energy in the z direction. So how did my time dependent potential turn into a z dependent potential? Especially since it is a two-dimensional problem
     
    Last edited by a moderator: May 4, 2017
  5. Apr 15, 2010 #4
    I tried a different approach not using vectors, is this valid?

    [PLAIN]http://img59.imageshack.us/img59/9499/45953816.jpg [Broken]

    [PLAIN]http://img411.imageshack.us/img411/7677/42701744.jpg [Broken]

    [PLAIN]http://img265.imageshack.us/img265/7100/33037853.jpg [Broken]

    and I then get the Euler-Lagrange equations for r and theta. When differentiating w.r.t r, am I correct in thinking that r1 and r2 and constants and the only thing I have to differentiate is r?
     
    Last edited by a moderator: May 4, 2017
Share this great discussion with others via Reddit, Google+, Twitter, or Facebook