(Página creada con «==Enunciado== Se tiene un sistema horizontal en el que una partícula P, de masa m, se encuentra unida a una varilla de longitud <math>\ell</math> cuyo otro extremo, A, se halla articulado a una segunda varilla, de longitud b, cuyo segundo extremo, O, está fijo. La varilla OA gira en torno a O con velocidad angular constante Ω, mientras que la varilla AP puede girar libremente en torno a A. Sea θ el ángulo que AP forma con la prolongación de OA. <center>Ar…»)
 
Línea 1: Línea 1:
==Enunciado==
==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
Se tiene un sistema horizontal en el que una partícula P, de masa m, se encuentra unida a una varilla de longitud <math>\ell</math> cuyo otro extremo, A, se halla articulado a una segunda varilla, de longitud b, cuyo segundo extremo, O, está fijo. La varilla OA gira en torno a O con velocidad angular constante Ω, mientras que la varilla AP puede girar libremente en torno a A. Sea &theta; el ángulo que AP forma con la prolongación de OA.
# ¿Qué vínculos ligan las posiciones y velocidades de las partículas?
<center>[[Archivo:masa-varilla-articulada.png]]</center>
# ¿Hacia dónde van dirigidas las fuerzas de reacción vincular?
# ¿Qué vínculo hay entre las coordenadas cartesianas de P? Escríbalo en forma geométrica, cinemática y pfaffiana.
# 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.
# Obtenga la ecuación de movimiento para el ángulo &theta;.
# Introduciendo las variables adecuadas, reduzca este problema a un sistema de ecuaciones de primer orden.
# ¿Qué puntos de equilibrio hay para el ángulo &theta;? ¿Son estables o inestables?
<center>[[Archivo:Dos-masas-cuchilla.png|400px]]</center>
==Ecuación de vínculo==
La ecuación de vínculo sobre P es que su distancia al punto A es constante


==Vínculos==
<center><math>|\overrightarrow{AP}|=\ell=\mathrm{cte}</math></center>
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
Se trata de escribir esta relación en términos de las coordenadas cartesianas de P.


<center><math>|\overrightarrow{AB}|=b\qquad\Rightarrow\qquad(x_B-x_A)^2+(y_B-y_A)^2 = b^2\,</math></center>
La posición de A es en todo momento


* La velocidad de la partícula 1 es paralela al vector de posición relativa
<center><math>\overrightarrow{OA}=b\cos(\Omega t)\vec{\imath}_1+b\,\mathrm{sen}(\Omega t)\vec{\jmath}_1</math></center>


<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>
con lo que la ecuación de vínculo se puede escribir en la forma


: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.
<center><math>(x-b\cos(\Omega t))^2+(y-b\,\mathrm{sen}(\Omega t))^2 = \ell^2</math></center>


Si usamos como coordenadas las indicadas en el enunciado
Esta sería la forma geométrica.


<center><math>x_A = x\qquad\qquad y_A = y\qquad\qquad x_B=x+bC\qquad\qquad y_B=y+bS</math></center>
Obtenemos la forma cinemática derivando esta respecto al tiempo.


(con C y S representando al coseno y el seno de &theta;), el primer vínculo se vuelve trivial
<center><math>2(x-b\cos(\Omega t))(\dot{x}+b\Omega\,\mathrm{sen}(\Omega t))+2(y-b\,\mathrm{sen}(\Omega t))(\dot{y}-b\Omega\cos(\Omega t))=0</math></center>


<center><math>b=\mathrm{cte.}\,</math></center>
Simplificando por 2 y agrupando términos queda


y el segundo se escribe
<center><math>\dot{x}(x-b\cos(\Omega t))+\dot{y}(y-b\,\mathrm{sen}(\Omega t))+b\Omega(x\,\mathrm{sen}(\Omega t)-y\cos(\Omega t))=0</math></center>


<center><math>-S\dot{x}+C\dot{y}=0</math></center>
La forma pfaffiana se obteien multiplicando la condición cinemática por <math>\mathrm{d}t</math>
==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>\mathrm{d}x(x-b\cos(\Omega t))+\mathrm{d}y(y-b\,\mathrm{sen}(\Omega t))+b\Omega\,\mathrm{d}t(x\,\mathrm{sen}(\Omega t)-y\cos(\Omega t))=0</math></center>


<center><math>\vec{F}_{T1}=F_T(C\vec{\imath}+S\vec{\jmath})</math></center>
==Ecuación de movimiento==
Siguiendo los cálculos y la notación del problema &ldquo;[[Cinem%C3%A1tica_de_dos_barras_articuladas_(CMR)|Dos barras articuladas]]&rdquo; la posición, velocidad y aceleración de A son


:Para el segundo, es perpendicular a la dirección de movimiento posible y por tanto, perpendicular a la varilla
<center><math>\overrightarrow{OA}=b\vec{\imath}_2\qquad\qquad \vec{v}^A_{21}=b\Omega\vec{\jmath}_2\qquad\qquad \vec{a}^A_{21}=-b\Omega^2\vec{\imath}_2</math></center>


<center><math>\vec{F}_{n1}=F_n(-S\vec{\imath}+C\vec{\jmath})</math></center>
Por ser A una articulación entre los sólidos 2 y 3, estos valores corresponden también al movimiento {31}.


* 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.
Para el punto P tenemos la posición


<center><math>\vec{F}_{T2}=-F_T(C\vec{\imath}+S\vec{\jmath})</math></center>
<center><math>\overrightarrow{OP}=b\vec{\imath}_2+\ell\vec{\imath}_3</math></center>


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


<center><math>\begin{array}{rcl}
<center><math>\vec{v}^P_{31}=b\Omega\vec{\jmath}_2+\ell(\Omega+\dot{\theta})\vec{\jmath}_3</math></center>
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>
y la aceleración


<center><math>\begin{array}{rcl}
<center><math>\vec{a}^P_{31}=-b\Omega^2\vec{\imath}_2+\ell\ddot{\theta}\vec{\jmath}_3-\ell(\Omega+\dot{\theta})^2\vec{\jmath}_3</math></center>
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
La única fuerza que actúa sobre la masa en P es la tensión de la varilla, de manera que la segunda Ley de Newton queda


<center><math>\begin{array}{rcl}
<center><math>m\vec{a}^P_{31}=-F_{T}\vec{\imath}_3</math></center>
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
Sustituimos la aceleración de P


<center><math>mb\dot{\theta}^2 = 2F_T</math></center>
<center><math>-mb\Omega^2\vec{\imath}_2+m\ell\ddot{\theta}\vec{\jmath}_3-m\ell(\Omega+\dot{\theta})^2\vec{\jmath}_3=-F_{T}\vec{\imath}_3</math></center>


y si multiplicamos la primera por el seno, la segunda por el coseno y restamos
Para eliminar la tensión de los cálculos proyectamos sobre la dirección ortogonal, que en este caso es la de <math>\vec{\jmath}_3</math>. Queda


<center><math>mb\ddot{\theta}=-F_n</math></center>
<center><math>mb\Omega^2\,\mathrm{sen}(\theta)+m\ell\ddot{\theta}=0</math></center>


Esto nos permite eliminar las fuerzas de reacción vincular
o, equivalentemente,


<center><math>\begin{array}{rcl}
<center><math>\ddot{\theta}=-\frac{b\Omega^2}{\ell}\mathrm{sen}(\theta)</math></center>
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
Esta es la ecuación de un péndulo simple. Empleando fuerzas ficticias, diríamos que aquí el papel de la gravedad lo desempeña la fuerza centrífuga. Así es como se llegaría a este mismo resultado empleando un sistema de referencia en rotación.


<center><math>-S\dot{x}+C\dot{y}=0</math></center>
==Puntos de equilibrio==
Al ser la ecuación de movimiento equivalente a un péndulo, al análisis del equilibrio es idéntico.


Hay que remarcar que aunque tengamos tres variables, el sistema tiene solo dos grados de libertad.
* <math>\theta=0</math> representa un punto de equilibrio estable.
* <math>\theta=\pi</math> representa un punto de equilibrio inestable.


==Ecuaciones de movimiento. Forma alternativa==
Es decir, el sistema sería estable con la varilla AP completamente extendida y sería inestable con AP plegada sobre OA.
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)]]

Revisión del 22:13 26 nov 2023

Enunciado

Se tiene un sistema horizontal en el que una partícula P, de masa m, se encuentra unida a una varilla de longitud cuyo otro extremo, A, se halla articulado a una segunda varilla, de longitud b, cuyo segundo extremo, O, está fijo. La varilla OA gira en torno a O con velocidad angular constante Ω, mientras que la varilla AP puede girar libremente en torno a A. Sea θ el ángulo que AP forma con la prolongación de OA.

  1. ¿Qué vínculo hay entre las coordenadas cartesianas de P? Escríbalo en forma geométrica, cinemática y pfaffiana.
  2. Obtenga la ecuación de movimiento para el ángulo θ.
  3. ¿Qué puntos de equilibrio hay para el ángulo θ? ¿Son estables o inestables?

Ecuación de vínculo

La ecuación de vínculo sobre P es que su distancia al punto A es constante

Se trata de escribir esta relación en términos de las coordenadas cartesianas de P.

La posición de A es en todo momento

con lo que la ecuación de vínculo se puede escribir en la forma

Esta sería la forma geométrica.

Obtenemos la forma cinemática derivando esta respecto al tiempo.

Simplificando por 2 y agrupando términos queda

La forma pfaffiana se obteien multiplicando la condición cinemática por

Ecuación de movimiento

Siguiendo los cálculos y la notación del problema “Dos barras articuladas” la posición, velocidad y aceleración de A son

Por ser A una articulación entre los sólidos 2 y 3, estos valores corresponden también al movimiento {31}.

Para el punto P tenemos la posición

la velocidad

y la aceleración

La única fuerza que actúa sobre la masa en P es la tensión de la varilla, de manera que la segunda Ley de Newton queda

Sustituimos la aceleración de P

Para eliminar la tensión de los cálculos proyectamos sobre la dirección ortogonal, que en este caso es la de . Queda

o, equivalentemente,

Esta es la ecuación de un péndulo simple. Empleando fuerzas ficticias, diríamos que aquí el papel de la gravedad lo desempeña la fuerza centrífuga. Así es como se llegaría a este mismo resultado empleando un sistema de referencia en rotación.

Puntos de equilibrio

Al ser la ecuación de movimiento equivalente a un péndulo, al análisis del equilibrio es idéntico.

  • representa un punto de equilibrio estable.
  • representa un punto de equilibrio inestable.

Es decir, el sistema sería estable con la varilla AP completamente extendida y sería inestable con AP plegada sobre OA.