Eigen vals/vecs in MDOF Dynamics Problems

[JohnSimpson]

Say I have the following simple situation

||
||
||---/\/\/\---m1---/\/\/\----m2
---------------------------------

fixed walled and two masses connected by springs

When you write Newton's 2nd for each mass you can assume a sinusoid solution for the position as a function of time and plugging back into your equations you can solve for the two fundamental frequencies of the system. My question is this, If i write both equations in a matrix form

<br /> \frac{d^2 \, \vec{x}}{dt^2} = A \, \vec{x} <br />

with A being a matrix of combinations of spring constants divided by appropriate masses, do the eigenvalues of the matrix correspond to the squares of the natural angular frequencies? If so, why physically is this true? Additionally, what do the eigenvectors mean? What is the physical connection there?

-John

 

[kanato]

So, for this situation, if you write down the eigenvalue equation it looks like this:

<br /> \left(\begin{array}{cc}m_1^{-1} &amp; 0 \\ 0 &amp; m_2^{-1}\end{array}\right)<br /> \left(\begin{array}{cc}-k_1-k_2 &amp; k_2 \\ k_2 &amp; -k_2\end{array}\right)<br /> \left(\begin{array}{c}x_1\\x_2\end{array}\right) = <br /> \left(\begin{array}{c}\ddot x_1\\\ddot x_2\end{array}\right)<br />

(that's what I got, anyway. If it's not quite right, it doesn't really matter because the rest of this discussion doesn't focust on the form of the matrix.)

Now we assume there are harmonic solutions similar to a single harmonic oscillator, so that our trial solution would be \ddot x_1 = \omega^2 x_1 and \ddot x_2 = \omega^2 x_2. (This is done without loss of generality; \omega may be complex, so picking \omega^2 is really just handy for looking at the end result. But it's based on the fact that we've taken two derivatives wrt time, so we already know that omega has dimensions of frequency.) If you put that in:

<br /> \left(\begin{array}{cc}-(k_1+k_2)/m_1 &amp; k_2/m_1 \\ k_2/m_2 &amp; -k_2/m_2\end{array}\right)<br /> \left(\begin{array}{c}x_1\\x_2\end{array}\right) = <br /> \left(\begin{array}{cc}\omega^2 &amp; 0 \\ 0 &amp; \omega^2\end{array}\right)<br /> \left(\begin{array}{c} x_1\\ x_2\end{array}\right)<br />

which is exactly the eigenvalue problem. Solving it will give omega in terms of the k's and m's. Then we go back to our differential equation and look at our trial solution where we assumed \ddot x_i = \omega^2 x_i and solve that, knowing that we have a solution for omega. Of course that will be a harmonic solution; we don't need to know omega for that. But the omegas that we got out of the eigenvalue equation are the ones that go into our trial solution.

The eigenvectors have significance; a matrix made out of the eigenvectors can be used to diagonalize our k/m matrix. So this tells us what kind of x_1 and x_2 displacements give us the harmonic motion. If m1 = m2 and k1 = k2 then the eigenvectors, ignoring normalization, will look something like (1 1) and (1 -1). This tells us that one type of motion is where both masses move the same direction at the same time (1 1) and another type where they move opposite directions at the same time (1 -1). If the parameters are different, it gets more complicated, but my understanding is that it will be true that the ratio of the components of an eigenvector indicate the ratio of the amplitude of the motions of each mass. The eigenvectors could also be complex, which should indicate a phase difference in the motion of the two masses. This probably won't be the case for two masses, but more complicated systems could give this type of result I suppose.

 

[JohnSimpson]

So if i have this right

Assuming an SHM solution for each component, \omega_i^2 = \lambda_i

and the ratios of components in any eigenvector x_i give the ratios of he amplitudes of oscillation while in normal mode i