1. Not finding help here? Sign up for a free 30min 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!

Coupled harmonic oscillators

  1. Dec 10, 2007 #1
    1. The problem statement, all variables and given/known data

    A particle of mass m1 is attached to a wall by a spring of constant k. A second particle of mass m2 is attached to a different wall by another sping of constant k. The two masses are attached to eachother by a third spring of constant k. Let [tex]x_1[/tex] and [tex]x_2[/tex] be the displacement of the respective particles from their equilibrium positions. Find the normal coordinates for the system.

    This problem either with spings or pendula appears in many classical mechanics textbooks (I am using the book by Thornton & Marion). The wrinkle is that the two particles have different masses in this case.


    2. Relevant equations

    for the first mass:
    [tex]m_{1} \ddot{x_1} = -k x_1 - k \left(x_1 - x_2 \right)[/tex]
    for the second mass:
    [tex]m_{2} \ddot{x_2} = -k x_2 - k \left(x_2 - x_1 \right)[/tex]

    the eigenfrequencies arise from taking [tex]\texttt{det} \left( A_{j k} - \omega^2 m_{j k} \right) = 0 [/tex]
    where [tex]m_{j k} = ( (m_1, 0), (0, m_2) ) [/tex] and [tex] A_{j k} = ( (2 k, -k), (-k, 2 k) ) [/tex] are both 2 x 2 matrices.

    These are used to find conditions on the eigenvectors [tex]a_{j r}[/tex] which satisfy
    [tex]\left( A_{j k} - \omega_r^2 m_{j k} \right) a_{j r} = 0[/tex]
    where [tex]\omega_r^2[/tex] is the rth eigenfrequency.

    The eigenvectors transform the original coordinates [tex]x_j[/tex] into the normal coordinates [tex]\eta_j[/tex] via
    [tex]x_j = \sum_r a_{j r} \eta_r[/tex].





    3. The attempt at a solution

    The eigenfrequencies are:
    [tex]\omega^2_{\pm} = \frac{k}{m_1 m_2} \left( m_1 + m_2 \pm \sqrt{m_1^2 - m_1 m_2 + m_2^2} \right) [/tex]

    The (unnormalized) eigenvectors are
    [tex] a_{j \pm} = \left(1, \frac{1}{m_2} \left( - m_1 + m_2 \pm \sqrt{m_1^2 - m_1 m_2 + m_2^2} \right) \right) [/tex]

    Then the new coordinates are found by
    [tex]x_1 = \eta_1 + \eta_2 [/tex]
    [tex]x_2 = \left( m_1 + m_2 + \sqrt{m_1^2 - m_1 m_2 + m_2^2} \right) \eta_2 + \left( m_1 + m_2 - \sqrt{m_1^2 - m_1 m_2 + m_2^2} \right) \eta_2 [/tex]

    Note that in the limit [tex]m_1 = m_2[/tex] we obtain the textbook result
    [tex]x_1 = \eta_1 + \eta_2 [/tex]
    [tex]x_2 = \eta_2 - \eta_1[/tex]

    So it seems like I have solved the problem - but I haven't. If I check the solution by substituting the new coordinates into the equations of motion I find that the equations have not been decoupled, so these cannot be the normal coordinates.

    I wish anyone who attempts this problem good luck. I have attached a Mathematica program which shows this result.

    Chris
     

    Attached Files:

    Last edited: Dec 11, 2007
  2. jcsd
  3. Dec 10, 2007 #2
    Well, this is a pretty long problem, so let's start towards the beginning. When you take the [tex]\texttt{det} \left( A_{j k} - \omega^2 m_{j k} \right) = 0 [/tex] what is the quadratic equation you get for omega squared? I have done the problem up to this point but want to make sure we're on the same page before I proceed with the ridiculous amount of algebra/simplifying. Yes, sorry for the laziness, haha.
     
  4. Dec 11, 2007 #3
    The eigenvalue equation is
    [tex]\left( 2 k - \omega^2 m_1 \right) \left( 2 k - \omega^2 m_2 \right) - k^2 = 0 [/tex]
    which has two roots in the variable [tex]\omega^2[/tex]:
    [tex]\omega^2_{\pm} = \frac{k}{m_1 m_2} \left( m_1 + m_2 \pm \sqrt{m_1^2 - m_1 m_2 + m_2^2} \right) [/tex]

    You are right, there is a rediculous amount of algebra which is why I opted for a computer algebra package (see the file attached to my original post). I believe my program is correct and it does give the correct normal coordinates for [tex]m_1 = m_2 [/tex].
     
  5. Dec 11, 2007 #4
    Ok, I got the same characteristic frequencies that you did.

    Given the formulas you've posted, your normalized coordinates cannot be right. When I plug in m1=m2 into x2, I get something like 4*m*n2 which cannot be right. I think one of your etas needs to switch an index.

    When I computed the eigenfrequencies, I multiplied that matrix out and just set the first component equal to 1, and then solved for the 2nd component, is this what you did?
     
  6. Dec 11, 2007 #5
    Yes, this sounds similar to what I did. I made an error in my original post regarding the eigenvectors :blushing:. The correct eigenvectors are:
    [tex] a_{j (1)} = \left(1, \frac{1}{m_2} \left( - m_1 + m_2 + \sqrt{m_1^2 - m_1 m_2 + m_2^2} \right) \right) [/tex]
    [tex] a_{j (2)} = \left(1, \frac{1}{m_2} \left( - m_1 + m_2 - \sqrt{m_1^2 - m_1 m_2 + m_2^2} \right) \right) [/tex]

    Then
    [tex] x_1 = a_{1 (1)} \eta_1 + a_{1 (2)} \eta_2 = \eta_1 + \eta_2 [/tex]

    [tex] x_2 = a_{2 (1)} \eta_1 + a_{2 (2)} \eta_2[/tex]
    [tex] = \frac{1}{m_2} \left( - m_1 + m_2 + \sqrt{m_1^2 - m_1 m_2 + m_2^2} \right) \eta_1 + \frac{1}{m_2} \left( - m_1 + m_2 - \sqrt{m_1^2 - m_1 m_2 + m_2^2} \right) \eta_2 [/tex]

    When [tex]m_1 = m_2 [/tex] this gives the right result. Sorry for the error. Thanks for picking it up.

    Chris
     
  7. Dec 12, 2007 #6
    So, wait. Did this allow you to solve the problem or was that a typo?
     
  8. Dec 12, 2007 #7
    It was a typo. If you sub the "normal coordinates" just found into the equations of motion you will find they do not separate as they should. Or I don't see how they separate.
     
  9. Dec 13, 2007 #8
    I'm getting different eigenvectors than you because I don't know how you went about solving for yours, so I just want to check if yours are ok. Here's my method: multiplying out the your first equation and letting the components of my eigenvectors be 'a' and 'b' I get:

    [tex]a(-m_1\omega^2+2k)-bk=0[/tex]
    [tex]-ak+b(-m_2\omega^2+2k)=0[/tex]

    Letting a=1, and solving for b from the first equation

    [tex]b=2-\frac{m_1\omega^2}{k}[/tex]

    Plugging in for [tex]\omega^2[/tex], I get eigenvectors:

    [tex](1,2-\frac{1}{m_1}(m_1+m_2\pm\sqrt{m_1^2-m_1m_2+m_2^2}))[/tex]
    which also checks with the limit
    [tex]m_1=m_2[/tex]

    If your method matches mine for the eigenvectors, I suggest tinkering with the final equation when you sub the normal coordinates back into the original differential equation and just keep playing with it until you get the right answer or proof that you messed up somewhere.
     
  10. Dec 13, 2007 #9
    It seems like we arrive at the same eigenvectors which is reassuring. Thanks for all your hard work on this problem.

    I am becoming skeptical that normal modes exist for cases when the masses are different or the sping constant connecting m1 to the wall is different than the spring constant connecting m2 to the wall. In these cases the two masses will naturally oscillate at different frequencies. Consequently it is impossible for their motion to be completely in phase or completely out of phase as is the case in the textbook treatment of normal modes (i.e. m1 = m2).
     
    Last edited: Dec 13, 2007
Know someone interested in this topic? Share this thread via Reddit, Google+, Twitter, or Facebook

Have something to add?



Similar Discussions: Coupled harmonic oscillators
Loading...