Ver el código de «Dos masas unidas sobre una cuchilla»

==Enunciado==
Dos masas iguales m están unidas por una varilla rígida ideal de longitud b. La varilla reposa sobre un plano horizontal. una de masas (la “2) puede deslizar sin rozamiento sobre el plano, pero la “1” está montada sobre una cuchilla que la obliga a desplazarse solo en la dirección paralela a la propia varilla 
# ¿Qué vínculos ligan las posiciones y velocidades de las partículas? 
# ¿Hacia dónde van dirigidas las fuerzas de reacción vincular?
# Escriba el sistema de ecuaciones de movimiento y de vínculos para este sistema empleando como variables las coordenadas cartesianas de la masa 1 y el ángulo que la varilla forma con el eje OX. Sugerencia: emplee tanto el sistema de referencia fijo “1” como el “2” ligado a la varilla.
# Introduciendo las variables adecuadas, reduzca este problema a un sistema de ecuaciones de primer orden.
<center>[[Archivo:Dos-masas-cuchilla.png|400px]]</center>

==Vínculos==
Si identificamos la posición de cada partícula en el plano por sus coordenadas cartesianas, los vínculos sobre el sistema son tales que:

* La distancia entre las partículas es constante

<center><math>|\overrightarrow{AB}|=b\qquad\Rightarrow\qquad(x_B-x_A)^2+(y_B-y_A)^2 = b^2\,</math></center>

* La velocidad de la partícula 1 es paralela al vector de posición relativa

<center><math>\vec{v}_A\times\overrightarrow{AB}=\vec{0}\qquad\Rightarrow\qquad (x_B-x_A)\dot{y}_A-(y_B-y_A)\dot{x}_A=0</math></center>

:Este vínculo no puede integrarse para dar un vínculo geométrico. No se puede integrar porque aparecen las derivadas de las coordenadas de la partícula 1 multiplicadas por las coordenadas de la partícula 2, que no dependen de la 1. Por ello no se puede escribir como la derivada de una sola función.

Si usamos como coordenadas las indicadas en el enunciado

<center><math>x_A = x\qquad\qquad y_A = y\qquad\qquad x_B=x+bC\qquad\qquad y_B=y+bS</math></center>

(con C y S representando al coseno y el seno de &theta;), el primer vínculo se vuelve trivial

<center><math>b=\mathrm{cte.}\,</math></center>

y el segundo se escribe

<center><math>-S\dot{x}+C\dot{y}=0</math></center>
==Fuerzas de reacción==
Asociadas a cada vínculo se producen fuerzas de reacción vincular.

* La partícula 1 se ve sometida a una fuerza por cada vínculo. Para el primero, la fuerza va en la dirección de la varilla. 

<center><math>\vec{F}_{T1}=F_T(C\vec{\imath}+S\vec{\jmath})</math></center>

:Para el segundo, es perpendicular a la dirección de movimiento posible y por tanto, perpendicular a la varilla 

<center><math>\vec{F}_{n1}=F_n(-S\vec{\imath}+C\vec{\jmath})</math></center>

* Sobre la segunda partícula solo actúa la tensión de la varilla, en sentido opuesto a la que actúa sobre la 1.

<center><math>\vec{F}_{T2}=-F_T(C\vec{\imath}+S\vec{\jmath})</math></center>

==Ecuaciones de movimiento==
Si separamos por componentes, las ecuaciones de movimiento son

<center><math>\begin{array}{rcl}
m\ddot{x}_1 & = & F_TC-F_nS\\
m\ddot{y}_1 & = & F_TS+F_nC\\
m\ddot{x}_2 & = & -F_TC\\
m\ddot{y}_2 & = & -F_TS
\end{array}</math></center>

En términos de las coordenadas sugeridas, derivamos dos veces <math>x_2</math> e <math>y_2</math>

<center><math>\begin{array}{rcl}
m\ddot{x} & = & F_TC-F_nS\\
m\ddot{y} & = & F_TS+F_nC\\
m(\ddot{x}-b\ddot{\theta}S-b\dot{\theta}^2C) & = & -F_TC\\
m(\ddot{y}+b\ddot{\theta}C-b\dot{\theta}^2S) & = & -F_TS
\end{array}</math></center>

Podemos reducir este sistema restando la tercera de la primera y la segunda de la cuarta para obtener

<center><math>\begin{array}{rcl}
mb(\ddot{\theta}S+\dot{\theta}^2C) & = & 2F_TC-F_nS\\
mb(-\ddot{\theta}C+\dot{\theta}^2S) & = & 2F_TS+F_nC
\end{array}</math></center>

Si ahora multiplicamos la primera por el coseno, la segunda por el seno y sumamos

<center><math>mb\dot{\theta}^2 = 2F_T</math></center>

y si multiplicamos la primera por el seno, la segunda por el coseno y restamos

<center><math>mb\ddot{\theta}=-F_n</math></center>

Esto nos permite eliminar las fuerzas de reacción vincular

<center><math>\begin{array}{rcl}
m\ddot{x} & = & mb\left(\dfrac{\dot{\theta}^2C}{2}+\ddot{\theta}S\right)\\
m\ddot{y} & = & mb\left(\dfrac{\dot{\theta}^2S}{2}-\ddot{\theta}C\right)
\end{array}</math></center>

Tenemos aqui dos ecuaciones con tres variables. Para completar el sistema hay que añadir la ecuación del vínculo no holónomo

<center><math>-S\dot{x}+C\dot{y}=0</math></center>

Hay que remarcar que aunque tengamos tres variables, el sistema tiene solo dos grados de libertad.

==Ecuaciones de movimiento. Forma alternativa==
Estas ecuaciones, o unas equivalentes, pueden obtenerse también empleando dos sistemas de referencia. Definimos un sistema de referencia ligado con origen en A y cuyo eje <math>OX_2</math> va en la dirección de la varilla. 
===Posiciones===
Empleando los dos sistemas, la posición de A es

<center><math>\overrightarrow{OA}=x\vec{\imath}_1+y\vec{\jmath}_1</math></center>
y la de B

<center><math>\overrightarrow{AB}=b\vec{\imath}_2 \qquad\Rightarrow\qquad \overrightarrow{OB}=x\vec{\imath}_1+y\vec{\jmath}_1+b\vec{\imath}_2</math></center>
===Velocidades===
La velocidad de A es

<center><math>\vec{v}^A_{21}=\dot{x}\vec{\imath}_1+\dot{y}\vec{\jmath}_1</math></center>

Pero esta velocidad solo puede apuntar en la dirección de la varilla, por lo que debe ser de la forma

<center><math>\vec{v}^A_{21}=v\vec{\imath}_2</math></center>

cumpliéndose las relaciones

<center><math>\dot{x}=v\cos(\theta)\qquad\qquad \dot{y}=v\,\mathrm{sen}(\theta)</math></center>

La velocidad de B cumple

<center><math>\vec{v}^B_{21}=\vec{v}^A_{21}+\vec{\omega}_{21}\times\overrightarrow{AB}=v\vec{\imath}+\omega b\vec{\jmath}_2</math></center>

siendo <math>\omega=\dot{\theta}</math> la velocidad angular con que gira la varilla. Se cumple

===Aceleraciones===
Derivando la expresión de la velocidad de A obtenemos

<center><math>\vec{a}^A_{21}=\dot{v}\vec{\imath}_2+v\frac{\mathrm{d}\vec{\imath}_2}{\mathrm{d}t}=\dot{v}\vec{\imath}_2+v\omega\vec{\jmath}_2</math></center>

ya que

<center><math>\frac{\mathrm{d}\vec{\imath}_2}{\mathrm{d}t}=\vec{\omega}_{21}\times\vec{\imath}_2=\omega\vec{\jmath}_2</math></center>

La aceleración de B cumple

<center><math>\vec{a}^B_{21}=\vec{a}^A_{21}+\alpha_{21}\vec{k}\times\overrightarrow{AB}-\omega_{21}^2\overrightarrow{AB}=
\dot{v}\vec{\imath}_2+v\omega\vec{\jmath}_2+\dot{\omega}b\vec{\jmath}_2-\omega^2b\vec{\imath}_2</math></center>

Agrupando términos queda

<center><math>\vec{a}^B_{21}=(\dot{v}-\omega^2b)\vec{\imath}_2+(\omega v+\dot{\omega}b)\vec{\jmath}_2</math></center>

===Ecuaciones de movimiento===
La única fuerza que actúa sobre la masa B es la tensión de la varilla, por lo que

<center><math>m(\dot{v}-\omega^2b)\vec{\imath}_2+m(\omega v+\dot{\omega}b)\vec{\jmath}_2=-F_T\vec{\imath}_2</math></center>

Igualando componente a componete

<center><math>m(\dot{v}-\omega^2b) = -F_T \qquad\qquad m(\omega v+\dot{\omega}b)=0</math></center>

Sobre la masa A actúan la tensión y la fuerza normal a la cuchilla

<center><math>m\dot{v}\vec{\imath}_2+mv\omega\vec{\jmath}_2=F_T\vec{\imath}_2+F_n\vec{\jmath}_2</math></center>

Igualando componente a componente

<center><math>m\dot{v}=F_T\qquad\qquad m\omega v = F_n</math></center>

Si eliminamos la tensión entre la ecuación para A y la ecuación para B llegamos a

<center><math>m(2\dot{v}-\omega^2b)=0</math></center>

que junto con la ecuación

<center><math>m(\omega v+\dot{\omega}b)=0</math></center>

constituyen las ecuaciones de movimiento del sistema.

==Sistema de ecuaciones de primer orden==
El sistema anterior posee solución analítica, aunque es laboriosa. Una forma más práctica se abordarlo sería mediante una integración numérica. Para ésta, conviene que no aparezcan funciones trigonométricas, que son costosas computacionalmente.

La solución alternativa nos da directamente un sistema de ecuaciones de primer orden, al que también se puede llegar con la primera forma de obtener las ecuaciones de movimiento.

Despejando, tenemos el sistema de ecuaciones

<center><math>\begin{array}{rcl}
\dot{v}&=&\dfrac{\omega^2b}{2}\\ && \\ \dot{\omega}&=&-\dfrac{\omega v}{b} \end{array}
</math></center>

Este sistema nos permite hallar v y &omega;, dadas las condiciones iniciales.

Una vez que hemos integrado este sistema, la variación del ángulo de la varilla con el tiempo resulta de integrar

<center><math>\dot{\theta}=\omega</math></center>

y, una vez que hayamos calculado &theta; hallamos x e y integrando las ecuaciones

<center><math>\begin{array}{rcl}
\dot{x}&=&v\cos(\theta)\\ \dot{y}&=&v\,\mathrm{sen}(\theta) \end{array}
</math></center>

Reuniendo todas las ecuaciones queda el sistema

<center><math>\begin{array}{rcl}
\dot{v}&=&\dfrac{\omega^2b}{2}\\ && \\ \dot{\omega}&=&-\dfrac{\omega v}{b} \\
\dot{\theta}&=&\omega\\ 
\dot{x}&=&v\cos(\theta)\\  \dot{y}&=&v\,\mathrm{sen}(\theta)
\end{array}
</math></center>

Los valores iniciales de estas variables salen de las posiciones y velocidades iniciales de las dos partículas, de las cuales se pueden despejar estas 5 variables.

Podemos eliminar completamente las funciones trigonométricas, añadiendo una variable más. Si consideramos S y C como variables independientes se cumple

<center><math>\dot{C}=-\omega S\qquad\qquad \dot{S}=\omega C</math></center>

lo que daría el sistema de 6 ecuaciones
<center><math>\begin{array}{rcl}
\dot{v}&=&\dfrac{\omega^2b}{2}\\ && \\ \dot{\omega}&=&-\dfrac{\omega v}{b} \\
\dot{x}&=&vC\\  
\dot{y}&=&vS\\
\dot{C}&=&-\omega S\\
\dot{S}&=&\omega C
\end{array}
</math></center>

que estaría listo para implementarlo en un programa de ordenador.

[[Categoría:Problemas de dinámica vectorial (CMR)]]