Ver el código de «Rapidez y tensión de un péndulo»

==Enunciado==
Empleando la ley de conservación de la energía, determine la velocidad con la que un péndulo simple de masa <math>m</math> y longitud <math>\ell</math> pasa por su punto más bajo, como función del ángulo máximo <math>\theta_0</math> con el que se separa de la vertical.

Compare este resultado con el que se obtiene empleando la aproximación lineal. Determine el error relativo cometido con esta aproximación para <math>\theta_0=10^\circ</math>, <math>\theta_0=20^\circ</math>,&hellip; <math>\theta_0=90^\circ</math>

Determine la tensión de la cuerda en el punto más bajo y en el punto de máxima separación de la vertical. en función del ángulo <math>\theta_0</math>

==Valor exacto==
[[Archivo:pendulo-Fvatan.png|300px|left]] 

La lenteja del péndulo está sometida a dos fuerzas: su peso y la tensión de la cuerda.

La tensión va siempre a lo largo de la propia cuerda y es por tanto perpendicular al desplazamiento en todo momento. Esto quiere decir que no realiza trabajo y no interviene en la ley de conservación de la energía.

La única energía potencial que debemos considerar es entonces la debida al peso. Tomando como origen de potencial el punto más bajo del péndulo, el valor de la energía potencial para una separación arbitraria <math>\theta</math> vale

<center><math>U = mgz = mg\ell(1-\cos(\theta))\,</math></center>

Cuando consideramos oscilaciones del péndulo, el punto de máxima amplitud es uno de velocidad nula (alcanza un estado de reposo instantáneo a partir del cual se da la vuelta). Igualando la energía mecánica de este punto a la que tiene en elpunto mas bajo obtenemos la igualdad

<center><math>E = 0 + mg\ell(1-\cos\theta_0)=\frac{1}{2}m |\vec{v}|_\mathrm{max}^2+ 0\,</math></center>

Obtenemos la velocidad máxima

<center><math>|\vec{v}|_\mathrm{max} = \sqrt{2g\ell(1-\cos\theta_0)}</math></center>

Esta expresión se puede simplificar un poco observando que

<center><math>\frac{1-\cos(\theta_0)}{2} = \mathrm{sen}^2\left(\frac{\theta_0}{2}\right)</math></center>

lo que da

<center><math>|\vec{v}|_\mathrm{max} = 2\sqrt{g\ell}\,\mathrm{sen}\left(\frac{\theta_0}{2}\right)</math></center>

==Valor aproximado==
Un péndulo obedece la ecuación de movimiento

<center><math>\frac{\mathrm{d}^2\theta}{\mathrm{d}t^2}=-\frac{g}{\ell}\,\mathrm{sen}\,\theta</math></center>

siendo <math>\theta</math> la inclinación respecto a la vertical (medida en radianes). Cuando esta separación es pequeña, se puede usar la aproximación lineal

<center><math>\mathrm{sen}\,\theta\simeq\theta \qquad\theta\ll 1</math></center>

lo que reduce la ecuación del péndulo a la de un oscilador armónico

<center><math>\frac{\mathrm{d}^2\theta}{\mathrm{d}t^2}\simeq-\frac{g}{\ell}\theta=-\omega^2\theta</math></center>

Cuando parte del reposo, desde una cierta separación <math>\theta_0</math>, el ángulo sigue una ley cosenoidal

<center><math>\theta = \theta_0\cos(\omega t)\,\qquad\omega = \sqrt{\frac{g}{\ell}}</math></center>

La velocidad lineal de la lenteja del péndulo es

<center><math>v = l\frac{\mathrm{d}\theta}{\mathrm{d}t} = -l\omega\theta_0\,\mathrm{sen}(\omega t)</math></center>

El valor máximo (en módulo) de esta velocidad lo alcanza en el momento en que se encuentra en el punto más bajo

<center><math>|\vec{v}|_\mathrm{max}=l\omega\theta_0= \sqrt{g\ell}\,\theta_0</math></center>

==Comparación==
El cociente entre el valor aproximado y el exacto es

<center><math>\frac{|\vec{v}|_\mathrm{ap}}{|\vec{v}|_\mathrm{ex}} = \frac{\theta_0}{\sqrt{2(1-\cos\theta_0)}} = \frac{\theta_0/2}{\mathrm{sen}\,(\theta_0/2)}</math></center>

donde hemos empleado la fórmula del ángulo mitad. Dado que

<center><math>\mathrm{sen}\,(x) < x\qquad \forall x > 0</math></center>

esto quiere decir que la aproximación de oscilador armónico predice una velocidad mayor que la real. El periodo calculado con esta aproximación será entonces más pequeño que el exacto.

<center>[[Archivo:velocidades-pendulo.png]]{{qquad}}{{qquad}}[[Archivo:error-pendulo.png]]</center>

El error relativo cometido en la aproximación es

<center><math>\epsilon = \frac{\left||\vec{v}|_\mathrm{ap}-|\vec{v}|_\mathrm{ex}\right|}{|\vec{v}|_\mathrm{ex}} = \frac{\theta_0/2-\,\mathrm{sen}\,(\theta_0/2)}{\mathrm{sen}\,(\theta/2)}</math></center>

Aplicando esta fórmula a los ángulos del enunciado

<center>
{| class="bordeado"
|-
! <math>\theta_0</math> (&deg;)
! <math>\theta_0</math> (rad)
! <math>\epsilon\ (\%)</math>
! <math>\theta_0</math> (&deg;)
! <math>\theta_0</math> (rad)
! <math>\epsilon\ (\%)</math>
|-
| 10
| 0.1745
| 0.127
| 60
| 1.047
| 4.720
|-
| 20
| 0.3491
| 0.510
| 70
| 1.222
| 6.501
|-
| 30
| 0.5236
| 1.152
| 80
| 1.396
| 8.610
|-
| 40
| 0.6981
| 2.060
| 90
| 1.571
| 11.07
|-
| 50
| 0.8727
| 3.245
|
|
|
|}
</center>

Vemos que, en general la aproximación es bastante buena y que incluso para ángulos tan grandes como 60&deg; el error es inferior al 5&thinsp;%.

==Tensión==
===Punto inferior===
[[Archivo:pendulo-tension-abajo.png|right]]

En el punto más bajo las dos fuerzas son normales a la trayectoria y por tanto la aceleración tangencial es nula (por ello la rapidez tiene un máximo). La partícula, no obstante, está acelerada, pues posee aceleración normal. La segunda ley de Newton nos da

<center><math>m\vec{a} = \vec{F}_T+m\vec{g}</math></center>

siendo <math>\vec{F}_T</math> la fuerza de reacción vincular, en este caso, la tensión de la cuerda. Despejando de aquí

<center><math>\vec{F}_T = m\vec{a}-m\vec{g}</math></center>

El término debido al peso vale

<center><math>-m\vec{g}=mg\vec{k}</math></center>

mientras que el debido a la aceleración normal es

<center><math>m\vec{a}=m\frac{|\vec{v}|^2}{R}\vec{N} = m\frac{4gl\mathrm{sen}^2(\theta_0/2)}{l}\vec{k} = 4mg\,\mathrm{sen}^2\left(\frac{\theta_0}{2}\right)\vec{k}</math></center>

lo que nos da la tensión

<center><math>\vec{F}_T= mg\left(1+4\,\mathrm{sen}^2\left(\frac{\theta}{2}\right)\right)\vec{k} = mg(3-2\cos(\theta_0))\vec{k}</math></center>

Esta tensión es siempre superior al peso. Sólo iguala al peso cuando la lenteja se encuentra en equilibrio en el punto inferior (<math>\theta_0=0</math>). 

El valor máximo de la tensión se da para el caso que la partícula parta del extremo superior (<math>\theta_0=\pi</math>), en cuyo caso iguala a 5 veces el peso. Este resultado es coincidente con el que se obtiene para una [[Anilla ensartada en aro|anilla ensartada en un aro]], ya que el movimiento de dicha anilla es el mismo que el de un péndulo (la condición de vínculo es la misma: distancia constante al centro de la circunferencia).

===Ángulo máximo===
[[Archivo:pendulo-tension-arriba.png|right]]

En la posición de máxima separación de la vertical la aceleración normal es nula, por serlo la celeridad. No así la aceleración tangencial. Puesto que la fuerza de reacción vincular es puramente normal a la trayectoria tenemos las ecuaciones

<center><math>m\vec{a}=m\vec{g}+\vec{F}_T\quad\Rightarrow\quad\begin{cases}m\vec{a}_t = m\vec{a} = m\vec{g}_\parallel & \\ m\vec{a}_n = \vec{0} = m\vec{g}_\perp + \vec{F}_T & \end{cases}</math></center>

siendo <math>\vec{g}_\parallel</math> y <math>\vec{g}_\perp</math> las componentes de la gravedad en la dirección tangente a la trayectoria y en la perpendicular a ella. 

La tensión de la cuerda compensa entonces a la componente perpendicular de la gravedad
<center>
<math>\vec{F}_T=-m\vec{g}_\perp</math></center>

y, en módulo,

<center><math>F_T = |m\vec{g}_\perp|=mg\left|\cos(\theta_0)\right|</math></center>

En forma vectorial multiplicamos por el vector unitario en la dirección de la cuerda y hacia adentro

<center><math>\vec{F}_T=mg\cos(\theta_0)\left(-\mathrm{sen}(\theta_0)\vec{\imath}+\cos(\theta_0)\vec{k}\right)</math></center>

<center>[[Archivo:tension-funcion-angulo.png]]</center>

En módulo, esta tensión es siempre menor que el peso, salvo en el caso en que la partícula esté en equilibrio en el punto más bajo, en que se reduce al caso anterior.

La tensión se anula si la masa alcanza una posición horizontal (<math>\theta_0=\pi/2</math>). Para ángulos máximos superiores, la tensión debe ir dirigida ''hacia afuera'' lo que es posible si se trata de una barra rígida (vínculo bilateral), pero no si se trata de una cuerda flexible (vínculo unilateral). Para una barra rígida y un ángulo <math>\theta_0=\pi</math> (extremo superior), la tensión vuelve a igualar al peso en módulo. Este caso está de nuevo de acuerdo con el sistema de una [[Anilla ensartada en aro|anilla ensartada en un aro]].

==Curva de potencial==
La lenteja del péndulo está sometida a dos fuerzas: el peso y la tensión de la barra. 
<center>[[Archivo:Pendulo-Fvatan.gif]]</center>

La tensión es siempre perpendicular al desplazamiento (ya que el desplazamiento es a lo largo de una circunferencia, mientras que la tensión es radial), por ello no realiza ningún trabajo y no influye en la conservación de la energía.

Por tanto, la única energía potencial es la gravitatoria

<center><math>U = mgz\,</math></center>

La altura z la podemos poner en función del ángulo que el péndulo forma con la vertical. Si tomamos como referencia de alturas el punto más bajo del péndulo, la altura en el momento que forma un ángulo \theta con la vertical es

<center><math>z = \ell - \ell\cos(\theta)\,</math></center>

por lo que la energía potencial puede escribirse en la forma

<center><math>U(\theta) = mg\ell(1-\cos(\theta))\,</math></center>
La gráfica de esta función es una sinusoide con el mínimo en <math>\theta=0</math>. En principio, la gráfica se extiende solo al rango <math>[-\pi,\pi]</math>. Sin embargo, para incluir la posibilidad de que el péndulo pueda superar la posición del máximo y dar vueltas completas, extendemos la gráfica a valores mayores de <math>\theta</math> por ambos extremos.

<center>[[Archivo:energia-pendulo-02.png]]</center>

Una forma más correcta, pero más difícil de ver, sería &ldquo;enrollar&rdquo; la gráfica de forma que el ángulo &minus;&pi; coincida con &pi;. En esta gráfica, la curva de potencial tiene forma de una elipse trazada sobre un cilindro.

<center>[[Archivo:energia-pendulo-03.png]]</center>

==Puntos de equilibrio==
Los puntos de equilibrio son los extremos de la curva de potencial. Estos son dos:

* <math>\theta = 0</math> en el cual la energía tiene un mínimo <math>U=0</math> y por tanto es un punto de equilibrio estable.
* <math>\theta = \pm \pi</math> (los dos valores son el mismo punto) en el cual la energía tiene un máximo <math>U=2mgL</math> y por tanto es un punto de equilibrio estable.

{| class="bordeado"
|-
| [[Archivo:pendulo-estable.png]]
| [[Archivo:pendulo-inestable.png]]
|-
! Estable
! Inestable
|-
|}

Alrededor del mínimo, se puede desarrollar la energía potencial en serie de Taylor

<center><math>U(\theta) = U(0) + \theta U'(0) + \frac{1}{2}\theta^2 U''(0)+\cdots</math></center>

siendo

<center><math>U(0) = mg\ell(1-\cos(0)) = 0\qquad\qquad U'(0) = mg\ell\,\mathrm{sen}(0) = 0\qquad\qquad U''(0) = mg\ell\cos(0)=mgL</math></center>

lo que permite aproximar la energía por una parábola

<center><math>U(\theta)\simeq \frac{1}{2}mg\ell\theta^2</math></center>

==Casos particulares==
Cuando la energía mecánica vale 

<center><math>E=mg\ell\,</math></center>

la energía es menor que la del máximo. Por tanto el movimiento es acotado, con dos puntos de retorno. Físicamente esto implica que el péndulo oscila periódicamente.

En el segundo caso, 

<center><math>E=3mg\ell\,</math></center>

la energía es mayor que la del máximo, por lo que el movimiento es no acotado. Físicamente, esto quiere decir que el péndulo llega hasta el punto más alto con una cierta velocidad, y continúa su marcha, describiendo vueltas completas, más rápidas en la parte inferior y más lentas en la superior.

<center>[[Archivo:energia-pendulo-04.png]]{{qquad}}{{qquad}}[[Archivo:energia-pendulo-05.png]]</center>