Ver el código de «Tres masas unidas por resortes (CMR)»

==Enunciado==
Se tienen tres masas, de valores <math>m_A=5\,\mathrm{kg}</math>, <math>m_B=5\,\mathrm{kg}</math> y <math>m_C=6\,\mathrm{kg}</math> se hallan unidas por resortes. Entre la A y la B se encuentra uno de constante <math>k_1=15\,\mathrm{N}/\mathrm{m}</math> y longitud natural 10cm, y entre la B y la C uno de constante <math>k_2=60\,\mathrm{N}/\mathrm{m}</math> y longitud natural 10cm. Inicialmente se encuentran las tres masas en reposo y los muelles en su longitud natural. Entonces se le comunica a la masa A una velocidad inicial <math>v_{A0}=+16\,\mathrm{cm}/\mathrm{s}</math>.

Determine la posición de cada masa como función del tiempo.
==Ecuaciones de movimiento==
En lo que sigue se usarán todas las magnitudes en las unidades fundamentales del SI.

Dado un sistema de referencia fijo, las posiciones de las tres partículas serán, respectivamente <math>X_A(t)</math>, <math>X_B(t)</math> y <math>X_C(t)</math>. Por comodidad, definimos las variables
<center><math>
x_A=X_A\qquad\qquad x_B = X_B-\ell_0\qquad\qquad x_C=X_C-2\ell_0</math></center>

siendo <math>\ell_0</math> la longitud natural de ambos resortes. De esta manera podemos describir el movimiento como las desviaciones respecto a posiciones de equlibrio y podemos tratar los resortes como si tuvieran longitud natural nula.

Las ecuaciones de movimiento las da la segunda ley de Newton. Para la primera masa

<center><math>m_A\ddot{x}_A=-k_1(x_A-x_B)\qquad\Rightarrow\qquad 5\ddot{x}_A=-15(x_A-x_B)\qquad\Rightarrow\qquad \ddot{x}_A=-3x_A+3x_B</math></center>

Para la segunda
<center><math>m_B\ddot{x}_B=-k_1(x_B-x_A)-k_2(x_B-x_C)\qquad\Rightarrow\qquad 5\ddot{x}_B=-15(x_B-x_A)-60(x_B-x_C)\qquad\Rightarrow\qquad \ddot{x}_B=3x_A-15x_B+12x_C</math></center>

y para la tercera

<center><math>m_C\ddot{x}_C=-k_2(x_C-x_B)\qquad\Rightarrow\qquad 6\ddot{x}_C=-60(x_C-x_B)\qquad\Rightarrow\qquad \ddot{x}_C=10x_B-10x_C</math></center>

Por tanto, el sistema de ecuaciones de movimiento es

<center><math>\begin{array}{rcl}
\ddot{x}_A&=&-3x_A+3x_B\\
\ddot{x}_B&=&3x_A-15x_B+12x_C\\
\ddot{x}_C&=&10x_B-10x_C
\end{array}</math></center>

con las condiciones iniciales

<center><math>\begin{array}{rcl}
x_{A0}&=&0\\
x_{B0}&=&0\\
x_{C0}&=&0
\end{array}\qquad\qquad \begin{array}{rcl}
\dot{x}_{A0}&=&+0.16\\
\dot{x}_{B0}&=&0\\
\dot{x}_{C0}&=&0
\end{array}</math></center>

==Modos normales==
El sistema de ecuaciones de movimiento puede escribirse en forma matricial como

<center><math>\frac{\mathrm{d}\ }{\mathrm{d}t^2}\begin{pmatrix}x_A\\x_B\\x_C\end{pmatrix}=\begin{pmatrix}-3&3&0\\3&-15&12\\0&10&-10\end{pmatrix}\cdot \begin{pmatrix}x_A\\x_B\\x_C\end{pmatrix}</math></center>

Existen diferentes formas de resolver este sistema de ecuaciones diferenciales de coeficientes constantes. 

Una de ellas consiste en buscar los llamados modos normales, que son oscilaciones del sistema con una frecuencia completa, de forma que la solución puede escribirse en la forma

<center><math>\begin{pmatrix}x_A\\x_B\\x_C\end{pmatrix} = \begin{pmatrix}\hat{x}_A\\\hat{x}_B\\\hat{x}_C\end{pmatrix}\mathrm{e}^{\mathrm{j}\omega t}</math></center>

Una solución de este tipo no cumplirá, en general, las condiciones iniciales. Sin embargo, puede demostrarse que la solución general es una combinación lineal de los diferentes modos normales. El problema se traslada entonces a encontrar los coeficientes adecuados para cada caso concreto.

Veamos primero cuáles son los modos normales en este caso. Sustituyendo la solución propuesta la ecuación diferencial se transforma en la ecuación algebraica

<center><math>-\omega^2\begin{pmatrix}\hat{x}_A\\\hat{x}_B\\\hat{x}_C\end{pmatrix}=\begin{pmatrix}-3&3&0\\3&-15&12\\0&10&-10\end{pmatrix}\cdot \begin{pmatrix}\hat{x}_A\\\hat{x}_B\\\hat{x}_C\end{pmatrix}</math></center>

o, agrupando términos,

<center><math>\begin{pmatrix}3-\omega^2&-3&0\\-3&15-\omega^2&-12\\0&-10&10-\omega^2\end{pmatrix}\cdot \begin{pmatrix}\hat{x}_A\\\hat{x}_B\\\hat{x}_C\end{pmatrix}=\begin{pmatrix}0\\0\\0\end{pmatrix}</math></center>

Este sistema es homogéneo. Para que haya soluciones distinta de l trivial (0,0,0) debe cumplirse que el determinante de la matriz de coeficientes se anule, en cuyo caso resulta un sistema indeterminado. Es decir, los valores de &omega; posibles vienen dados por la ,
<center><math>\left|
\begin{matrix}3-\omega^2&-3&0\\-3&15-\omega^2&-12\\0&-10&10-\omega^2\end{matrix}\right|=0</math></center>

Desarrollamos el determinante y resulta la ecuación

<center><math>-\omega^6 +28\omega^4 -96\omega^2=0\,</math></center>

Vemos que resulta una ecuación de sexto grado, por lo que hay 6 posibles modos normales. Una ecuación de sexto grado no tiene solución general, pero en este caso podemos obser var que en realidad se trata de una ecuación cúbica en <math>\omega^2</math>, que sí la tiene. Es más, en este caso podemos sacar factor común y escribirla como

<center><math>-\omega^2(\omega^4 -28\omega^2 +96)=0\,</math></center>

con lo que problema se reduce a resolver una ecuación se segundo grado en &omega;&sup2;. Esto nos da las soluciones para &omega;&sup2;

<center><math>\omega^2 = 0,\ 4,\ 24</math></center>

y para &omega;

<center><math>\omega = 0\ (\mathrm{doble}),\ \pm 2,\ \pm 2\sqrt{6}</math></center>

Dejaremos el caso de la solución <math>\omega=0</math> para el final y consideraremos antes el resto de casos.

===Caso &omega;&sup2; = 4 ===
Una vez que tenemos la frecuencia debemos buscar la forma de las amplitudes <center><math>\{\hat{x}_A,\hat{x}_B,\hat{x}_C\}</math></center> que no son arbitrarias, sino que deben cumplir la ecuación matricial anterior.

En el caso &omega;&sup2; = 4, la ecuación se convierte en 

<center><math>\begin{pmatrix}-1&-3&0\\-3&11&-12\\0&-10&6\end{pmatrix}\cdot \begin{pmatrix}\hat{x}_A\\\hat{x}_B\\\hat{x}_C\end{pmatrix}=\begin{pmatrix}0\\0\\0\end{pmatrix}</math></center>

lo que nos da las ecuaciones
<center><math>-\hat{x}_A-3\hat{x}_B=0\qquad -3\hat{x}_A+11\hat{x}_B-12\hat{x}_C=0\qquad -10\hat{x}_B+6\hat{x}_C=0
</math></center>

Estas ecuaciones no son las tres independientes, al ser nulo el determinante. Podemos determinar dos variables en función de la tercera. Por ejemplo

<center><math>\hat{x}_A=-3\hat{x}_B\qquad\qquad \hat{x}_C=\frac{5}{3}\hat{x}_B</math></center>

Dada una solución, cualquier múltiplo de ella es también solución. Podemos elegir arbitrariamente <math>\hat{x}_B=3</math> y tener la solución

<center><math>\hat{x}_A=-9\qquad\qquad \hat{x}_B=3\qquad\qquad \hat{x}_C=5</math></center>

Por tanto, los dos primeros modos normales son

<center><math>\begin{pmatrix}-9\\3\\5\end{pmatrix}\mathrm{e}^{2\mathrm{j}t}\qquad\qquad \begin{pmatrix}-9\\3\\5\end{pmatrix}\mathrm{e}^{-2\mathrm{j}t}</math></center>
===Caso &omega;&sup2; = 24 ===
En el segundo caso, operamos de manera análoga. Sustituyendo &omega;&sup2; = 24 resulta la ecuación

<center><math>\begin{pmatrix}-21&-3&0\\-3&-6&-12\\0&-10&-14\end{pmatrix}\cdot \begin{pmatrix}\hat{x}_A\\\hat{x}_B\\\hat{x}_C\end{pmatrix}=\begin{pmatrix}0\\0\\0\end{pmatrix}</math></center>

cuyas soluciones son los múltiplos de

<center><math>\hat{x}_A=1\qquad\qquad\hat{x}_B=-7\qquad\qquad\hat{x}_C=5</math></center>

Por tnto, los siguientes dos modos normales son de la forma

<center><math>\begin{pmatrix}1\\-7\\5\end{pmatrix}\mathrm{e}^{2\sqrt{6}\mathrm{j}t}\qquad\qquad \begin{pmatrix}1\\-7\\5\end{pmatrix}\mathrm{e}^{-2\sqrt{6}\mathrm{j}t}</math></center>

En realidad, tanto en este caso como en el anterior, las dos soluciones corresponden al mismo modo de oscilación, entendido como forma en la que se mueven las partículas. 

===Caso &omega;&sup2; = 0 ===
El caso &omega;&sup2; = 0 es un poco especial ya que no se trata de una solución oscilante. Si la frecuencia es nula, lo que se cumple en realidad es que 

<center><math>\frac{\mathrm{d}\ }{\mathrm{d}t^2}\begin{pmatrix}x_A\\x_B\\x_C\end{pmatrix}=\begin{pmatrix}0\\0\\0\end{pmatrix}</math></center>

y si la segunda derivada es nula, eso quiere decir que cada variable esde la forma <math>x=A + B t</math>, es decir corresponde a un movimiento uniforme. No obstante, las tres variables no son independientes, ya que deben seguir cumpliendo el sistema de ecuaciones

<center><math>\begin{pmatrix}3&-3&0\\-3&15&-12\\0&-10&10\end{pmatrix}\cdot \begin{pmatrix}\hat{x}_A\\\hat{x}_B\\\hat{x}_C\end{pmatrix}=\begin{pmatrix}0\\0\\0\end{pmatrix}</math></center>

cuya solución es simplemente

<center><math>\hat{x}_A=\hat{x}_B=\hat{x}_C=1</math></center>

Por tanto, las dos últimas soluciones independientes son de la forma

<center><math>\begin{pmatrix}1\\1\\1\end{pmatrix}\qquad\qquad \begin{pmatrix}1\\1\\1\end{pmatrix}t</math></center>

El significado físico de estas dos soluciones es sencillo. 

La primera corresponde a que estando las tres masas a las distancias de equiibrio las desplazamos una cierta cantidad fija e igual para las tres. El sistema sigue estando en equilibrio.

La segunda corresponde a un movimiento uniforme en que las tres masas se desplazan con la misma velocidad, manteniendo las distancias respectivas.
==Solución particular==

Una vez que tenemos las 6 soluciones independientes, la solución general es de la forma

<center><math>\begin{pmatrix}x_A\\x_B\\x_C\end{pmatrix}= c_1 \begin{pmatrix}-9\\3\\5\end{pmatrix}\mathrm{e}^{2\mathrm{j}t}+c_2 \begin{pmatrix}-9\\3\\5\end{pmatrix}\mathrm{e}^{-2\mathrm{j}t}+c_3\begin{pmatrix}1\\-7\\5\end{pmatrix}\mathrm{e}^{2\sqrt{6}\mathrm{j}t}+c_4\begin{pmatrix}1\\-7\\5\end{pmatrix}\mathrm{e}^{-2\sqrt{6}\mathrm{j}t}+c_5\begin{pmatrix}1\\1\\1\end{pmatrix}+c_6\begin{pmatrix}1\\1\\1\end{pmatrix}t</math></center>

Siendo las velocidades

<center><math>\begin{pmatrix}\dot{x}_A\\\dot{x}_B\\\dot{x}_C\end{pmatrix}=2\mathrm{j} c_1 \begin{pmatrix}-9\\3\\5\end{pmatrix}\mathrm{e}^{2\mathrm{j}t}-2\mathrm{j}c_2 \begin{pmatrix}-9\\3\\5\end{pmatrix}\mathrm{e}^{-2\mathrm{j}t}+2\sqrt{6}\mathrm{j}c_3\begin{pmatrix}1\\-7\\5\end{pmatrix}\mathrm{e}^{2\sqrt{6}\mathrm{j}t}-2\sqrt{6}\mathrm{j}c_4\begin{pmatrix}1\\-7\\5\end{pmatrix}\mathrm{e}^{-2\sqrt{6}\mathrm{j}t}+c_6\begin{pmatrix}1\\1\\1\end{pmatrix}</math></center>

Para hallar los coeficientes imponemos las condiciones iniciales. Hacemos <math>t =0</math> y queda

<center><math>\begin{pmatrix}0\\0\\0\end{pmatrix}= c_1 \begin{pmatrix}-9\\3\\5\end{pmatrix}+c_2 \begin{pmatrix}-9\\3\\5\end{pmatrix}+c_3\begin{pmatrix}1\\-7\\5\end{pmatrix}+c_4\begin{pmatrix}1\\-7\\5\end{pmatrix}+c_5\begin{pmatrix}1\\1\\1\end{pmatrix}</math></center>

y

<center><math>\begin{pmatrix}0.16\\0\\0\end{pmatrix}=2\mathrm{j} c_1 \begin{pmatrix}-9\\3\\5\end{pmatrix}-2\mathrm{j}c_2 \begin{pmatrix}-9\\3\\5\end{pmatrix}+2\sqrt{6}\mathrm{j}c_3\begin{pmatrix}1\\-7\\5\end{pmatrix}-2\sqrt{6}\mathrm{j}c_4\begin{pmatrix}1\\-7\\5\end{pmatrix}+c_6\begin{pmatrix}1\\1\\1\end{pmatrix}</math></center>

Esto nos da el sistema de ecuaciones

<center><math>\begin{array}{rcl}
-9(c_1+c_2)+(c_3+c_4)+c_5&=&0\\
3(c_1+c_2)-7(c_3+c_4)+c_5&=&0\\
5(c_1+c_2)+5(c_3+c_4)+c_5&=&0\\
-18\mathrm{j}(c_1-c_2)+2\sqrt{6}\mathrm{j}(c_3-c_4)+c_6&=&0.16\\
6\mathrm{j}(c_1-c_2)-14\sqrt{6}\mathrm{j}(c_3-c_4)+c_6&=&0\\
10\mathrm{j}(c_1-c_2)+10\sqrt{6}\mathrm{j}(c_3-c_4)+c_6&=&0
\end{array}</math></center>

Con solución

<center><math>c_1=0.003\mathrm{j}\qquad c_2=-0.003\mathrm{j}\qquad c_3=-\frac{0.0005\mathrm{j}}{\sqrt{6}}\qquad c_4=+\frac{0.0005\mathrm{j}}{\sqrt{6}}\qquad\qquad c_5=0\qquad\qquad c6=0.05</math></center>

Llevando esto a la expresión general queda

<center><math>\begin{pmatrix}x_A\\x_B\\x_C\end{pmatrix}= -0.003\mathrm{j}\begin{pmatrix}-9\\3\\5\end{pmatrix}\mathrm{e}^{2\mathrm{j}t}+0.003\mathrm{j}\begin{pmatrix}-9\\3\\5\end{pmatrix}\mathrm{e}^{-2\mathrm{j}t}-\frac{0.0005\mathrm{j}}{\sqrt{6}}\begin{pmatrix}1\\-7\\5\end{pmatrix}\mathrm{e}^{2\sqrt{6}\mathrm{j}t}+\frac{0.0005\mathrm{j}}{\sqrt{6}}\begin{pmatrix}1\\-7\\5\end{pmatrix}\mathrm{e}^{-2\sqrt{6}\mathrm{j}t}+0,05\begin{pmatrix}1\\1\\1\end{pmatrix}t</math></center>

Aplicando aquí la fórmula de Euler
	
<center><math>\mathrm{e}^{\mathrm{j}\alpha}=\cos(\alpha)+\mathrm{j}\,\mathrm{sen}(\alpha)</math></center>

resulta

<center><math>\begin{pmatrix}x_A\\x_B\\x_C\end{pmatrix}= -0.006\mathrm{j}\begin{pmatrix}-9\\3\\5\end{pmatrix}\mathrm{sen}(2t)+\frac{0.001}{\sqrt{6}}\begin{pmatrix}1\\-7\\5\end{pmatrix}\mathrm{sen}(2\sqrt{6}t)+0,05\begin{pmatrix}1\\1\\1\end{pmatrix}t</math></center>

o separando por componentes

<center><math>\begin{array}{lcr}
x_A&=&0.054\,\mathrm{sen}(2t)+\dfrac{0.001}{\sqrt{6}}\mathrm{sen}(2\sqrt{6}t)+0.05t\\&&\\
x_B&=&-0.018\,\mathrm{sen}(2t)+\dfrac{0.007}{\sqrt{6}}\mathrm{sen}(2\sqrt{6}t)+0.05t\\&&\\
x_A&=&-0.030\,\mathrm{sen}(2t)+\dfrac{0.005}{\sqrt{6}}\mathrm{sen}(2\sqrt{6}t)+0.05t
\end{array}</math></center>

==Solución en el sistema centro de masas==
Una forma alternativa de abordar este problema consiste en identificar en primer lugar el movimiento del CM y a continuación analizar como se mueven las masas respecto al CM.

Por ser todas las fuerzas internas, el movimiento del CM es rectilíneo y uniforme, siendo su posición y velocidad inicial

<center><math>x_{G0}=\frac{m_Ax_{A0}+m_Bx_{B0}+m_Cx_{C0}}{m_A+m_B+m_C}=0\qquad\qquad v_{G0}=\frac{m_Av_{A0}+m_Bv_{B0}+m_Cv_{C0}}{m_A+m_B+m_C}=0.05 </math></center>   

(posteriormente a <math>x_G</math> habría que sumarle <math>11\ell_0/16</math>, que es su posición teniendo en cuenta las longitudes de los muelles)
Si ahora definimos coordenadas de cada masa respecto al CM

<center><math>y_A=x_A-x_G\qquad\qquad y_B=x_B-x_G\qquad\qquad y_C=x_C-x_G</math></center>

Estas tres coordenadas cumplen el vínculo

<center><math>m_Ay_A+m_By_B+m_Cy_C=0\qquad\qquad 5y_A+5y_B+6y_C=0</math></center>

Las ecuaciones diferenciales para las y's son las mismas para las x's

<center><math>\begin{array}{rcl}
\ddot{y}_A&=&-3y_A+3y_B\\
\ddot{y}_B&=&3y_A-15y_B+12y_C\\
\ddot{y}_C&=&10y_B-10y_C
\end{array}</math></center>

pero sustituyendo la ecuación de vínculo se puede reducir a solo dos

<center><math>\begin{array}{rcl}
\ddot{y}_A&=&-3y_A+3y_B\\
\ddot{y}_B&=&3y_A-15y_B+2(-5y_A-5y_B)=-7y_A-25y_B
\end{array}</math></center>

Con esto reduce el análisis a una matriz 2&times;2 y una ecuación de segundo grado en &omega;&sup2;. El resultado final es, lógicamente, el mismo que antes.