¿Cuál es el problema de tener un tensor de inercia que no satisface la desigualdad del triángulo?

Es bien sabido que los tensores de inercia de cuerpo rígido son matrices semidefinidas positivas de 3 por 3, lo que es lo mismo que decir que sus valores propios son todos no negativos.

Un poco menos conocido es el hecho de que esos valores propios también satisfacen las desigualdades del triángulo , lo que significa que la suma de dos valores propios siempre es mayor o igual que el restante. Esto surge de la propia definición del tensor que involucra integrales de una distribución de masa siempre no negativa. Según https://physics.stackexchange.com/a/48273/116038 :

En otras palabras, si un real simétrico definido semipositivo 3 × 3 matriz con valores propios no negativos [...] no satisface la desigualdad triangular (1), no representa una distribución de masa físicamente posible.

Al estimar tensores de inercia usando, por ejemplo, técnicas de regresión, se pueden obtener matrices que no representan cuerpos rígidos físicamente posibles (si no restringimos la solución de regresión a valores físicamente consistentes).

Usando un semidefinido no positivo 3 × 3 matriz simétrica en lugares donde se espera un tensor de inercia, por ejemplo, en algunas formulaciones de cuerpo rígido, en simulación, en esquemas de control, etc., implicaría problemas importantes. Por ejemplo, una matriz semidefinida no positiva utilizada como tensor de inercia daría lugar a energías cinéticas negativas, mi k = 1 2 ω I ω .

Ahora, mi pregunta es: ¿Qué problemas (en formulaciones, simulación, control, etc.), si los hay, pueden surgir si usamos como tensor de inercia alguna matriz que, si bien satisface la condición semidefinida positiva, no satisface las condiciones de desigualdad del triángulo en los valores propios?


Actualización : he comprobado que los "tensores de inercia" que son semidefinidos positivos pero que no satisfacen las desigualdades triangulares aún verifican la ley de conservación de energía con respecto a la energía cinética y el trabajo.

Respuestas (3)

Los principales momentos de inercia de un cuerpo rígido se obtienen a partir de estas ecuaciones:

I 1 = V ( X 2 2 + X 3 2 ) ρ d V
I 2 = V ( X 3 2 + X 1 2 ) ρ d V
I 3 = V ( X 1 2 + X 2 2 ) ρ d V

dónde X = X 1   , y = X 2   , z = X 3 y el tensor de inercia es:

I = [ I 1 0 0 0 I 2 0 0 0 I 3 ]

con

i α = V X α 2 ρ d V > 0   , α = 1 , 2 , 3
de este modo:
I 1 = i 2 + i 3
I 2 = i 3 + i 1
I 3 = i 1 + i 2

y

I 1 + I 2 = i 1 + i 2 + 2 i 3 = I 3 + 2 i 3 > I 3
I 2 + I 3 = 2 i 1 + i 2 + i 3 = I 1 + 2 i 1 > I 1
I 3 + I 1 = i 1 + 2 i 2 + i 3 = I 2 + 2 i 3 > I 2

por lo tanto, la desigualdad del triángulo es una característica física de un tensor de inercia de cuerpo rígido. Si el cuerpo rígido es simétrico, entonces los ejes de simetría son ejes principales y el momento principal de inercia debe obedecer a la desigualdad del triángulo, de lo contrario, no describe el cuerpo rígido que desea describir.

La desigualdad del triángulo para el momento del tensor de inercia en última instancia se deriva de la positividad de la masa de inercia (como se muestra en la respuesta de Eli ). Un cuerpo cuyo tensor de inercia no satisfaga las desigualdades del triángulo necesariamente tendría alguna región donde ρ < 0 .

Entonces podríamos considerar este cuerpo como compuesto por dos sub-cuerpos, uno con masa de inercia estrictamente positiva metro 1 > 0 y uno con masa de inercia estrictamente negativa metro 2 < 0 . Cualquier fuerza interna entre estos cuerpos daría como resultado que estos dos subcuerpos aceleraran en la misma dirección espontáneamente : las fuerzas internas en metro 1 y metro 2 actuarían en direcciones opuestas, según la Tercera Ley de Newton, pero sus vectores de aceleración apuntarían en la misma dirección debido a sus diferentes signos de masa.

Estas "soluciones desbocadas" generalmente se consideran malas. Las únicas lagunas que veo son (a) abandonar la Tercera Ley de Newton, pero eso entristecería a Emmy Noether; o (b) exigir que estos subcuerpos nunca ejerzan ninguna fuerza neta entre sí, lo que sería problemático si queremos que el cuerpo sea rígido.

Es útil observar los tensores de inercia en la base esférica, en los que los tensores de base son tensores propios de rotaciones sobre el eje z. Así que en lugar de considerar 3 x 3 = 9 T i j para yo, j ( X , y , z ) , on tiene 1 + 3 + 5 = 9 estados básicos de la siguiente manera:

T ( 0 ) = 1 3 T i i   I

es esféricamente simétrica y proporcional a la identidad. Como tensor de inercia, corresponde a una esfera.

La parte antisimétrica de un tensor:

T ( 1 ) = A i j 1 2 ( T i j T j i )

se transforma como un vector (muy parecido a un producto cruzado), y NO TIENE LUGAR en los tensores de inercia. Hay 3 componentes.

Finalmente, la parte simétrica libre de trazas tiene 5 componentes:

T ( 2 ) = S i j 1 2 ( T i j + T j i ) ) T ( 0 )

que se transforman como armónicos esféricos: Y yo = 2 metro     ( θ , ϕ ) para metro   ( 2 , 1 , 0 , 1 , 2 ) .

El | metro | = 1 componentes:

T ( 2 , metro = ± 1 )         = S X z ± i S y z

se puede diagonalizar mediante una rotación de coordenadas adecuada dejando:

yo = 2 , metro = 0 :

T ( 2 , 0 )     = 1 2 S z z

Esta parte corresponde a la ablación/alargamiento del objeto; entonces, para un elipsoide, eso sería con la excentricidad con el eje de simetría alineado con el eje z.

yo = 2 , | metro | = 2 :

T ( 2 , metro = ± 2 )           = 1 2 ( S X X S y y ± 2 i S X y )

Para los tensores de inercia, m=2 y m=-2 son iguales, por lo que las componentes cartesianas son reales. Esta parte de los tensores corresponde a cualquier asimetría cilíndrica , por lo que es cero para un elipsoide de revolución.

Gracias por esta otra forma de ver los tensores de inercia, pero no pude obtener más información sobre mi pregunta.