2.8. Ejemplo de movimiento helicoidal
De Laplace
Contenido |
1 Enunciado
Una partícula se mueve a lo largo de la hélice descrita por la ecuación paramétrica
![\vec{r}(\theta)=A\cos(\theta)\vec{\imath}+A\,\mathrm{sen}(\theta)\vec{\jmath}+\frac{b\theta}{2\pi}\vec{k}](/wiki/images/math/4/a/b/4abb860f31f21403ba0ca901d502cd17.png)
donde A y b son constantes conocidas. El movimiento de la partícula sigue la ley horaria
![\theta(t) = \Omega_0 t + \beta t^2\,](/wiki/images/math/b/8/4/b845cee57f7114c5e1253d37a99c0a69.png)
donde Ω0 y β son constantes conocidas.
- Determine el parámetro arco de la hélice descrita, como función del parámetro θ y como función del tiempo.
- Halle la rapidez del movimiento.
- Calcule la componente tangencial de la aceleración de la partícula en todo instante.
- Para el instante t = 0 calcule la velocidad y la aceleración de la partícula.
- Para el mismo instante, halle los vectores del triedro de Frenet, así como el radio de curvatura de la partícula y su aceleración normal.
2 Parámetro arco
Podemos determinar el parámetro arco empleando la variable θ según la relación
![\frac{\mathrm{d}s}{\mathrm{d}\theta}=\left|\frac{\mathrm{d}\vec{r}}{\mathrm{d}\theta}\right|](/wiki/images/math/b/1/0/b100bb87115878279a2c60a5f9978eaa.png)
Derivando y calculando el módulo
![\frac{\mathrm{d}\vec{r}}{\mathrm{d}\theta}=-A\,\mathrm{sen}(\theta)\vec{\imath}+A\cos(\theta)\vec{\jmath}+\frac{b}{2\pi}\vec{k}](/wiki/images/math/5/a/5/5a52f2001f0471bcda2e80ab796ff7ba.png)
El módulo de este vector vale
![\left|\frac{\mathrm{d}\vec{r}}{\mathrm{d}\theta}\right|=\sqrt{A^2+\frac{b^2}{4\pi^2}}](/wiki/images/math/8/d/3/8d3fa1ec41001798262cd349369ad11c.png)
Puesto que el módulo es independiente de θ, su integración es inmediata
![\frac{\mathrm{d}s}{\mathrm{d}\theta}=\sqrt{A^2+\frac{b^2}{4\pi^2}}](/wiki/images/math/1/c/5/1c5325fe443270d5a5276c619022bfc0.png)
![\Rightarrow](/wiki/images/math/d/f/0/df09aea884019cb88a2957126faba316.png)
![s=\theta\sqrt{A^2+\frac{b^2}{4\pi^2}}](/wiki/images/math/7/f/6/7f60fc4dd8bf4c6134f361de4a9063a0.png)
Y para obtener el parámetro arco en función del tiempo, basta sustituir en s(θ) la ley horaria θ(t) del enunciado
![s=(\Omega_0 t + \beta t^2)\sqrt{A^2+\frac{b^2}{4\pi^2}}](/wiki/images/math/4/0/8/40833071352d4b74f7f39c8618990917.png)
3 Celeridad
El cálculo de la rapidez es inmediato por derivación respecto al tiempo del parámetro arco
![v = \dot{s} = \left(\Omega_0+2\beta t\right)\sqrt{A^2+\frac{b^2}{4\pi^2}}](/wiki/images/math/b/e/c/bec18f6bae2793bf26550119e318bbca.png)
4 Aceleración tangencial
Obtenemos la componente tangencial de la aceleración derivando la celeridad respecto al tiempo
![a_t = \dot{v}= \ddot{s}= 2\beta \sqrt{A^2+\frac{b^2}{4\pi^2}}](/wiki/images/math/b/5/5/b5510d24dcdc881377c2652efbd10d26.png)
Puesto que esta aceleración tangencial es constante, el movimiento a lo largo de la hélice es uniformemente acelerado.
5 Velocidad y aceleración iniciales
Hallamos la ecuación horaria sustituyendo la ley horaria en la ecuación de la trayectoria
![\vec{r}(t)=A\cos(\Omega_0 t + \beta t^2)\vec{\imath}+A\,\mathrm{sen}(\Omega_0 t + \beta t^2)\vec{\jmath}+\frac{b(\Omega_0 t + \beta t^2)}{2\pi}\vec{k}](/wiki/images/math/9/5/f/95fb259e900a67c013c748a511f4151e.png)
Derivando en esta expresión respecto al tiempo
![\vec{v}(t)=\frac{\mathrm{d}\vec{r}}{\mathrm{d}t}=\left(\Omega_0+2\beta t\right)\left(-A\,\mathrm{sen}(\Omega_0 t + \beta t^2)\vec{\imath}+A\cos(\Omega_0 t + \beta t^2)\vec{\jmath}+\frac{b}{2\pi}\vec{k}\right)](/wiki/images/math/f/3/e/f3ec1751889dda058ddb974a8338780b.png)
Haciendo t = 0
![\vec{v}(0)=\Omega_0\left(A\vec{\jmath}+\frac{b}{2\pi}\vec{k}\right)](/wiki/images/math/f/c/5/fc51d2355a3201ed9175f26224829232.png)
Obtenemos la aceleración derivando la velocidad respecto al tiempo
![\vec{a}(t)=\frac{\mathrm{d}\vec{v}}{\mathrm{d}t}=2\beta\left(-A\,\mathrm{sen}(\Omega_0 t + \beta t^2)\vec{\imath}+A\cos(\Omega_0 t + \beta t^2)\vec{\jmath}+\frac{b}{2\pi}\vec{k}\right)](/wiki/images/math/4/c/7/4c7f761c85b3b72a89ba6a15dd4d4d25.png)
![+\left(\Omega_0+2\beta t\right)^2\left(-A\cos(\Omega_0 t + \beta t^2)\vec{\imath}-A\,\mathrm{sen}(\Omega_0 t + \beta t^2)\vec{\jmath}\right)](/wiki/images/math/7/4/c/74cd3ffc21525beb3011fb752f66ba6a.png)
Haciendo aquí t = 0
![\vec{a}(0)=-\Omega_0^2A\vec{\imath} + 2\beta A\vec{\jmath}+\frac{b\beta}{\pi}\vec{k}](/wiki/images/math/3/c/0/3c04bb3f8ac4c10cccb619c2aa042092.png)
6 Triedro de Frenet
Para el instante inicial hallamos el vector tangente a la circunferencia normalizando la velocidad
![\vec{T}(0)=\frac{\vec{v}(0)}{v(0)}=\frac{1}{\sqrt{A^2+{b^2}/{4\pi^2}}}\left(A\vec{\jmath}+\frac{b}{2\pi}\vec{k}\right)](/wiki/images/math/a/b/2/ab2471590349023952ae17ce1a841c9a.png)
Si lo escribimos en función de la inclinación de la hélice, α, definida en la teoría, queda:
![\vec{T}(0)=\cos(\alpha)\vec{\jmath}+\mathrm{sen}(\alpha)\vec{k}](/wiki/images/math/b/c/b/bcb256a06d9e8c140796339ee19328c7.png)
siendo evidente a partir de la figura que:
![\cos(\alpha)=\frac{A}{\sqrt{A^2+{b^2}/{4\pi^2}}}](/wiki/images/math/3/e/4/3e440060424e5675e8c046913937b11d.png)
![\mathrm{sen}(\alpha)=\frac{b/2\pi}{\sqrt{A^2+{b^2}/{4\pi^2}}}](/wiki/images/math/2/f/5/2f5cd767d27edcc2fc03456761329ac2.png)
El vector binormal lo hallamos normalizando el producto vectorial de la velocidad y la aceleración:
![\vec{v}(0)\times\vec{a}(0)=\left|\begin{matrix}\vec{\imath} & \vec{\jmath} & \vec{k} \\ 0 & \Omega_0A & \Omega_0b/2\pi \\-\Omega_0^2A &2\beta A & b\beta/\pi\end{matrix}\right|=A\Omega_0^3\left(-\frac{b}{2 \pi}\vec{\jmath}+A\vec{k}\right)](/wiki/images/math/c/7/0/c70edb4c9713aebfe28736e22f569f8a.png)
Resulta el vector
![\vec{B}(0) = \frac{1}{\sqrt{A^2+{b^2}/{4\pi^2}}}\left(-\frac{b}{2\pi}\vec{\jmath}+A\vec{k}\right)](/wiki/images/math/d/9/c/d9c2d5f36eaaa3f50512f86dba592ad3.png)
que es claramente ortogonal al vector tangente.
Escrito en función de α:
![\vec{B}(0)=-\mathrm{sen}(\alpha)\vec{\jmath}+\cos(\alpha)\vec{k}](/wiki/images/math/a/b/2/ab2ad836e74e06f19142869ca00ace84.png)
Multiplicando estos dos hallamos el vector normal
![\vec{N}(0)=\vec{B}(0)\times\vec{T}(0) = \frac{1}{A^2+b^2/4\pi^2}\left|\begin{matrix}\vec{\imath} & \vec{\jmath} & \vec{k} \\0 & -b\beta/\pi & A \\ 0 & A & b/2\pi \end{matrix}\right|=-\vec{\imath}](/wiki/images/math/8/8/0/880e0e1ad75942b2f061df532f5dcb71.png)
que es ortogonal a los dos anteriores.
De aquí tenemos que la aceleración normal en el instante inicial es igual a
![a_n(0) = \vec{a}(0)\cdot(-\vec{\imath})= A\Omega_0^2](/wiki/images/math/f/2/0/f209918a319ef46f38631fd0fed9812f.png)
y el radio de curvatura inicial vale
![R(0) = \frac{v^2(0)}{a_n(0)} = \frac{\Omega_0^2\left(A^2 + b^2/4\pi^2\right)}{A\Omega_0^2} = A + \frac{b^2}{4\pi^2A}](/wiki/images/math/3/5/6/356deb97b5806b6cdf37e5421399521c.png)
Puede demostrarse que este radio de curvatura es constante a lo largo de todo el movimiento.
Método alternativo: Tras calcular , podríamos haber calculado el vector aceleración normal como:
![\vec{a}_n(0) = \vec{a}(0)-a_t\vec{T}(0) = -A\Omega_0^2 \vec{\imath}](/wiki/images/math/6/d/f/6dfdad0bda02ee047e36eeeee273e309.png)
Tomando módulo y normalizando, obtendríamos la aceleración normal y el vector normal principal, respectivamente:
![a_n(0)=|\vec{a}_n(0)| = A\Omega_0^2 \,\,\,;\,\,\,\,\,\,\,\,\,\,\,\,\, \vec{N}(0)=\frac{\vec{a}_n(0)}{a_n(0)} = -\vec{\imath}](/wiki/images/math/9/4/e/94ee546c378a83c0188042a487c3429f.png)
Y el vector binormal se obtendría entonces como