Partícula en anilla giratoria

De Laplace

Revisión a fecha de 20:41 2 dic 2017; Antonio (Discusión | contribuciones)

(dif) ← Revisión anterior | Revisión actual (dif) | Revisión siguiente → (dif)

Contenido

1 Enunciado
2 Ecuación del vínculo
3 Ecuaciones de movimiento
4 Puntos de equilibrio
- 4.1 Estabilidad del equilibrio
  - 4.1.1 A partir de la ecuación de movimiento
  - 4.1.2 A partir de la energía
5 Equilibrio con rozamiento

1 Enunciado

Una pequeña anilla de masa $m$ está ensartada sin rozamiento en un aro circular de masa $M$ y radio $b$ situado verticalmente que puede girar en torno a un diámetro vertical. Suponga que este aro se hace girar con velocidad angular constante Ω alrededor de este diámetro.

Escriba la ecuación del vínculo sobre la partícula. ¿De qué tipo es?
Determine la posición de los posibles puntos de equilibrio en la vertical, así como la estabilidad de éstos, en función del valor de Ω.
Suponga ahora que entre la anilla y el aro existe un coeficiente de rozamiento seco μ, ¿cómo queda en ese caso la ecuación de movimiento para la anilla? Para un valor de Ω dado, ¿cuál es en ese caso el rango de posiciones verticales de equilibrio?

2 Ecuación del vínculo

La partícula está obligada a moverse sobre una línea, la cual la podemos considerar como una intersección de dos superficies:

Una esfera de radio b

$x^2+y^2+z^2=b^2\,$

o, en cilíndricas,

$\rho^2+z^2 = b^2\,$

o, en esféricas,

$r=b\,$

Un plano giratorio. Este plano fija la coordenada θ de cilíndricas, en

$\theta=\Omega t\,$

(que sería $\varphi=\Omega t$ en esféricas). En cartesianas, sería

$\frac{y}{x}=\mathrm{tg}(\Omega t)\qquad\Rightarrow\qquad -x\,\mathrm{sen}(\Omega t)+y\cos(\Omega t)=0$

Este vínculo es liso, bilateral, reónomo (pues depende del tiempo) y geométrico. Al ser geométrico, es también holónomo.

3 Ecuaciones de movimiento

3.1 Coordenadas

Para describir el movimiento emplearemos la base de cilíndricas y parametrizamos la ecuación de la curva como

$\rho = b\,\mathrm{sen}(\phi)\qquad\qquad \theta=\Omega t\qquad\qquad z = -b\cos(\phi)$

siendo $φ$ el ángulo con la vertical descendente. En realidad, esto equivale a usar coordenadas esféricas, pero lo haremos así para aprovechar las fórmulas de coordenadas cilíndricas, más fáciles de recordar.

Empleando la abreviatura, que usaremos en lo sucesivo

$C=\cos(\phi)\qquad\qquad S=\mathrm{sen}(\phi)$

estas coordenadas se escriben

$\rho = b\,S\qquad\qquad \theta=\Omega t\qquad\qquad z = -b\,C$

Las derivadas respecto al tiempo valen, por la regla de la cadena

$\dot{\rho} = b\,\dot{\phi}\,C\qquad\qquad \dot{\theta}=\Omega \qquad\qquad \dot{z} = b\,\dot{\phi}\,S$

y las segundas derivadas

$\ddot{\rho} = b(\ddot{\phi}\,C-\dot{\phi}^2S)\qquad\qquad \ddot{\theta}=0 \qquad\qquad \ddot{z} = b(\ddot{\phi}\,S+\dot{\phi}^2\,C)$

3.2 Velocidad y aceleración

3.2.1 Velocidad

Sustituímos en la expresión de la velocidad en cilíndricas

$\vec{v}=\dot{\rho}\vec{u}_\rho+\rho\vec{u}_\theta+\dot{z}\vec{u}_z=b\left(\dot{\phi}\,C\vec{u}_\rho+\Omega S\vec{u}_\theta+\dot{\phi}\,S\vec{u}_z\right)$

3.2.2 Aceleración

De la misma manera obtenemos la aceleración. Si separamos por componentes queda, para la componente cilíndrica radial

$a_\rho=\ddot{\rho}-\rho\dot{\theta}^2=b(\ddot{\phi}\,C-\dot{\phi}^2S-\Omega^2 S)$

para la acimutal

$a_\theta=2\dot{\rho}\dot{\theta}+\rho\ddot{\theta}=2b\dot{\phi}\,\Omega\,C$

y para la vertical

$a_z=\ddot{z}=b(\ddot{\phi}\,S+\dot{\phi}^2\,C)$

3.3 Fuerzas

En ausencia de rozamiento, la partícula está sometida a dos fuerzas: el peso y la reacción del aro.

A su vez, la reacción la podemos descomponer en dos partes, cada una correspondiente a una ecuación del vínculo. La primera impide que la partícula se mueva radialmente, por lo que tendrá una dirección radial respecto al centro del aro

$\vec{F}_n=F_n(S\vec{u}_\rho-C\vec{u}_z)$

Se puede hallar esta dirección calculando el gradiente de la ecuación del vínculo y escribiendo un vector paralelo a este gradiente.

Para la segunda ecuación tenemos una fuerza acimutal que es la que obliga a la anilla a dar vueltas con el aro

$\vec{F}_\theta=F_\theta \vec{u}_\theta$

Las dos cantidades $F n$ y $F θ$ son también incógnitas que habrá que determinar.

3.4 Segunda ley de Newton

Sustituyendo la aceleración y las fuerzas en la segunda ley de Newton obtenemos las ecuaciones de movimiento siguientes. Para la componente cilíndrica radial

$mb(\ddot{\phi}\,C-\dot{\phi}^2S-\Omega^2 S)=F_nS$

para la acimutal

$2mb\dot{\phi}\,\Omega\,C=F_\theta$

y para la vertical

$mb(\ddot{\phi}\,S+\dot{\phi}^2\,C)=-mg-F_nC$

Si multiplicamos la primera por C y la tercera por S y sumamos eliminamos la fuerza normal

$mb(\ddot{\phi}\,C^2-\dot{\phi}^2SC-\Omega^2 SC)+mb(\ddot{\phi}\,S^2+\dot{\phi}^2\,CS)=-mgS$

Simplificamos aquí y queda

$mb(\ddot{\phi}-\Omega^2 SC)=-mgS$

o, equivalentemente

$\ddot{\phi}=-S\left(\frac{g}{b}-\Omega^2 C\right)$

que es la ecuación de movimiento para la coordenada φ. Podemos comprobar que, en ausencia de rotación, se reduce a la ecuación del péndulo simple

$(\Omega = 0)\qquad\qquad \ddot{\phi}=-\frac{g}{b}\mathrm{sen}(\phi)$

4 Puntos de equilibrio

===Posiciones de equilibrio===& Si el aro no girara, habría solo dos posiciones de equilibrio, una en el punto más bajo y otra en el punto más alto. La primera sería estable y la segunda inestable.

La rotación del aro introduce nuevas posiciones de equilibrio. Podemos visualizarlo imaginando la acción de una fuerza centrífuga que apunta hacia el exterior, posibilitando que la anilla se quede estacionaria a una cierta altura.

Hay que señalar que en ese caso, estas posiciones serían estacionarias respecto al aro, pero no respecto a un observador exterior, el cual vería la anilla dar vueltas alrededor del eje, como el resto del aro en que se halla ensartada.

Determinamos las posiciones de equilibrio a partir de la ecuación de movimiento. En el equilibrio la segunda derivada de $φ$ debe ser nula, ya que si no la partícula saldría de la posición. Esto nos da la ecuación

$-S\left(\frac{g}{b}-\Omega^2 C\right)=0\,$

que ofrece las siguientes posibilidades:

que se anule el primer factor

0 = S = sen(φ)

con dos soluciones

- $\phi=0\,$ . El punto más bajo del aro.
- $\phi=\pi\,$ . El punto más alto.
Que se anule el segundo factor

$\frac{g}{b}-\Omega^2C = 0\qquad\Rightarrow\qquad \cos(\phi)=\frac{g}{b\Omega^2}$

Esta ecuación no siempre tiene solución real.

- Si $\Omega < \sqrt{g/b}$ resultaría un coseno mayor que la unidad y por tanto no hay solución-
- Si $\Omega\geq \sqrt{g/b}$ aparecen dos nuevas soluciones:

$\phi=\pm \arccos\left(\frac{g}{\Omega^2 b}\right)$

Estas corresponden a posiciones simétricas alrededor del eje de rotación.

Hay que destacar que estas posiciones no son estacionarias. Cuando el aro gira, la anilla se mueve con él y por tanto su velocidad y su aceleración no son nulas.

4.1 Estabilidad del equilibrio

La estabilidad del equilibrio la podemos determinar a partir de la propia ecuación de movimiento o a partir de la energía.

4.1.1 A partir de la ecuación de movimiento

El procedimiento consiste en linealizar la ecuación de movimiento alrededor de la posición de equilibrio. Si resulta una ecuación de tipo del oscilador armónico

$\ddot{x}=-Kx$

la posición de equilibrio es estable. Si por el contrario queda de la forma

$\ddot{x}=+Kx$

la posición es repulsiva y el equilibrio es inestable.

Analizamos los casos que nos han aparecido:

Punto más bajo del aro

Cuando el ángulo $φ$ es muy pequeño podemos hacer la aproximación

$\mathrm{sen}(\phi)\simeq \phi\qquad\qquad \cos(\phi)\simeq 1$

lo que nos permite aproximar la ecuación de movimiento por

$\ddot{\phi}\simeq -\left(\frac{g}{b}-\Omega^2\right)\phi$

que nos da dos posibilidades:

Si $\Omega < \sqrt{g/b}$ este punto de equilibrio es estable.
Si $\Omega > \sqrt{g/b}$ el punto es inestable.

El caso de igualdad es lo que se conoce como una bifurcación, pasando de un estado al otro.

Punto más alto del aro

Cuando el ángulo $\phi\simeq \pi$ podemos hacer la sustitución

$\phi=\pi-\beta\qquad\qquad \beta \ll \pi$

con lo que se cumple

$\mathrm{sen}(\phi)=\mathrm{sen}(\pi-\beta)=\mathrm{sen}(\beta)\simeq \beta\qquad\qquad \cos(\phi)\simeq -1$

y

$\ddot{\beta}=-\ddot{\phi}\simeq\left(\frac{g}{b}+\Omega^2\right)\beta$

El coeficiente es siempre positivo, por lo que este punto es siempre inestable, como cabe esperar.

Puntos laterales: Suponiendo que Ω supera el valor crítico, tenemos los puntos definidos por la ecuación

$\cos(\phi_0)=\frac{g}{b\Omega^2}$

Si suponemos un ángulo próximo a este valor

$\phi=\phi_0+\gamma\qquad\qquad \gamma\ll \phi_0\,$

queda

$\begin{array}{rcl} \cos(\phi) &= & \cos(\phi_0)\cos(\gamma)-\mathrm{sen}(\phi_0)\mathrm{sen}(\gamma)\simeq \cos(\phi_0)-\mathrm{sen}(\phi_0)\gamma \\ \mathrm{sen}(\phi) &= & \mathrm{sen}(\phi_0)\cos(\gamma)+\mathrm{sen}(\phi_0)\cos(\gamma)\simeq \mathrm{sen}(\phi_0)+\cos{sen}(\phi_0)\gamma \end{array}$

Llevamos esto a la ecuación de movimiento

$\ddot{\gamma}=\ddot{\phi}\simeq-\left(\mathrm{sen}(\phi_0)+\cos{sen}(\phi_0)\gamma\right)\left(\frac{g}{b}-\Omega^2\cos(\phi_0)+\Omega^2\mathrm{sen}(\phi_0)\gamma \right)$

Ahora bien, por la ecuación que define a $φ 0$ se cancelan términos y esta expresión se reduce a

$\ddot{\gamma}\simeq -\mathrm{sen}^2(\phi_0)\Omega^2 \gamma=-\left(\Omega^2-\frac{g}{b}\right)\gamma$

El coeficiente es siempre negativo y por tanto los dos puntos simétricos son de equilibrio estable.

Resumiendo tenemos las siguientes posibilidades:

Si $\Omega < \sqrt{g/b}$ hay un punto estable en el punto más bajo del aro y uno inestable en el más alto.
Si $\Omega > \sqrt{g/b}$ los dos anteriores son inestables y aparecen dos estables situados simétricamente.

Partícula en anilla giratoria

De Laplace

Contenido

1 Enunciado

2 Ecuación del vínculo

3 Ecuaciones de movimiento

3.1 Coordenadas

3.2 Velocidad y aceleración

3.2.1 Velocidad

3.2.2 Aceleración

3.3 Fuerzas

3.4 Segunda ley de Newton

4 Puntos de equilibrio

4.1 Estabilidad del equilibrio

4.1.1 A partir de la ecuación de movimiento

4.1.2 A partir de la energía

5 Equilibrio con rozamiento

Herramientas:

Buscar

Vistas

Herramientas

Herramientas personales

Navegación

SEARCH

TOOLBOX

LANGUAGES