Ver el código de «Partícula en el interior de un tubo»

==Enunciado==
Una partícula de masa <math>m</math> se encuentra en el interior de un tubo estrecho, el cual se halla en todo momento contenido en el plano horizontal OXY girando con velocidad angular <math>\omega</math> constante alrededor del eje OZ, de forma que la posición de la partícula puede escribirse como

<center><math>x = \rho\,\mathrm{cos}(\omega t)\,\,</math>;{{qquad}}{{qquad}}<math>y= \rho\,\mathrm{sen}(\omega t)</math></center> 

donde <math>\rho = \rho(t)\,</math>, función que hay que determinar, define la posición de la partícula a lo largo del tubo.

<ol>
<li> Halle la ecuación diferencial que debe satisfacer <math>\rho(t)\,</math> sabiendo que el tubo no puede ejercer fuerza en la dirección longitudinal (no hay rozamiento).</li>
<li> Compruebe que <math>\rho(t) =A\mathrm{e}^{\omega t}\,</math> es una solución de dicha ecuación diferencial.</li>
<li>Para esta solución particular
<ol><li>Calcule la fuerza ejercida por el tubo en cada instante.</li>
<li> Halle la potencia desarrollada por el tubo sobre la partícula.</li>
<li>Calcule el trabajo realizado sobre la partícula durante el tiempo que emplea en pasar de <math>\rho = b\,</math> a <math>\rho = 2b\,</math>.
<li>Evalúe el incremento de energía cinética de la partícula en el mismo intervalo y compruebe que se verifica el teorema de la energía cinética.</li>
</ol>
</ol>

==Ecuación diferencial==
La expresión dada en el enunciado para las ecuaciones horarias cartesianas de la partícula nos invita a replantear el problema en coordenadas polares. Como se recordará, la relación cartesianas-polares viene dada por

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

Por tanto, de la comparación con el enunciado se deduce que las ecuaciones horarias polares de la partícula son

<center><math> \rho=\rho(t)\,</math>{{qquad}}{{qquad}}<math>\theta=\omega t \,</math></center>

donde se desconoce en principio la dependencia de la coordenada radial con respecto al tiempo <math>\rho(t)\,</math>, razón por la cual tiene sentido que este primer apartado del problema nos pida precisamente que deduzcamos la ecuación diferencial que debe satisfacer dicha función.

Recordemos las expresiones generales de los vectores de posición, velocidad y aceleración de una partícula en coordenadas polares, y sustituyamos en ellas los valores particulares de <math>\dot{\theta}\,</math> y <math>\ddot{\theta}\,</math> que se obtienen de derivar la función <math>\theta(t)=\omega t\,</math>, es decir,   <math>\dot{\theta}=\omega\,\,\,</math> y <math>\,\,\ddot{\theta}=0\,</math>

<math> \vec{r} = \rho\,\vec{u}_{\rho}\,\,;</math>

<math>
\vec{v} = \dot{\rho}\,\vec{u}_{\rho}+\rho\,\dot{\theta}\,\vec{u}_{\theta} = \dot{\rho}\,\vec{u}_{\rho}+\omega\rho\,\vec{u}_{\theta}\,\,;</math>

<math> \vec{a} = (\ddot{\rho}-\rho\,\dot{\theta}^{2})\,\vec{u}_{\rho}+(2\dot{\rho}\,\dot{\theta}+\rho\,\ddot{\theta})\,\vec{u}_{\theta} = (\ddot{\rho}-\omega^{2}\rho)\,\vec{u}_{\rho}+2\omega\dot{\rho}\,\vec{u}_{\theta}
</math>

Teniendo presente que la única fuerza que actúa sobre partícula es la fuerza de reacción vincular ejercida por el tubo, de la cual se nos indica expresamente que no tiene componente longitudinal o radial (es decir, <math>\vec{\Phi}=\Phi\,\vec{u}_{\theta}\,</math>), podemos proyectar la segunda ley de Newton en las direcciones radial y acimutal obteniendo 

<center><math>m\vec{a}=\vec{\Phi}</math>{{qquad}}<math>\Longrightarrow</math>{{qquad}}<math>\left\{\begin{array}{rcl} m(\ddot{\rho}-\omega^{2}\rho) & = & 0 \\ 2m\omega\dot{\rho} & = & \Phi \end{array}\right.</math></center>

Prescindiendo de la masa (factor no nulo) en la primera de estas ecuaciones (proyección radial de la segunda ley de Newton), se obtiene la ecuación diferencial buscada 

<center><math> \ddot{\rho}-\omega^{2}\rho = 0
</math></center>

==Solución particular==
[[Archivo:particula-tubo.gif|right]]

Derivando dos veces

<center><math>
\begin{array}{rcl}
\rho & = & A\mathrm{e}^{\omega t}\\
\dot{\rho} & = & \omega A\mathrm{e}^{\omega t}\\
\ddot{\rho} & = & \omega^2 A\mathrm{e}^{\omega t}
\end{array}
</math></center>

y sustituyendo

<center><math>\ddot{\rho}-\omega^2\rho = \omega^2 A\mathrm{e}^{\omega t} - \omega^2 A\mathrm{e}^{\omega t} = 0
</math></center>

luego efectivamente es una solución particular de esta ecuación.

===Fuerza de reacción===
Obtenemos la fuerza de reacción vincular a partir de la ecuación que se obtuvo al proyectar la segunda ley de Newton en dirección acimutal

<center><math>\vec{\Phi}=\Phi\,\vec{u}_{\theta} = 2m\omega\dot{\rho}\,\vec{u}_{\theta} = 2m\omega^2A\mathrm{e}^{\omega t}\,\vec{u}_{\theta}
</math></center>

===Potencia desarrollada===
La potencia desarrollada por la fuerza de reacción vincular es el producto de ésta por la
velocidad de la partícula. Esta velocidad vale

<center><math>\vec{v} = \dot{\rho}\,\vec{u}_{\rho}+\omega\rho\,\vec{u}_{\theta} = 
\omega A\mathrm{e}^{\omega t}(\vec{u}_{\rho}+\vec{u}_{\theta})</math></center>

Obsérvese que esta velocidad no es puramente radial, sino que forma un ángulo de 45&deg; con el tubo.

Hallamos el producto escalar y resulta

<center><math>P=\vec{\Phi}\cdot\vec{v}= 2m\omega^3 A^2\mathrm{e}^{2\omega t}</math></center>

===Trabajo en un intervalo===
El trabajo realizado en un intervalo es

<center><math>W = \int_{t_1}^{t_2}P\,\mathrm{d}t = \int_{t_1}^{t_2} 2m\omega^3 A^2\mathrm{e}^{2\omega
t}\mathrm{d}t = m\omega^2 A^2\left(\mathrm{e}^{2\omega t_2}-\mathrm{e}^{2\omega
t_1}\right)</math></center>

En este caso, <math>t_1</math> y <math>t_2</math> son tales que

<center><math>b = A \mathrm{e}^{\omega t_1}\qquad 2b = A \mathrm{e}^{\omega t_2}</math></center>

por lo que el trabajo en este desplazamiento vale
<center><math>
W = m\omega^2((2b)^2-b^2) = 3m\omega^2b^2\,</math></center>

==Energía cinética==
La energía cinética en cada instante para esta solución es, sustituyendo,

<center><math>K = \frac{1}{2}m|\vec{v}|^2 = m\left(\dot{\rho}^2+\rho^2\dot{\theta}^2\right)</math></center>

Dado que

<center><math>\dot{\rho}=A\omega\mathrm{e}^{\omega t}=\omega\rho\qquad\qquad\dot{\theta}=\omega</math></center>

esto se reduce a

<center><math>K=\frac{1}{2}m\left(\omega^2\rho^2+\rho^2\omega^2\right)=m\omega^2\rho^2</math></center>

Por tanto, el incremento entre <math>\rho=b\,</math> y <math>\rho=2b\,</math> vale
<center>
<math>\Delta K = m\omega^2(4b^2-b^2) = 3m\omega^2b^2\,</math></center>

que coincide con el trabajo antes obtenido, tal como predice el teorema de las fuerzas vivas.

Tenemos entonces que la partícula se va acelerando y ganando energía cinética. El trabajo necesario para este incremento lo realiza la fuerza de reacción vincular, que en este caso no es normal a la trayectoria, sino que forma un ángulo de 45º con ella. Para entender esto, cabe recordar que el tubo tiene un movimiento forzado (rota con velocidad angular constante) y por tanto supone para la partícula un vínculo liso pero reónomo.

[[Categoría:Problemas de Dinámica del Punto (GITI)]]