Ver el código de «Movimiento bajo fuerza central en polares»

==Enunciado==
Sea una partícula <math>P</math> de masa <math>m</math> cuyo movimiento en el plano OXY se describe mediante coordenadas polares.

# Deduzca la expresión de su velocidad areolar respecto al origen de coordenadas <math>O</math>, y compruebe que la misma es constante en el tiempo si el movimiento transcurre bajo la acción de una fuerza central en <math>O</math>.
# Sabiendo que la partícula recorre la espiral <math>\rho=\rho_{0}e^{\theta}\,</math> sometida a una fuerza central en <math>O</math> y con condiciones iniciales <math>\theta(0)=0\,</math> y <math>\dot{\theta}(0)=\omega_0</math>, determine las ecuaciones horarias <math>\rho = \rho(t)\,</math> y <math>\theta = \theta(t)\,</math>.
<!--
, así como el valor de dicha fuerza en función de la coordenada radial <math>\vec{F} = \vec{F}(\rho)</math>.
-->

==Expresión general de la velocidad areolar==

Sustituyendo las expresiones generales (en coordenadas polares) de los vectores de posición y velocidad de una partícula

<center><math>\begin{array}{rcl} \vec{r} & = & \rho\,\vec{u}_{\rho} \\ \vec{v} & = & \dot{\rho}\,\vec{u}_{\rho}+\rho\,\dot{\theta}\,\vec{u}_{\theta}\end{array}</math></center>

en la definición de la velocidad areolar de la partícula respecto al origen de coordenadas <math>O\,</math>, se obtiene

<center><math>\vec{V}_A=\frac{1}{2}\,\vec{r}\times\vec{v}=\frac{1}{2}\,\rho\,\vec{u}_{\rho}\times(\dot{\rho}\,\vec{u}_{\rho}+\rho\,\dot{\theta}\,\vec{u}_{\theta})=\frac{1}{2}\,\rho^2\dot{\theta}\,\vec{k}</math></center>

==Conservación de la velocidad areolar==

Teniendo en cuenta que la dirección del vector velocidad areolar es constante (dirección perpendicular al plano del movimiento), es evidente que la citada magnitud vectorial se conservará constante en el tiempo si y sólo si es constante la velocidad areolar escalar <math>V_A</math>

<center><math>\vec{V}_A=V_A\,\vec{k}=\frac{1}{2}\,\rho^2\dot{\theta}\,\vec{k}=\overrightarrow{\mathrm{cte}}\,\,\,\,\,\,\Longleftrightarrow\,\,\,\,\,\,V_A=\frac{1}{2}\,\rho^2\dot{\theta}=\mathrm{cte} </math></center>

Ahora bien, derivando respecto al tiempo <math>V_{\! A}</math> se obtiene

<center><math>\frac{\mathrm{d}V_A}{\mathrm{d}t}=\frac{1}{2}\,\frac{\mathrm{d}(\rho^2\dot{\theta})}{\mathrm{d}t}=\frac{1}{2}\,(2\rho\dot{\rho}\,\dot{\theta}+\rho^2\ddot{\theta})=\frac{\rho}{2}\,(2\dot{\rho}\,\dot{\theta}+\rho\,\ddot{\theta})=\frac{\rho}{2}\,a_{\theta} </math></center>

Comprobamos, por tanto, que la conservación de la velocidad areolar como una magnitud constante a lo largo del tiempo se produce si y sólo si la componente acimutal de la aceleración es permanentemente nula, lo cual conlleva de forma obvia la centralidad de la fuerza respecto al origen de coordenadas:

<center><math>V_A=\mathrm{cte}\,\,\,\,\,\,\Longleftrightarrow\,\,\,\,\,\,\frac{\mathrm{d}V_A}{\mathrm{d}t}=0\,\,\,\,\,\,\Longleftrightarrow\,\,\,\,\,\,a_{\theta}=0\,\,\,\,\,\,\Longleftrightarrow\,\,\,\,\,\,\vec{F}=m\vec{a}=ma_{\rho}\,\vec{u}_{\rho}=F_{\rho}\,\vec{u}_{\rho}</math></center>

Obsérvese que la conservación de la velocidad areolar está intímamente ligada a la conservación del momento cinético (respecto al mismo punto <math>O\,</math>), ya que

<center><math>\vec{L}_O=\vec{r}\times m\vec{v}=2m\vec{V}_A</math></center>

siendo la propiedad de paralelismo <math>\vec{F}\parallel\vec{r}</math>, propia de las fuerzas centrales en <math>O\,</math>, la que provoca que tanto el momento cinético como la velocidad areolar respecto a <math>O\,</math> sean constantes.

==Ecuaciones horarias==
Como nos dicen que la partícula recorre la espiral bajo la acción de una fuerza central en <math>O\,</math>, tenemos garantizado que la velocidad areolar escalar es una constante, pudiéndose evaluar el valor de dicha constante a partir de las condiciones iniciales dadas

<center><math>\left.\begin{array}{l}\dot{\theta}(0)=\omega_0 \\ \rho(0)=\rho_0\,e^{\theta(0)}=\rho_0\,e^{0}=\rho_0 \end{array}\right|\,\,\,\,\,\,\longrightarrow\,\,\,\,\,\,   V_A=\frac{1}{2}\,\rho^2\dot{\theta}=\mathrm{cte}=\frac{1}{2}\,[\rho(0)]^2\dot{\theta}(0)=\frac{1}{2}\,\rho_0^2\omega_0</math></center>

Prescindiendo del factor <math>1/2</math> en ambos miembros de la ecuación anterior, expresando <math>\rho\,</math> como función de <math>\theta\,</math> (ecuación de la trayectoria), separando las variables <math>\theta\,</math> y <math>t</math>, y finalmente integrando entre el instante inicial y un instante genérico, se obtiene:

<center><math>\rho^2\dot{\theta}=\rho_0^2\omega_0\,\,\,\,\rightarrow\,\,\,\,\rho_0^2e^{2\theta}\frac{\mathrm{d}\theta}{\mathrm{d}t}=\rho_0^2\omega_0\,\,\,\,\rightarrow\,\,\,\,e^{2\theta}\mathrm{d}\theta=\omega_0\mathrm{d}t\,\,\,\,\rightarrow\,\,\,\,\int_{0}^{\theta}e^{2\theta}\mathrm{d}\theta=\omega_0\int_{0}^{t}\mathrm{d}t\,\,\,\,\rightarrow\,\,\,\,\frac{e^{2\theta}-1}{2}=\omega_0t</math></center>

de donde, despejando <math>\theta\,</math>, se llega a la ecuación horaria acimutal:

<center><math>\theta(t)=\frac{1}{2}\,\mathrm{ln}(1+2\omega_0t)</math></center>

y, sustituyendo ésta en la ecuación de la espiral (trayectoria), se obtiene la ecuación horaria radial:

<center><math>\rho(t)=\rho_0e^{\mathrm{ln}(\sqrt{1+2\omega_0t}\,)}=\rho_0\sqrt{1+2\omega_0t}</math></center>
<!--
==Fuerza en función de la coordenada radial==
Sabemos que la expresión del vector aceleración de la partícula en coordenadas polares viene dada por

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

Pero, en el movimiento que nos ocupa, la fuerza que actúa sobre la partícula es central en <math>O\,</math>. En consecuencia, tanto la fuerza como la aceleración sólo tienen componente radial:

<center><math>\vec{F}=F_{\rho}\,\vec{u}_{\rho}=m\,a_{\rho}\,\vec{u}_{\rho}=m\,(\ddot{\rho}-\rho\,\dot{\theta}^2)\,\vec{u}_{\rho}</math></center>

Como nos piden que expresemos la fuerza en función de la coordenada radial <math>\rho\,</math>, conviene calcular <math>\ddot{\rho}\,</math> y <math>\dot{\theta}\,</math> en función de <math>\rho\,</math>. Es inmediato expresar <math>\dot{\theta}\,</math> en función de <math>\rho\,</math> a partir de la relación que obtuvimos en el apartado anterior al evaluar la velocidad areolar constante:

<center><math>\rho^2\dot{\theta}=\rho_0^2\omega_0\,\,\,\,\,\,\longrightarrow\,\,\,\,\,\,\dot{\theta}=\frac{\rho_0^2\omega_0}{\rho^2}</math></center>

Y en cuanto a <math>\ddot{\rho}\,</math>, hemos de calcular antes <math>\dot{\rho}\,</math>. Aplicando la regla de la cadena:

<center><math>\dot{\rho}=\frac{\mathrm{d}\rho}{\mathrm{d}\theta}\dot{\theta}=\rho\frac{\rho_0^2\omega_0}{\rho^2}=\frac{\rho_0^2\omega_0}{\rho}</math></center>

donde se ha tenido en cuenta que

<center><math>\rho=\rho_0e^{\theta}\,\,\,\,\,\,\longrightarrow\,\,\,\,\,\,\frac{\mathrm{d}\rho}{\mathrm{d}\theta}=\rho_0e^{\theta}=\rho</math></center>

Y, aplicando de nuevo la regla de la cadena, calculamos ya <math>\ddot{\rho}\,</math>

<center><math>\ddot{\rho}=\frac{\mathrm{d}\dot{\rho}}{\mathrm{d}\rho}\dot{\rho}=-\frac{\rho_0^2\omega_0}{\rho^2}\frac{\rho_0^2\omega_0}{\rho}=-\frac{\rho_0^4\omega_0^2}{\rho^3}</math></center>

Finalmente, sustituyendo en la expresión de la fuerza, se obtiene:

<center><math>\vec{F}=m\left(-\displaystyle\frac{\rho_0^4\omega_0^2}{\rho^3}-\rho\frac{\rho_0^4\omega_0^2}{\rho^4}\right)\,\vec{u}_{\rho}=-\frac{2m\rho_0^4\omega_0^2}{\rho^3}\,\vec{u}_{\rho}</math></center>
-->
[[Categoría:Problemas de Dinámica del Punto (GITI)]]