Articulo de referencia

Tensor de velocidad angular

El tensor de velocidad angular es una matriz antisimétrica definida por: Ω = ( 0 − ω z ω y ω z 0 − ω incógnita − ω y ω incógnita 0 ) {\displaystyle \Omega ={\begin{pmatrix}0&-\o...

El tensor de velocidad angular es una matriz antisimétrica definida por:

Ω=(0ωzωyωz0ωincógnitaωyωincógnita0){\displaystyle \Omega ={\begin{pmatrix}0&-\omega _{z}&\omega _{y}\\\omega _{z}&0&-\omega _{x}\\-\omega _{y}&\omega _{x}&0\\\end{pmatrix}}}

Los elementos escalares anteriores corresponden a las componentes del vector de velocidad angular.ω=(ωincógnita,ωy,ωz){\displaystyle {\boldsymbol {\omega }}=(\omega _ {x}, \ omega _ {y}, \ omega _ {z})}.

Esta es una matriz de rotación infinitesimal . La transformación lineal Ω actúa como un producto vectorial.(ω×){\displaystyle ({\boldsymbol {\omega }}\times )}:

ω×r=Ωr{\displaystyle {\boldsymbol {\omega }}\times {\boldsymbol {r}}=\Omega {\boldsymbol {r}}}

dónder{\displaystyle {\boldsymbol {r}}}es un vector de posición .

Cuando se multiplica por una diferencia de tiempo, da como resultado el tensor de desplazamiento angular .

Cálculo del tensor de velocidad angular de un sistema de referencia giratorio.

Un vectorr{\displaystyle \mathbf {r} }que experimenta un movimiento circular uniforme alrededor de un eje fijo satisface:

drdt=ω×r=Ωr{\displaystyle {\frac {d\mathbf {r} }{dt}}={\boldsymbol {\omega }}\times \mathbf {r} =\Omega \mathbf {r} }

DejarA(t)=[mi1(t) mi2(t) mi3(t)]{\displaystyle A(t)=[\mathbf {e} _{1}(t)\ \mathbf {e} _{2}(t)\ \mathbf {e} _{3}(t)]} sea ​​la matriz de orientación de un marco, cuyas columnasmi1{\displaystyle \mathbf {e} _{1}},mi2{\displaystyle \mathbf {e} _{2}}, ymi3{\displaystyle \mathbf {e} _{3}}son los vectores de coordenadas ortonormales móviles del sistema de referencia. Podemos obtener el tensor de velocidad angular Ω ( t ) de A ( t ) de la siguiente manera:

La velocidad angularω{\displaystyle \omega }debe ser el mismo para cada uno de los vectores columna.mii{\displaystyle \mathbf {e} _ {i}}, entonces tenemos:

dAdt=[dmi1dtdmi2dtdmi3dt]=[ω×mi1ω×mi2ω×mi3]=[Ωmi1Ωmi2Ωmi3]=ΩA,{\displaystyle {\begin{aligned}{\frac {dA}{dt}}&={\begin{bmatrix}{\dfrac {d\mathbf {e} _{1}}{dt}}&{\dfrac {d\mathbf {e} _{2}}{dt}}&{\dfrac {d\mathbf {e} _{3}}{dt}}\end{bmatrix}}\\&={\begin{bmatrix}\omega \times \mathbf {e} _{1}&\omega \times \mathbf {e} _{2}&\omega \times \mathbf {e} _{3}\end{bmatrix}}\\&={\begin{bmatrix}\Omega \mathbf {e} _{1}&\Omega \mathbf {e} _{2}&\Omega \mathbf {e} _{3}\end{bmatrix}}\\&=\Omega A,\end{aligned}}}

lo cual se cumple incluso si A ( t ) no gira uniformemente. Por lo tanto, el tensor de velocidad angular es:

Ω=dAdtA1=dAdtAT,{\displaystyle \Omega ={\frac {dA}{dt}}A^{-1}={\frac {dA}{dt}}A^{\mathsf {T}},}

ya que la inversa de una matriz ortogonalA{\displaystyle A}es su transpuestaAT{\displaystyle A^{\mathsf {T}}}.

Propiedades

En general, la velocidad angular en un espacio n -dimensional es la derivada temporal del tensor de desplazamiento angular, que es un tensor antisimétrico de segundo rango .

Este tensor Ω tendrá n ( n −1)/2 componentes independientes, que es la dimensión del álgebra de Lie del grupo de Lie de rotaciones de un espacio de producto interno n -dimensional . [ 1 ]

Dualidad con respecto al vector velocidad

En tres dimensiones, la velocidad angular puede representarse mediante un pseudovector porque los tensores de segundo rango son duales a los pseudovectores en tres dimensiones. Dado que el tensor de velocidad angular Ω = Ω ( t ) es una matriz antisimétrica :

Ω=(0ωzωyωz0ωincógnitaωyωincógnita0),{\displaystyle \Omega ={\begin{pmatrix}0&-\omega _{z}&\omega _{y}\\\omega _{z}&0&-\omega _{x}\\-\omega _{y}&\omega _{x}&0\\\end{pmatrix}},}

Su dual de Hodge es un vector, que es precisamente el vector de velocidad angular anterior.ω=[ωincógnita,ωy,ωz]{\displaystyle {\boldsymbol {\omega }}=[\omega _{x},\omega _{y},\omega _{z}]}.

Exponencial de Ω

Si conocemos un sistema de referencia inicial A (0) y se nos da un tensor de velocidad angular constante Ω , podemos obtener A ( t ) para cualquier t dado . Recordemos la ecuación diferencial matricial:

dAdt=ΩA.{\displaystyle {\frac {dA}{dt}}=\Omega \cdot A.}

Esta ecuación se puede integrar para obtener:

A(t)=miWtA(0),{\displaystyle A(t)=e^{Wt}A(0),}

lo cual muestra una conexión con el grupo de rotaciones de Lie.

Ω es antisimétrico

Demostramos que el tensor de velocidad angular es antisimétrico , es decirΩ=dA(t)dtAT{\displaystyle \Omega ={\frac {dA(t)}{dt}}\cdot A^{\text{T}}}SatisfaceΩT=Ω{\displaystyle \Omega ^{\text{T}}=-\Omega }.

Una matriz de rotación A es ortogonal, inversa a su transpuesta, por lo que tenemosI=AAT{\displaystyle I=A\cdot A^{\text{T}}}. ParaA=A(t){\displaystyle A=A(t)}una matriz de marco, tomando la derivada temporal de la ecuación da como resultado:

0=dAdtAT+AdATdt{\displaystyle 0={\frac {dA}{dt}}A^{\text{T}}+A{\frac {dA^{\text{T}}}{dt}}}

Aplicando la fórmula(AB)T=BTAT{\displaystyle (AB)^{\text{T}}=B^{\text{T}}A^{\text{T}}},

0=dAdtAT+(dAdtAT)T=Ω+ΩT{\displaystyle 0={\frac {dA}{dt}}A^{\text{T}}+\left({\frac {dA}{dt}}A^{\text{T}}\right)^{\text{T}}=\Omega +\Omega ^{\text{T}}}

Por lo tanto, Ω es el negativo de su transpuesta, lo que implica que es antisimétrico.

Descripción sin coordenadas

En cualquier instantet{\displaystyle t}, el tensor de velocidad angular representa un mapeo lineal entre el vector de posiciónr(t){\displaystyle \mathbf {r} (t)}y los vectores de velocidadv(t){\displaystyle \mathbf {v} (t)}de un punto en un cuerpo rígido que gira alrededor del origen:

v=Ωr.{\displaystyle \mathbf {v} =\Omega \mathbf {r} .}

La relación entre este mapa lineal y el pseudovector de velocidad angularω{\displaystyle {\boldsymbol {\omega }}}es lo siguiente.

Debido a que Ω es la derivada de una transformación ortogonal , la forma bilineal

B(r,s)=(Ωr)s{\displaystyle B(\mathbf {r} ,\mathbf {s} )=(\Omega \mathbf {r} )\cdot \mathbf {s} }

es antisimétrica . Por lo tanto, podemos aplicar el hecho del álgebra exterior de que existe una única forma lineal.L{\displaystyle L}enΛ2V{\displaystyle \Lambda ^{2}V}eso

L(rs)=B(r,s){\displaystyle L(\mathbf {r} \wedge \mathbf {s} )=B(\mathbf {r} ,\mathbf {s} )}

dóndersΛ2V{\displaystyle \mathbf {r} \wedge \mathbf {s} \in \Lambda ^{2}V}es el producto exterior der{\displaystyle \mathbf {r} }ys{\displaystyle \mathbf {s} }.

Tomando la L sostenida de L obtenemos

(Ωr)s=L(rs){\displaystyle (\Omega \mathbf {r} )\cdot \mathbf {s} =L^{\sharp }\cdot (\mathbf {r} \wedge \mathbf {s} )}

Presentandoω:=(L){\displaystyle {\boldsymbol {\omega }}:={\star }(L^{\sharp })}, como el dual de Hodge de L , y aplicando la definición del dual de Hodge dos veces suponiendo que el 3-vector unitario preferido es1{\displaystyle \star 1}

(Ωr)s=((L)rs)=(ωrs)=(ωr)s=(ω×r)s,{\displaystyle (\Omega \mathbf {r} )\cdot \mathbf {s} ={\star }({\star }(L^{\sharp })\wedge \mathbf {r} \wedge \mathbf {s} )={\star }({\boldsymbol {\omega }}\wedge \mathbf {r} \wedge \mathbf {s} )={\star }({\boldsymbol {\omega }}\wedge \mathbf {r} )\cdot \mathbf {s} =({\boldsymbol {\omega }}\times \mathbf {r} )\cdot \mathbf {s} ,}

dónde

ω×r:=(ωr){\displaystyle {\boldsymbol {\omega }}\times \mathbf {r} :={\star }({\boldsymbol {\omega }}\wedge \mathbf {r} )}

por definición.

Porques{\displaystyle \mathbf {s} }es un vector arbitrario, de la no degeneración del producto escalar se deduce

Ωr=ω×r{\displaystyle \Omega \mathbf {r} ={\boldsymbol {\omega }}\times \mathbf {r} }

Velocidad angular como un campo vectorial

Dado que el tensor de velocidad angular de giro de un cuerpo rígido (en su sistema de referencia en reposo ) es una transformación lineal que asigna posiciones a velocidades (dentro del cuerpo rígido), puede considerarse como un campo vectorial constante . En particular, la velocidad angular de giro es un campo vectorial de Killing perteneciente a un elemento del álgebra de Lie SO(3) del grupo de rotación tridimensional SO(3) .

Además, se puede demostrar que el campo vectorial de velocidad angular de giro es exactamente la mitad del rotacional del campo vectorial de velocidad lineal v ( r ) del cuerpo rígido. En símbolos,

ω=12×v{\displaystyle {\boldsymbol {\omega }}={\frac {1}{2}}\nabla \times \mathbf {v} }

Consideraciones sobre cuerpos rígidos

Posición del punto P ubicado en el cuerpo rígido (mostrado en azul). R i es la posición con respecto al sistema de referencia del laboratorio, centrado en O , y r i es la posición con respecto al sistema de referencia del cuerpo rígido, centrado en O . El origen del sistema de referencia del cuerpo rígido se encuentra en la posición vectorial R desde el sistema de referencia del laboratorio.

Las mismas ecuaciones para la velocidad angular se pueden obtener considerando un cuerpo rígido en rotación . Aquí no se supone que el cuerpo rígido gire alrededor del origen, sino que puede suponerse que gira alrededor de un punto arbitrario que se mueve con una velocidad lineal V ( t ) en cada instante.

Para obtener las ecuaciones, conviene imaginar un cuerpo rígido unido a los marcos y considerar un sistema de coordenadas fijo respecto a dicho cuerpo. A continuación, estudiaremos las transformaciones de coordenadas entre este sistema y el marco fijo del laboratorio .

Como se muestra en la figura de la derecha, el origen del sistema de laboratorio está en el punto O , el origen del sistema de cuerpo rígido está en O y el vector de O a O es R. Una partícula ( i ) en el cuerpo rígido se encuentra en el punto P y la posición vectorial de esta partícula es R i en el sistema de referencia del laboratorio, y en la posición r i en el sistema de referencia del cuerpo. Se observa que la posición de la partícula se puede escribir como:

Ri=R+ri{\displaystyle \mathbf {R} _{i}=\mathbf {R} +\mathbf {r} _{i}}

La característica definitoria de un cuerpo rígido es que la distancia entre dos puntos cualesquiera de un cuerpo rígido no cambia en el tiempo. Esto significa que la longitud del vectorri{\displaystyle \mathbf {r} _{i}}es inmutable. Por el teorema de rotación de Euler , podemos reemplazar el vectorri{\displaystyle \mathbf {r} _{i}}conRrio{\displaystyle {\mathcal {R}}\mathbf {r} _{io}}dóndeR{\displaystyle {\mathcal {R}}}es una matriz de rotación de 3×3 yrio{\displaystyle \mathbf {r} _{io}}es la posición de la partícula en un punto fijo en el tiempo, digamos t = 0. Este reemplazo es útil, porque ahora es solo la matriz de rotación.R{\displaystyle {\mathcal {R}}}que está cambiando en el tiempo y no el vector de referenciario{\displaystyle \mathbf {r} _{io}}, a medida que el cuerpo rígido gira alrededor del punto O . Además, dado que las tres columnas de la matriz de rotación representan los tres versores de un sistema de referencia que gira junto con el cuerpo rígido, cualquier rotación alrededor de cualquier eje se vuelve ahora visible, mientras que el vectorri{\displaystyle \mathbf {r} _{i}}no rotaría si el eje de rotación fuera paralelo a él, y por lo tanto solo describiría una rotación alrededor de un eje perpendicular a él (es decir, no vería la componente del pseudovector de velocidad angular paralela a él, y solo permitiría el cálculo de la componente perpendicular a él). La posición de la partícula ahora se escribe como:

Ri=R+Rrio{\displaystyle \mathbf {R} _{i}=\mathbf {R} +{\mathcal {R}}\mathbf {r} _{io}}

Al tomar la derivada con respecto al tiempo se obtiene la velocidad de la partícula:

Vi=V+dRdtrio{\displaystyle \mathbf {V} _{i}=\mathbf {V} +{\frac {d{\mathcal {R}}}{dt}}\mathbf {r} _{io}}

donde V i es la velocidad de la partícula (en el sistema de referencia del laboratorio) y V es la velocidad de O (el origen del sistema de referencia del cuerpo rígido). Dado queR{\displaystyle {\mathcal {R}}}es una matriz de rotación, su inversa es su transpuesta. Entonces sustituimosI=RTR{\displaystyle {\mathcal {I}}={\mathcal {R}}^{\text{T}}{\mathcal {R}}}:

Vi=V+dRdtIrio{\displaystyle \mathbf {V} _{i}=\mathbf {V} +{\frac {d{\mathcal {R}}}{dt}}{\mathcal {I}}\mathbf {r} _{io}}
Vi=V+dRdtRTRrio{\displaystyle \mathbf {V} _{i}=\mathbf {V} +{\frac {d{\mathcal {R}}}{dt}}{\mathcal {R}}^{\text{T}}{\mathcal {R}}\mathbf {r} _{io}}
Vi=V+dRdtRTri{\displaystyle \mathbf {V} _{i}=\mathbf {V} +{\frac {d{\mathcal {R}}}{dt}}{\mathcal {R}}^{\text{T}}\mathbf {r} _{i}}

o

Vi=V+Ωri{\displaystyle \mathbf {V} _{i}=\mathbf {V} +\Omega \mathbf {r} _{i}}

dóndeΩ=dRdtRT{\displaystyle \Omega ={\frac {d{\mathcal {R}}}{dt}}{\mathcal {R}}^{\text{T}}}es el tensor de velocidad angular anterior.

Se puede demostrar que esta es una matriz antisimétrica , por lo que podemos tomar su dual para obtener un pseudovector tridimensional que es precisamente el vector de velocidad angular anterior.ω{\displaystyle {\boldsymbol {\omega }}}:

ω=[ωincógnita,ωy,ωz]{\displaystyle {\boldsymbol {\omega }}=[\omega _{x},\omega _{y},\omega _{z}]}

Sustituyendo ω por Ω en la expresión de velocidad anterior y reemplazando la multiplicación de matrices por un producto vectorial equivalente:

Vi=V+ω×ri{\displaystyle \mathbf {V} _{i}=\mathbf {V} +{\boldsymbol {\omega }}\times \mathbf {r} _{i}}

Se puede observar que la velocidad de un punto en un cuerpo rígido se puede dividir en dos términos: la velocidad de un punto de referencia fijo en el cuerpo rígido más el término de producto vectorial que involucra la velocidad angular orbital de la partícula con respecto al punto de referencia. Esta velocidad angular es lo que los físicos llaman la "velocidad angular de giro" del cuerpo rígido, en contraposición a la velocidad angular orbital del punto de referencia O ' alrededor del origen O.

Consistencia

Hemos supuesto que el cuerpo rígido gira alrededor de un punto arbitrario. Debemos demostrar que la velocidad angular de giro previamente definida es independiente del origen, lo que significa que es una propiedad intrínseca del cuerpo rígido en rotación. (Nótese el marcado contraste con la velocidad angular orbital de una partícula puntual , que sí depende del origen).

Demostrar la independencia de la velocidad angular de giro con respecto a la elección del origen.

Vea el gráfico de la derecha: El origen del sistema de referencia del laboratorio es O , mientras que O1 y O2 son dos puntos fijos en el cuerpo rígido, cuya velocidad esv1{\displaystyle \mathbf {v} _{1}}yv2{\displaystyle \mathbf {v} _{2}}respectivamente. Supongamos que la velocidad angular con respecto a O 1 y O 2 esω1{\displaystyle {\boldsymbol {\omega }}_{1}}yω2{\displaystyle {\boldsymbol {\omega }}_{2}}respectivamente. Dado que los puntos P y O2 tienen solo una velocidad,

v1+ω1×r1=v2+ω2×r2{\displaystyle \mathbf {v} _{1}+{\boldsymbol {\omega }}_{1}\times \mathbf {r} _{1}=\mathbf {v} _{2}+{\boldsymbol {\omega }}_{2}\times \mathbf {r} _{2}}
v2=v1+ω1×r=v1+ω1×(r1r2){\displaystyle \mathbf {v} _{2}=\mathbf {v} _{1}+{\boldsymbol {\omega }}_{1}\times \mathbf {r} =\mathbf {v} _{1}+{\boldsymbol {\omega }}_{1}\times (\mathbf {r} _{1}-\mathbf {r} _{2})}

Los dos resultados anteriores que

(ω2ω1)×r2=0{\displaystyle ({\boldsymbol {\omega }}_{2}-{\boldsymbol {\omega }}_{1})\times \mathbf {r} _{2}=0}

Dado que el punto P (y por lo tantor2{\displaystyle \mathbf {r} _{2}}) es arbitrario, por lo tanto se deduce que

ω1=ω2{\displaystyle {\boldsymbol {\omega }}_{1}={\boldsymbol {\omega }}_{2}}

Si el punto de referencia es el eje instantáneo de rotación, la expresión de la velocidad de un punto en el cuerpo rígido solo incluirá el término de velocidad angular. Esto se debe a que la velocidad del eje instantáneo de rotación es cero. Un ejemplo de eje instantáneo de rotación es la bisagra de una puerta. Otro ejemplo es el punto de contacto de un cuerpo rígido esférico (o, más generalmente, convexo) que rueda libremente.

Referencias

  1. Rotaciones y momento angular en la página de Mecánica Clásica del sitio web de John Baez , especialmente las preguntas 1 y 2.