Ver el código de «Oscilador armónico tridimensional»

==Enunciado==
Una partícula se mueve en tres dimensiones de forma tal que verifica la ecuación del oscilador armónico

<center><math>\vec{a}=-\omega^2\vec{r}</math></center>

con <math>\omega = 2.0\,\mathrm{rad}/\mathrm{s}</math>. Su posición inicial es <math>\vec{r}_0=5\,\vec{\imath}\ (\mathrm{m})</math>.

# Para el caso <math>\vec{v}_0=\vec{0}\,\mathrm{m}/\mathrm{s}</math>. ¿Qué tipo de movimiento describe la partícula?
# Para el caso <math>\vec{v}_0=10.0\,\vec{\jmath}\,\mathrm{m}/\mathrm{s}</math>, ¿cómo es la trayectoria? ¿Qué tipo de movimiento describe la partícula?
# Suponga ahora que <math>\vec{v}_0=8.0\,\vec{\jmath}\,\mathrm{m}/\mathrm{s}</math>, ¿cómo es ahora la trayectoria de la partícula?
# Para los tres casos anteriores, determine
:# la rapidez,
:# las componentes intrínsecas de la aceleración,
:# los vectores tangente y normal,
:# el radio de curvatura y el centro de curvatura.
::para los instantes <math>t=0\,</math>, <math>t=0.25\pi\,\mathrm{s}</math> y <math>t = 0.125\pi\,\mathrm{s}</math>.

==Solución general==
La solución general de la ecuación del oscilador armónico en 3D es de la forma

<center><math>\vec{r}(t) = \vec{r}_0\cos(\omega t)+\frac{\vec{v}_0}{\omega}\mathrm{sen}(\omega t)</math></center>

siendo <math>\vec{r}_0</math> la posición inicial y <math>\vec{v}_0</math> la velocidad inicial de la partícula.

En este problema tenemos para todos los casos, empleando las unidades fundamentales del SI,

<center><math>\omega = 2\qquad\qquad\vec{r}_0=5\vec{\imath}\qquad\qquad\vec{v}_0=v_0\vec{\jmath}</math></center>

que al sustituir en la solución general nos dan la ecuación horaria

<center><math>\vec{r}=5\cos(2t)\vec{\imath}+\frac{v_0}{2}\mathrm{sen}(2t)\vec{\jmath}</math></center>

A partir de esta ecuación horaria pueden hallarse todas las características del movimiento.

==Caso v<sub>0</sub> = 0 m/s==
El primer caso tiene velocidad inicial nula, lo que reduce la ecuación horaria a

<center><math>\vec{r}=5\cos(2 t)\vec{\imath}</math></center>

Esta es la ecuación de movimiento rectilíneo a lo largo del eje OX. Por cumplir la ecuación del oscilador armónico y ser rectilíneo, el movimiento es armónico simple.

==Caso v<sub>0</sub> = 10 m/s==
En el segundo caso, la ecuación horario se reduce a

<center><math>\vec{r}=5\cos(2 t)\vec{\imath}+5\,\mathrm{sen}(2 t)\vec{\jmath}</math></center>

Este no es un movimiento rectilíneo. Sí es plano, porque tiene solo dos componentes, y verifica en todo momento que

<center><math>|\vec{r}|^2 = x^2 + y^2 = 5^2=\mathrm{cte}</math></center>

Por tanto se trata de un movimiento circular de radio 5&thinsp;m alrededor del origen.

Además, podemos hallar la velocidad y rapidez del movimiento

<center><math>\vec{v}=-10\,\mathrm{sen}(2t)\vec{\imath}+10\cos(2t)\vec{\jmath}\qquad\Rightarrow\qquad |\vec{v}|=10\,\frac{\mathrm{m}}{\mathrm{s}}</math></center>

Al ser la rapidez constante el movimiento es circular uniforme.

==Caso v<sub>0</sub> = 8 m/s==
El tercer caso tiene la ecuación horaria

<center><math>\vec{r}=5\cos(2 t)\vec{\imath}+4\,\mathrm{sen}(2 t)\vec{\jmath}</math></center>

Esta trayectoria no es ni rectilínea ni circular. Si separamos en componentes nos queda

<center><math>x = 5\cos(2t)\qquad\qquad y = 4\,\mathrm{sen}(2 t)\qquad\qquad z=0</math></center>

Podemos eliminar el tiempo de las dos primeras ecuaciones y combinarlas como

<center><math>\left(\frac{x}{5}\right)^2 + \left(\frac{y}{4}\right)^2 = 1\qquad\qquad z=0</math></center>

que es la ecuación de una elipse de semiejes 5m y 4m

{| class="bordeado"
|-
! <math>v_0=0</math>
! <math>v_0=10</math>
! <math>v_0=8</math>
|-
| [[Archivo:OC3D-0.gif]]
| [[Archivo:OC3D-5.gif]]
| [[Archivo:OC3D-4.gif]]
|}

==Cálculo de magnitudes==
===Rapidez===
Para cualquier movimiento, la velocidad la da la derivada del vector de posición respecto al tiempo

<center><math>\vec{v}=-10\,\mathrm{sen}(2t)\vec{\imath}+v_0\cos(2t)\vec{\jmath}</math></center>

Podemos comprobar que, como corresponde, la velocidad en t=0 coincide con la velocidad inicial.

La rapidez es el módulo de esta velocidad

<center><math>|\vec{v}| = \sqrt{10^2\mathrm{sen}^2(2t)+ v_0^2\cos^2(2t)}=\sqrt{100-(100-v_0^2)\cos^2(2 t)}</math></center>

* En el caso particular <math>v_0=0</math> esta expresión se reduce a

<center><math>|\vec{v}|=10|\mathrm{sen}(2t)|</math></center>

:que es una función oscilante, pero no sinusoidal, sino que va como el valor absoluto de un seno.

* Si <math>v_0=10</math>, la rapidez vale

<center><math>|\vec{v}| = 10</math></center>

:que es constante, lo que indica que en este caso el movimiento es uniforme.

* En el caso <math>v_0=8</math> resulta una función oscilante

<center><math>|\vec{v}| = \sqrt{100\mathrm{sen}^2(2t)+ 64\cos^2(2t)}=\sqrt{100-36\cos^2(2 t)}</math></center>

:Esta función varía desde 8m/s en el instante inicial (que es su rapidez mínima) a 10/m/s cuando el coseno se anula.

<center>[[Archivo:comparacion-rapidez.png]]</center>

En concreto, para los tres instantes del enunciado tenemos la siguiente tabla de valores

{| class="bordeado"
|-
! <math>|\vec{v}|</math>
! t=0
! t=&pi;/4
! t=&pi;/8
|-
| <math>v_0=0</math>
| 0.00
| 7.07 
| 10.0
|-
| <math>v_0=10</math>
| 10.0
| 10.0 
| 10.0
|-
| <math>v_0=8</math>
| 8.00
| 9.06 
| 10.0
|}

===Componentes de la aceleración===
Para hallar las componentes intrínsecas, en lugar de operar con expresiones dependientes del tiempo, vamos a calcular en primer lugar las velocidades y aceleraciones para los casos indicados y posteriormente calcularemos las componentes a partir de las proyecciones de los vectores.

Para la velocidad tenemos

<center><math>\vec{v}=-10\,\mathrm{sen}(2t)\vec{\imath}+v_0\cos(2t)\vec{\jmath}</math></center>

y para las aceleraciones

<center><math>\vec{a}=-20\,\mathrm{cos}(2t)\vec{\imath}-2v_0\mathrm{sen}(2t)\vec{\jmath}</math></center>

* En el caso <math>v_0=0</math> queda

{| class="bordeado"
|-
!
! t=0
! t=&pi;/4
! t=&pi;/8
|-
| v(m/s)
| <math>\vec{0}</math>
| <math>-7.07\vec{\imath}</math>
| <math>-10\vec{\imath}</math>
|-
| a(m/s&sup2;)
| <math>-20\vec{\imath}</math>
| <math>-14.1\vec{\imath}</math>
| <math>\vec{0}</math>
|}

* En el caso <math>v_0=10</math>

{| class="bordeado"
|-
!
! t=0
! t=&pi;/4
! t=&pi;/8
|-
| v(m/s)
| <math>10\vec{\jmath}</math>
| <math>-7.07\vec{\imath}+7.07\vec{\jmath}</math>
| <math>-10\vec{\imath}</math>
|-
| a(m/s&sup2;)
| <math>-20\vec{\imath}</math>
| <math>-14.1\vec{\imath}-14.1\vec{\jmath}</math>
| <math>-20\vec{\jmath}</math>
|}

* En el caso <math>v_0=8</math>

{| class="bordeado"
|-
!
! t=0
! t=&pi;/4
! t=&pi;/8
|-
| v(m/s)
| <math>8\vec{\jmath}</math>
| <math>-7.07\vec{\imath}+5.66\vec{\jmath}</math>
| <math>-10\vec{\imath}</math>
|-
| a(m/s&sup2;)
| <math>-20\vec{\imath}</math>
| <math>-14.1\vec{\imath}+11.3\vec{\jmath}</math>
| <math>-16\vec{\imath}</math>
|}

====Aceleración tangencial====
La aceleración tangencial la calculamos como la proyección de <math>\vec{a}</math> sobre <math>\vec{v}</math>

<center><math>a_t=\frac{\vec{a}\cdot\vec{v}}{|\vec{v}|}</math></center>

De las tablas anteriores podemos ver que en muchos de los casos solicitados esta cantidad es nula, por ser la aceleración ortogonal a la velocidad

{| class="bordeado"
|-
! <math>a_t</math>
! t=0
! t=&pi;/4
! t=&pi;/8
|-
| <math>v_0=0</math>
| ???
| 14.1
| 0.0
|-
| <math>v_0=10</math>
| 0.0
| 0.0
| 0.0
|-
| <math>v_0=8</math>
| 0.0
| 3.98
| 0.0
|}

En el caso <math>v_0=0</math> la velocidad es nula en t=0 y por ello la aceleración tangencial no está definida en ese instante.

====Aceleración normal====
El valor de la aceleración normal escalar lo podemos hallar de varias formas

<center><math>a_n=\frac{|\vec{v}\times\vec{a}|}{|\vec{v}|}=\sqrt{|\vec{a}|^2-a_t^2}</math></center>

De nuevo se anula en varios de los casos por ser la aceleración paralela a la velocidad

{| class="bordeado"
|-
!<math>a_n</math>
! t=0
! t=&pi;/4
! t=&pi;/8
|-
| <math>v_0=0</math>
| ???
| 0.0
| 0.0
|-
| <math>v_0=10</math>
| 20
| 20
| 20
|-
| <math>v_0=8</math>
| 20
| 17.7
| 16
|}

===Vectores tangente y normal===
====Vector tangente====
El vector tangente es el unitario en la dirección de la velocidad, por lo que simplemente debemos dividir cada una de las velocidades anteriores por la rapidez correspondiente.

{| class="bordeado"
|-
! <math>\vec{T}</math>
! t=0
! t=&pi;/4
! t=&pi;/8
|-
| <math>v_0=0</math>
| ???
| <math>-\vec{\imath}</math>
| <math>-\vec{\imath}</math>
|-
| <math>v_0=10</math>
| <math>\vec{\jmath}</math>
| <math>-0.707\vec{\imath}+0.707\vec{\jmath}</math>
| <math>-\vec{\imath}</math>
|-
| <math>v_0=8</math>
| <math>\vec{\jmath}</math>
| <math>-0.780\vec{\imath}+0.625\vec{\jmath}</math>
| <math>-\vec{\imath}</math>
|}

====Vector normal====
El vector normal es el unitario en la dirección y sentido de la aceleración normal, por lo que solo hay que normalizar la tabla correspondiente

{| class="bordeado"
|-
! <math>\vec{N}</math>
! t=0
! t=&pi;/4
! t=&pi;/8
|-
| <math>v_0=0</math>
| ???
| ???
| ???
|-
| <math>v_0=10</math>
| <math>-\vec{\imath}</math>
| <math>-0.707\vec{\imath}-0.707\vec{\jmath}</math>
| <math>-\vec{\jmath}</math>
|-
| <math>v_0=8</math>
| <math>-\vec{\imath}</math>
| <math>-0.625\vec{\imath}-0.780\vec{\jmath}</math>
| <math>-\vec{\jmath}</math>
|}

===Radio y centro de curvatura===
====Radio de curvatura====
El radio de curvatura puede calcularse como

<center><math>R = \frac{|\vec{v}|^2}{a_n}=\frac{|\vec{v}|^3}{|\vec{v}\times\vec{a}|}</math></center>

Lo que da la tabla

{| class="bordeado"
|-
! R
! t=0
! t=&pi;/4
! t=&pi;/8
|-
| <math>v_0=0</math>
| ???
| <math>\infty\,</math>
| <math>\infty</math>
|-
| <math>v_0=10</math>
| 5
| 5
| 5
|-
| <math>v_0=8</math>
| 3.2
| 4.64
| 6.25
|}

En el primer caso el radio de curvatura no está definido porque se trata de un punto en el que la partícula se da la vuelta en su movimiento. En el segundo caso resulta un valor constante como corresponde a un movimiento circular.

====Centro de curvatura====
Una vez que tenemos el radio de curvatura y el vector normal localizamos el centro de curvatura como

<center><math>\vec{r}_c = \vec{r}+R\vec{N}</math></center>

y resulta la tabla

{| class="bordeado"
|-
! <math>\vec{r}_c</math>
! t=0
! t=&pi;/4
! t=&pi;/8
|-
| <math>v_0=0</math>
| ???
| ???
| ???
|-
| <math>v_0=10</math>
| <math>\vec{0}</math>
| <math>\vec{0}</math>
| <math>\vec{0}</math>
|-
|-
| <math>v_0=8</math>
| <math>1.8\vec{\imath}</math>
| <math>0.636\vec{\imath}-0.795\vec{\jmath}</math>
| <math>-2.25\vec{\jmath}</math>
|}

Podemos ilustrar las diferentes situaciones para el caso <math>v_0=8</math>, que es el menos trivial

{| class="bordeado"
|-
! t=0
! t = &pi;/8
! t=&pi;4
|-
| [[Archivo:OA3D-8-01.png]]
| [[Archivo:OA3D-8-02.png]]
| [[Archivo:OA3D-8-03.png]]
|}

[[Categoría:Problemas de cinemática tridimensional (GIOI)]]