Oscilador armónico tridimensional
De Laplace
Contenido |
1 Enunciado
Una partícula se mueve en tres dimensiones de forma tal que verifica la ecuación del oscilador armónico
![\vec{a}=-\omega^2\vec{r}](/wiki/images/math/f/2/a/f2a22c44c67d467642985afdcbae2235.png)
con . Su posición inicial es
.
- Para el caso
. ¿Qué tipo de movimiento describe la partícula?
- Para el caso
, ¿cómo es la trayectoria? ¿Qué tipo de movimiento describe la partícula?
- Suponga ahora que
, ¿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
,
y
.
2 Solución general
La solución general de la ecuación del oscilador armónico en 3D es de la forma
![\vec{r}(t) = \vec{r}_0\cos(\omega t)+\frac{\vec{v}_0}\mathrm{sen}(\omega t)](/wiki/images/math/0/2/5/025edd0c23c8cbfc1834834e4e090d36.png)
siendo la posición inicial y
la velocidad inicial de la partícula.
En este problema tenemos para todos los casos, empleando las unidades fundamentales del SI,
![\omega = 2\qquad\qquad\vec{r}_0=5\vec{\imath}\qquad\qquad\vec{v}_0=v_0\vec{\jmath}](/wiki/images/math/5/b/0/5b0c94d9730c7984a1fcfadfb4a254a0.png)
que al sustituir en la solución general nos dan la ecuación horaria
![\vec{r}=5\cos(2t)\vec{\imath}+\frac{v_0}{2}\mathrm{sen}(2t)\vec{\jmath}](/wiki/images/math/5/e/c/5ec85a418e756201bed75f85fb0f9478.png)
A partir de esta ecuación horaria pueden hallarse todas las características del movimiento.
3 Caso v0 = 0 m/s
El primer caso tiene velocidad inicial nula, lo que reduce la ecuación horaria a
![\vec{r}=5\cos(\omega t)\vec{\imath}](/wiki/images/math/0/2/3/023415ffa4cb42e940e6d113c2e8c228.png)
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.
4 Caso v0 = 10 m/s
En el segundo caso, la ecuación horario se reduce a
![\vec{r}=5\cos(\omega t)\vec{\imath}+5\,\mathrm{sen}(2 t)\vec{\jmath}](/wiki/images/math/d/a/5/da50c16c5cbf017aa7c8785ce43df651.png)
Este no es un movimiento rectilíneo. Sí es plano, porque tiene solo dos componentes, y verifica en todo momento que
![|\vec{r}|^2 = x^2 + y^2 = 5^2=\mathrm{cte}](/wiki/images/math/4/4/9/4496b8c18bd6ffc48f03655f582bed2e.png)
Por tanto se trata de un movimiento circular de radio 5 m alrededor del origen.
Además, podemos hallar la velocidad y rapidez del movimiento
![\vec{v}=-10º,\mathrm{sen}(2t)\vec{\imath}+10\cos(2t)\vec{\jmath}\qquad\Rightarrow\qquad |\vec{v}|=10\,\frac{\mathrm{m}}{\mathrm{s}}](/wiki/images/math/e/0/f/e0f9c473c3cf0a29a1a1ca50f6ddc783.png)
Al ser la rapidez constante el movimiento es circular uniforme.
5 Caso v0 = 8 m/s
El tercer caso tiene la ecuación horaria
![\vec{r}=5\cos(\omega t)\vec{\imath}+5\,\mathrm{sen}(2 t)\vec{\jmath}](/wiki/images/math/d/a/5/da50c16c5cbf017aa7c8785ce43df651.png)
Esta trayectoria no es ni rectilínea ni circular. Si separamos en componentes nos queda
![x = 5\cos(2t)\qquad\qquad y = 4\,\mathrm{sen}(2 t)\qquad\qquad z=0](/wiki/images/math/0/7/a/07ae70ac875fb9821a908cb9606394bd.png)
Podemos eliminar el tiempo de las dos primeras ecuaciones y combinarlas como
![\left(\frac{x}{5}\right)^2 + \left(\frac{y}{4}\right)^2 = 1\qquad\qquad z=0](/wiki/images/math/9/b/0/9b0d1c5e8e8873bb2f1f39f404e707b6.png)
que es la ecuación de una elipse de semiejes 5m y 4m
v0 = 0 | v0 = 10 | v0 = 8 |
---|---|---|
![]() | ![]() | ![]() |
6 Cálculo de magnitudes
6.1 Rapidez
Para cualquier movimiento, la velocidad la da la derivada del vector de posición respecto al tiempo
![\vec{v}=-10\,\mathrm{sen}(2t)\vec{\imath}+v_0\cos(2t)\vec{\jmath}](/wiki/images/math/e/2/6/e26ef200e4d6f1a7110cffc65950098b.png)
Podemos comprobar que, como corresponde, la velocidad en t=0 coincide con la velocidad inicial.
La rapidez es el módulo de esta velocidad
![|\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)}](/wiki/images/math/f/6/5/f651d8b3416f6d3a9932de38b309066a.png)
- En el caso particular v0 = 0 esta expresión se reduce a
![|\vec{v}|=10|\mathrm{sen}(2t)|](/wiki/images/math/9/a/b/9ab9ee1d1d427a156b43c2b9fed27557.png)
- que es una función oscilante, pero no sinusoidal, sino que va como el valor absoluto de un seno.
- Si v0 = 10, la rapidez vale
![|\vec{v}| = 10](/wiki/images/math/2/4/7/24790b5c63190457a876549a7699514e.png)
- que es constante, lo que indica que en este caso el movimiento es uniforme.
- En el caso v0 = 8 resulta una función oscilante
![|\vec{v}| = \sqrt{100\mathrm{sen}^2(2t)+ 64\cos^2(2t)}=\sqrt{100-36\cos^2(2 t)}](/wiki/images/math/5/6/e/56e55704e350b1c141d6962d4209b6e3.png)
- 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.
![Archivo:comparacion-rapidez.png](/wiki/images/8/89/Comparacion-rapidez.png)
En concreto, para los tres instantes del enunciado tenemos la siguiente tabla de valores
t=0 | t=π/4 | t=π/8 | |
---|---|---|---|
v0 = 0 | 0.00 | 7.07 | 10.0 |
v0 = 10 | 10.0 | 10.0 | 10.0 |
v0 = 8 | 8.00 | 9.06 | 10.0 |
6.2 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
![\vec{v}=-10\,\mathrm{sen}(2t)\vec{\imath}+v_0\cos(2t)\vec{\jmath}](/wiki/images/math/e/2/6/e26ef200e4d6f1a7110cffc65950098b.png)
y para las aceleraciones
![\vec{a}=-20\,\mathrm{cos}(2t)\vec{\imath}-2v_0\mathrm{sen}(2t)\vec{\jmath}](/wiki/images/math/a/4/c/a4c916349a061f0655abc46f8a77361d.png)
- En el caso v0 = 0 queda
t=0 | t=π/4 | t=π/8 | |
---|---|---|---|
v(m/s) | ![]() | ![]() | ![]() |
a(m/s²) | ![]() | ![]() | ![]() |
- En el caso v0 = 10
t=0 | t=π/4 | t=π/8 | |
---|---|---|---|
v(m/s) | ![]() | ![]() | ![]() |
a(m/s²) | ![]() | ![]() | ![]() |
- En el caso v0 = 8
t=0 | t=π/4 | t=π/8 | |
---|---|---|---|
v(m/s) | ![]() | ![]() | ![]() |
a(m/s²) | ![]() | ![]() | ![]() |
6.2.1 Aceleración tangencial
La aceleración tangencial la calculamos como la proyección de sobre
![a_t=\frac{\vec{a}\cdot\vec{v}}{|\vec{v}|}](/wiki/images/math/6/1/b/61b26b1dc9ed1483ef7996b54d03c572.png)
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
t=0 | t=π/4 | t=π/8 | |
---|---|---|---|
v0 = 0 | ??? | 14.1 | 0.0 |
v0 = 10 | 0.0 | 0.0 | 0.0 |
v0 = 8 | 0.0 | 3.98 | 0.0 |
En el caso v0 = 0 la velocidad es nula en t=0 y por ello la aceleración tangencial no está definida en ese instante.
6.2.2 Aceleración normal
El valor de la aceleración normal escalar lo podemos hallar de varias formas
![a_n=\frac{|\vec{v}\times\vec{a}|}{|\vec{v}|}=\sqrt{|\vec{a}|^2-a_t^2}](/wiki/images/math/9/7/d/97d8c9ab2e13e586d6740060ae0c0982.png)
De nuevo se anula en varios de los casos por ser la aceleración paralela a la velocidad
t=0 | t=π/4 | t=π/8 | |
---|---|---|---|
v0 = 0 | ??? | 0.0 | 0.0 |
v0 = 10 | 20 | 20 | 20 |
v0 = 8 | 20 | 17.7 | 16 |
6.3 Vectores tangente y normal
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.
t=0 | t=π/4 | t=π/8 | |
---|---|---|---|
v = 0 | ![]() | ![]() | ![]() |
v0 = 10 | ![]() | ![]() | ![]() |
v0 = 8 | ![]() | ![]() | ![]() |