Ver el código de «Percusión sobre una barra articulada con muelle (MR-GIC)»

= Enunciado =
[[Archivo:Barra muelle articulada.png|sinmarco|derecha]]
El mecanismo de la figura está formado por una varilla delgada <math>OA</math> (sólido "2"), de masa <math>m</math> y longitud <math>L</math>, y un resorte ideal de constante recuperadora <math>k</math> y longitud natural nula. El extremo <math>O</math> de la varilla está unido mediante una rótula ideal al origen de un sistema de referencia fijo <math>OX_1Y_1Z_1</math> (sólido "1"). El otro extremo <math>A</math> de la varilla está conectado mediante el resorte a un pasador <math>C</math> de masa despreciable que puede deslizar libremente y sin rozamiento por el eje vertical <math>OZ_1</math>. En todo momento la orientación del eje del resorte es perpendicular a <math>OZ_1</math>. Todos los vínculos son lisos.  En el instante inicial <math>t=0</math>, el sistema se halla en reposo en la posición <math>\theta(0)=\pi/4</math>, <math>\phi(0)=0</math>. Recibe entonces una percusión <math>\vec{\hat{F}}=[\hat{F}_x,\, \hat{F}_y,\,0]_1</math> en el extremo <math>A</math>. Justo después de la percusión tenemos <math>\dot{\theta}(0^+)=3\Omega/\sqrt{2}</math>, <math>\dot{\phi}(0^+)=3\sqrt{2}\Omega</math>, donde <math>\Omega</math> es una constante conocida.  Calcula el valor de la percusión.

== Solución ==

== Cinemática y función de Lagrange ==
Este sistema es similar al que se estudio en [[Barra rotando con muelle horizontal (MRGIC) | barra rotando con muelle horizontal]], en el caso de las dos coordenadas <math>\{\theta, \phi \} </math>. Recuperamos de ese problema las expresiones de la reducción cinemática 
<center>
<math>
\vec{\omega}_{21} = \dot{\theta}\,\vec{\jmath}_0 + \dot{\phi}\,\vec{k}_0,
\qquad
\vec{v}^{\,O}_{21} = \vec{0}.</math>
</center>
y la función de Lagrange para el caso de las dos coordenadas libres
<center>
<math>
L = T - (U_g + U_m) = \dfrac{1}{2}I\,(\dot{\phi}^2\,\mathrm{sen}^2\theta + \dot{\theta}^2) - \dfrac{1}{2}mgL\cos\theta
- \dfrac{1}{2}kL^2\,\mathrm{sen}^2\theta.
</math>
</center>
Recordemos que <math> I = mL^2/3 </math>.

== Percusión ==
Como hay dos coordenadas independientes, tenemos dos ecuaciones percusivas analíticas
<center>
<math>
\begin{array}{l}
\Delta p_{\theta} = \hat{Q}^{\,NC}_{\theta}, \\
\Delta p_{\phi} = \hat{Q}^{\,NC}_{\phi}.
\end{array}
</math>
</center>
Los momentos generalizados son
<center>
<math>
\begin{array}{l}
p_{\theta} = \dfrac{\partial L}{\partial \dot{\theta}} = I\dot{\theta}, \\
p_{\phi} = \dfrac{\partial L}{\partial \dot{\phi}} = I\dot{\phi}\,\mathrm{sen}^2\theta.
\end{array}
</math>
</center>
Las percusiones generalizadas <math>\hat{Q}^{\,NC}_{\theta}, \hat{Q}^{\,NC}_{\phi}   </math> provienen de la acción de la percusión activa 
<center>
<math>
\vec{\hat{F}} = \hat{F}_x\,\vec{\imath}_1 + \hat{F}_y\,\vec{\jmath}_1.
</math>
</center>
Esta percusión está expresada en los vectores del sólido "1". Hay que tener cuidado en que base se proyectan los vectores. En este problema, la percusión ocurre en el instante <math>t=0 </math>. En este instante, <math>\phi(0)=0 </math>, por lo que los ejes <math>OX_1Y_1Z_1 </math> y <math>OX_0Y_0Z_0 </math> coinciden. Sin embargo, a la hora de calcular velocidades y aceleraciones, es mejor usar la base que sea válida durante todo el tiempo. En este caso, es la base del sólido "0".

En la formulación analítica las percusiones vinculares no aparecen. El peso y el muelle no intervienen porque no son fuerzas percusivas, es decir, su valor está acotado. 

Necesitamos calcular <math>\vec{v}^{\,A}_{21} </math>
<center>
<math>
\begin{array}{ll}
\vec{v}^{\,A}_{21} & = \vec{v}^{\,O}_{21} + \vec{\omega}_{21}\times\overrightarrow{OA} = 
L\dot{\theta}\cos\theta\,\vec{\imath}_0 + L\dot{\theta}\,\mathrm{sen}\,\theta\,\vec{\jmath}_0 -L\dot{\theta}\mathrm{sen}\,\theta\,\vec{k}_0. \\
& \overrightarrow{OA} = L\,\mathrm{sen}\,\theta\,\vec{\imath}_0 + L\cos\theta\,\vec{\jmath}_0.
\end{array}
</math>
</center>
Sus derivadas son
<center>
<math>
\begin{array}{l}
\dfrac{\partial \vec{v}^{\,A}_{21}}{\partial \dot{\theta}} = L\cos\theta\,\vec{\imath}_0 - L\,\mathrm{sen}\,\theta\,\vec{k}_0, \\
\dfrac{\partial \vec{v}^{\,A}_{21}}{\partial \dot{\phi}} = L\,\mathrm{sen}\,\theta\,\vec{\jmath}_0.
\end{array}
</math>
</center>
Las percusiones son
<center>
<math>
\begin{array}{l}
\hat{Q}^{\,NC}_{\theta} = \left.\vec{\hat{F}}\cdot\dfrac{\partial \vec{v}^{\,A}_{21}}{\partial \dot{\theta}}\right|_{\theta=\pi/4, \phi=0} = \dfrac{1}{\sqrt{2}}L\hat{F}_x, \\
\hat{Q}^{\,NC}_{\phi} = \left.\vec{\hat{F}}\cdot\dfrac{\partial \vec{v}^{\,A}_{21}}{\partial \dot{\phi}}\right|_{\theta=\pi/4, \phi=0} = \dfrac{1}{\sqrt{2}}L\hat{F}_y.
\end{array}
</math>
</center>
Recordemos que durante la percusión el valor de las coordenadas no cambia. 

Como se ha comentado mas arriba, durante la percusión los ejes de los sólidos "0" y "1" coinciden. Entonces
podemos escribir
<center>
<math>
\vec{\hat{F}} = \hat{F}_x\,\vec{\imath}_0 + \hat{F}_y\,\vec{\jmath}_0. 
</math>
</center>
Aplicando las ecuaciones percusivas analíticas tenemos
<center>
<math>
\begin{array}{l}
\Delta p_{\theta} = \dfrac{1}{\sqrt{2}}L\hat{F}_x, \\
\Delta p_{\phi} = \dfrac{1}{\sqrt{2}}L\hat{F}_y 
\end{array}
</math>
</center>
El problema dice que antes de la percusión el sólido estaba en reposo. Por tanto
<center>
<math>
p_{\theta}(0^-) = 0, \qquad p_{\phi}(0^-) = 0.
</math>
</center>
Después de la percusión tenemos
<center>
<math>
\begin{array}{l}
p_{\theta}(0^+) = I\dot{\theta}(0^+) = \dfrac{3}{\sqrt{2}}I\Omega, \\
p_{\phi}(0^+) = I\dot{\phi}(0^+)\,\mathrm{sen}^2(\theta(0)) = 3I\Omega.
\end{array}
</math>
</center>
Recordemos que durante la percusión el valor de las coordenadas no cambia.

En este problema los datos son el estado antes y después de la percusión, y las incógnitas son las componentes de la percusión. Es decir
<center>
<math>
\begin{array}{l}
\hat{F}_x = \dfrac{\sqrt{2}}{L} \Delta p_{\theta} = \dfrac{\sqrt{2}}{L}\,(p_{\theta}(0^+) - p_{\theta}(0^-)) = \dfrac{3I\Omega}{L} = mL\Omega, \\
\hat{F}_y = \dfrac{\sqrt{2}}{L} \Delta p_{\phi} = \dfrac{\sqrt{2}}{L}\,(p_{\phi}(0^+) - p_{\phi}(0^-)) = \dfrac{3I\Omega}{L} = mL\Omega.
\end{array}
</math>
</center>

== Resolución con ténicas vectoriales ==

La figura de la derecha muestra el digrama de percusiones libres. El peso y el muelle no intervienen en la percusión pues no son fuerzas percusivas. La percusión <math>\vec{\hat{O}}_{21} </math> garantiza que el extremo <math>O </math> de la barra no se mueve respecto al origen <math>O </math> del sólido "1". 

Observando la expresión de <math>\vec{\omega}_{21} </math>, podría parecer que, dado que no tiene componente en <math>\vec{\imath}_0 </math>, hay que añadir una componente de momento vincular en la dirección de <math>X_0 </math>. Sin embargo, no es así porque el sólido se modela con una barra de grosor nulo. Por tanto, sólo tiene dos grados de lbiertad de rotación, que sí están representados en <math>\vec{\omega}_{21} </math>. No hay que añdir ningún momento vincular.

Las expresiones de las percusiones son
<center>
<math>
\begin{array}{l}
\vec{\hat{F}} = \hat{F}_x\,\vec{\imath}_{0,1} + \hat{F}_y\,\vec{\jmath}_{0,1},\\
\vec{\hat{O}}_{21} = \hat{O}_x\,\vec{\imath}_{0,1} + \hat{O}_y\,\vec{\jmath}_{0,1}.
\end{array}
</math>
</center>
Hemos usado que, durante la percusión, los vectores de la base "0" y "1" coinciden.

Aplicamos ahora los teoremas fundamentales percusivos. Del T.C.M. percusivo (T.C.M.P.) tenemos
<center>
<math>
\Delta\vec{C} = \vec{\hat{F}} + \vec{\hat{O}}_{21}.
</math>
</center>
La cantidad de movimiento es
<center>
<math>
\begin{array}{ll}
\vec{C} & = m\vec{v}^{\,G}_{21} = m\,(\vec{v}^{\,O}_{21} + \vec{\omega}_{21}\times\overrightarrow{OG})
=
\dfrac{1}{2}mL\dot{\theta}\cos\theta\,\vec{\jmath}_0 + \dfrac{1}{2}mL\dot{\phi}\,\mathrm{sen}\,\theta\,\vec{\jmath}_0 - \dfrac{1}{2}mL\dot{\theta}\,\mathrm{sen}\,\theta\,\vec{k}_0.\\
& \overrightarrow{OG} = \dfrac{1}{2}L\,\mathrm{sen}\,\theta\,\vec{\imath}_0 + \dfrac{1}{2}L\cos\theta\,\vec{k}_0.
\end{array}
</math>
</center>
Es importante calcular las velocidades de modo que sean válidas para cualquier instante, no sólo durante la percusión.

La variación de la cantidad de movimiento durante la percusión es
<center>
<math>
\Delta\vec{C} = \left(\vec{C}^+ - \vec{C}^-\right)_{\theta=\pi/4, \phi=0}
=
\dfrac{mL}{2\sqrt{2}}\,\Delta\dot{\theta}\,\vec{\imath}_{0,1} + \dfrac{mL}{2\sqrt{2}}\Delta\dot{\phi}\,\vec{\jmath}_{0,1} - \dfrac{mL}{2\sqrt{2}}\Delta\dot{\theta}\,\vec{k}_{0,1}.
</math>
</center>
Aquí hemos vuelto a usar que durante la percusión los vectores de las bases "0" y "1" coinciden.

Aplicando el T.C.M.P. obtenemos tres ecuaciones
<center>
<math>
\begin{array}{lclr}
X_0) & \to & \dfrac{mL}{2\sqrt{2}}\Delta\dot{\theta} = \hat{O}_x + \hat{F}_x, & (1)\\
Y_0) & \to & \dfrac{mL}{2\sqrt{2}}\Delta\dot{\phi} = \hat{O}_y + \hat{F}_y, & (2), \\
Z_0) & \to & -\dfrac{mL}{2\sqrt{2}}\Delta\dot{\theta} = \hat{O}_z. & (3)
\end{array}
</math>
</center>

Ahora aplicamos el T.M.C. Percusivo (T.M.C.P.). Vamos a aplicarlo en el punto <math> O</math>. La única percusión que crea momento es <math>\vec{\hat{F}} </math>
<center>
<math>
\Delta\vec{L}_O = \overrightarrow{OA}\times\vec{\hat{F}}.
</math>
</center>
El momento angular es
<center>
<math>
\vec{L}_O = \overset\leftrightarrow{I}_O\cdot\vec{\omega}_{21}
</math>
</center>
El tensor de inercia debe ser expresado en ejes solidarios con el sólido "2". Vamos a utilizar los ejes "2" que usamos en [[Barra articulada rotando en el espacio(MR G.I.C.)| barra articulada rotando en el espacio]]. Para poder hacer el producto el tensor de inercia y el vector rotación tienen que estar expresados en la misma base. Tenemos
<center>
<math>
\begin{array}{l}
\vec{\imath}_0 = \cos\theta\,\vec{\imath}_2 + \mathrm{sen}\,\theta\,\vec{k}_2,\\
\vec{\jmath}_0 = \vec{\jmath}_2,\\
\vec{k}_0 = -\mathrm{sen}\,\theta\,\vec{\imath}_2 + \cos\theta\,\vec{k}_2.
\end{array}
</math>
</center>
con lo que el vector rotación expresado en la base "2" queda
<center>
<math>
\vec{\omega}_{21} =
\left[
\begin{array}{c}
-\dot{\phi}\,\mathrm{sen}\,\theta\\
\dot{\theta}\\
\dot{\phi}\cos\theta
\end{array}
\right]_2
</math>
</center>
El momento angular en <math>O </math> es
<center>
<math>
\vec{L}_O = \overset\leftrightarrow{I}_O\cdot\vec{\omega}_{21}
=
I
\left[
\begin{array}{ccc}
1  & 0 & 0 \\
0 & 1 & 0 \\
0 & 0 & 0
\end{array}
\right]_2
\left[
\begin{array}{c}
-\dot{\phi}\,\mathrm{sen}\,\theta\\
\dot{\theta}\\
\dot{\phi}\cos\theta
\end{array}
\right]_2
=
\left[
\begin{array}{c}
-I\dot{\phi}\,\mathrm{sen}\,\theta\\
I\dot{\theta}\\
 0
\end{array}
\right]_2
</math>
</center>
con <math>I = mL^2/3 </math>.

Podemos expresar el momento angular en la base "0" o, alternativamente, proyectar la percusión <math>\vec{\hat{F}} </math> en la base "2" 
<center>
<math>
\vec{\hat{F}} = \hat{F}_x\cos\theta\,\vec{\imath}_2 + \hat{F}_y\,\vec{\jmath}_2 + \hat{F}_x\,\mathrm{sen}\,\theta\,\vec{k}_2.
</math>
</center>
La variación del momento angular es
<center>
<math>
\Delta\vec{L}_O = \left(\vec{L}_O^+ - \vec{L}_O^-\right)_{\theta=\pi/4, \phi=0}
=
-\dfrac{I}{\sqrt{2}}\,\Delta\dot{\phi}\,\vec{\imath}_{2} + I\Delta\dot{\theta}\,\vec{\jmath}_{2}.
</math>
</center>
Y el momento percusivo es
<center>
<math>
\overrightarrow{OA}\times\vec{\hat{F}} =
\left|
\begin{array}{ccc}
\vec{\imath}_2 & \vec{\jmath}_2 & \vec{k}_2 \\
0 & 0 & L\\
\hat{F}_x\cos\theta & \hat{F}_y & \hat{F}_x\,\mathrm{sen}\,\theta
\end{array}
\right|_2
=
-L\hat{F}_y\,\vec{\imath}_2 + L\hat{F}_x\cos\theta\,\vec{\jmath}_2.
</math>
</center>
En el instante de la percusión es
<center>
<math>
\left.\overrightarrow{OA}\times\vec{\hat{F}}\right|_{\theta=\pi/4, \phi=0}
=
-L\hat{F}_y\,\vec{\imath}_2 + \dfrac{L\hat{F}_x}{\sqrt{2}}\,\vec{\jmath}_2.
</math>
</center>

Aplicando el T.M.C.P. obtenemos dos ecuaciones
<center>
<math>
\begin{array}{lclr}
X_2) & \to &\dfrac{I}{\sqrt{2}}\Delta\dot{\phi} = L\hat{F}_y, & (4)\\
Y_2) & \to &I\Delta\theta = \dfrac{L\hat{F}_x}{\sqrt{2}}. & (5)
\end{array}
</math>
</center>

Las incógnitas son <math>\{ \hat{F}_x, \hat{F}_y, \hat{O}_x, \hat{O}_y, \hat{O}_z \} </math>. Con las cinco ecuaciones obtenidas llegamos a 
<center>
<math>
\vec{\hat{F}} = mL\Omega\,(\vec{\imath}_{0,1} + \vec{\jmath}_{0,1}), 
\qquad
\vec{\hat{O}}_{21} = mL\Omega\left(-\dfrac{1}{4}\vec{\imath}_{0,1} + \dfrac{1}{2}\,\vec{\jmath}_{0,1} - \dfrac{3}{4}\,\vec{k}_{0,1}\right).
</math>
</center>















[[Categoría:Problemas de Dinámica Impulsiva]]