Determinante y adjunto de k−ω2mk−ω2mk-\omega^2m en términos de frecuencias naturales

Se da un sistema mecánico de múltiples grados de libertad descrito por las siguientes matrices y ecuaciones:

  • matriz de masa metro = [ metro 0 0 0 metro 0 0 0 metro / 2 ] ,

  • matriz de rigidez k = [ 2 k k 0 k 2 k k 0 k k ] ,

  • desplazamientos tu = [ tu 1 ( t ) tu 2 ( t ) tu 3 ( t ) ] ,

  • Fuerza externa pag = [ 0 0 pag 0 pecado ( ω t ) ] , y

  • ecuación de movimiento metro tu ¨ + k tu = pag .

Las frecuencias naturales ω i han sido derivados del problema de valores propios det ( k ω i 2 metro ) = 0 lo que llevó a:

ω 1 2 = ( 2 3 ) k metro , ω 2 2 = 2 k metro , ω 3 2 = ( 2 + 3 ) k metro .

Sé que la solución de estado estacionario para tu es

tu = 1 det ( k ω 2 metro ) a d j ( k ω 2 metro ) pag ,

así que tengo que calcular el determinante y el adjunto. De hecho, no tengo ningún problema en hacerlo, pero la solución del libro de texto se ve mucho más elegante que la mía y no sé cómo llegar allí.

Mi solución:

det ( k ω 2 metro ) = 1 2 ( 2 k ω 2 metro ) 3 3 k 3 ,

a d j ( k ω 2 metro ) = [ k 2 k ( 2 k ω 2 metro ) ( 2 k ω 2 metro ) 2 k 2 ] (solo la tercera columna es relevante para la solución)

Solución de libro de texto:

det ( k ω 2 metro ) = 1 2 metro 3 ( ω 1 2 ω 2 ) ( ω 2 2 ω 2 ) ( ω 3 2 ω 2 ) = k 3 ( 1 ω 2 ω 1 2 ) ( 1 ω 2 ω 2 2 ) ( 1 ω 2 ω 3 2 ) ,

a d j ( k ω 2 metro ) = [ 1 2 ( 1 ω 2 / ω 2 2 ) 4 ( 1 ω 2 / ω 2 2 ) 2 1 ] k 2

Mi pregunta es: ¿Cómo puedo obtener la solución del libro de texto que está escrita en términos de frecuencias naturales? Me parece que hay alguna manera de expresar k ω 2 metro principalmente en términos de las frecuencias naturales y luego podría ir desde allí, pero no sé cómo debería hacerlo.

Este es un problema de matemáticas, no de física. Marcar para migrar a math.stackexchange.com
No estoy de acuerdo. Esta es la mecánica interna de (las matemáticas que describen) una parte de la física.
cual libro de texto

Respuestas (2)

El determinante es bastante fácil de calcular. Ya conoces, esencialmente, los valores propios de la matriz de rigidez; con más precisión, conoce los valores propios de la matriz metro 1 k , porque el ω i son ceros de la ecuacion

0 = det ( metro 1 k ω 2 ) .
(Los más estéticos reemplazarían metro 1 k con metro 1 / 2 k metro 1 / 2 para obtener una matriz hermítica, pero no importa). Si expresas el segundo determinante en la base propia correspondiente, obtienes
det ( k metro ω 2 ) = det ( metro ) det ( metro 1 k ω 2 ) = metro 3 2 det ( ω 1 2 ω 2 0 0 0 ω 2 2 ω 2 0 0 0 ω 3 2 ω 2 ) ,
que da la expresión de su libro de texto. De manera más general, esta es una expresión del principio de que el determinante de una matriz es el producto de sus valores propios.


El adjunto, por otro lado, no satisface (que yo sepa) ninguna relación tan agradable; en cualquier caso, es una bestia desagradable con la que lidiar y creo que pocas personas lo sustituyen juiciosamente en la definición k = metro ω 2 2 / 2 de ω 2 en lugar de k .

Sospecho que está abrumado por la plétora de variables que oscurecen la simetría ciclométrica fundamental del problema. Puede escalar todos menos uno de ellos fuera de la matriz simétrica M , parte de cuya inversa (simétrica) está buscando, en realidad, redefiniendo

METRO k ω 2 metro = k [ 2 ( X 1 ) 1 0 1 2 ( X 1 ) 1 0 1 X 1 ] ,
dónde X ω 2 metro 2 k .

La simetría ciclométrica (Mercedes-Benz) de sus valores propios es entonces evidente a partir de su determinante simple, una vez que reconoce la fórmula del triple ángulo para el seno en él, det METRO = k 3 ( X 1 ) ( 4 ( X 1 ) 2 3 ) = k 3 ( X 1 ) ( X 1 3 2 ) ( X 1 + 3 2 ) , mostrando los senos de las 3 raíces de la unidad. Este resultado solo es útil para factorizar el polinomio cúbico a partir del inverso de M , lo que sabes que es posible ya que el adjunto, adj M = det M M 1 es un polinomio en x , y no una razón de polinomios.

Sin embargo, es más fácil encontrar la tercera columna de la inversa (que es todo lo que necesita para su solución) directamente, resolviendo las tres condiciones triviales de solo la tercera columna de METRO   METRO 1 = I . Sabes por la discusión anterior que tu respuesta debe tener un factor inverso de ( X 1 ) ( X 1 3 2 ) ( X 1 + 3 2 ) en él, resumiendo las resonancias del sistema, que se comprueba que tiene; y la transposición de esta tercera columna del adjunto es simplemente k 2 ( 1 , 2 ( 1 X ) , 4 ( 1 X ) 2 1 ) , que es solo la expresión de su libro de texto y, por supuesto, la suya, una vez que se reemplazan las raíces del determinante, los ceros. La última entrada es k 2 ( ( X 1 3 2 ) ( X 1 + 3 2 ) + 2 ) , por supuesto.

Aunque, como argumenta Kyle Kanos en el comentario, este es estrictamente un problema de álgebra lineal, las técnicas de simetría empleadas y la metodología utilizada son el pan y la mantequilla de la física, y bien se podría argumentar que los físicos normalmente son más rápidos, si no mejores, en el manejo aquellos.

Exceso académico PD : si, más allá del alcance de su problema, estaba interesado en el inverso completo, inútil aquí, solo necesita explotar la observación trigonométrica anterior y cambiar las variables por última vez, X 1 pecado ϕ , el dictado esencialmente por la ciclometricidad del problema. Entonces es sencillo observar que el adjunto es realmente elegante,

I = [ 2 pecado ϕ 1 0 1 2 pecado ϕ 1 0 1 pecado ϕ ] [ porque ( 2 ϕ ) pecado ϕ 1 pecado ϕ 2 pecado 2 ϕ 2 pecado ϕ 1 2 pecado ϕ 1 2 porque ( 2 ϕ ) ]   / pecado ( 3 ϕ ) ,
que por supuesto contiene la respuesta anterior.