Determinante und Adjunkte von k−ω2mk−ω2mk-\omega^2m in Bezug auf Eigenfrequenzen

Gegeben ist ein mechanisches System mit mehreren Freiheitsgraden, das durch die folgenden Matrizen und Gleichungen beschrieben wird:

  • Massenmatrix M = [ M 0 0 0 M 0 0 0 M / 2 ] ,

  • Steifigkeitsmatrix k = [ 2 k k 0 k 2 k k 0 k k ] ,

  • Verschiebungen u = [ u 1 ( T ) u 2 ( T ) u 3 ( T ) ] ,

  • äußere Kraft P = [ 0 0 P 0 Sünde ( ω T ) ] , Und

  • Bewegungsgleichung M u ¨ + k u = P .

Die Eigenfrequenzen ω ich wurden aus dem Eigenwertproblem abgeleitet det ( k ω ich 2 M ) = 0 Was zu ... führte:

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

Ich weiß, dass die Steady-State-Lösung für u Ist

u = 1 det ( k ω 2 M ) A D J ( k ω 2 M ) P ,

Also muss ich die Determinante berechnen und adjugieren. Ich habe eigentlich kein Problem damit, aber die Lehrbuchlösung sieht viel eleganter aus als meine und ich weiß nicht, wie ich da hinkomme.

Meine Lösung:

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

A D J ( k ω 2 M ) = [ k 2 k ( 2 k ω 2 M ) ( 2 k ω 2 M ) 2 k 2 ] (nur die dritte Spalte ist lösungsrelevant)

Lehrbuchlösung:

det ( k ω 2 M ) = 1 2 M 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 M ) = [ 1 2 ( 1 ω 2 / ω 2 2 ) 4 ( 1 ω 2 / ω 2 2 ) 2 1 ] k 2

Meine Frage ist: Wie bekomme ich die Lehrbuchlösung, die in Bezug auf die Eigenfrequenzen geschrieben ist? Es scheint mir, als gäbe es eine Möglichkeit, sich auszudrücken k ω 2 M hauptsächlich in Bezug auf die natürlichen Frequenzen und dann könnte ich von dort aus gehen, aber ich weiß nicht, wie ich das machen soll.

Das ist ein mathematisches Problem, kein physikalisches Problem. Markierung zur Migration auf math.stackexchange.com
Ich bin nicht einverstanden. Dies ist die interne Mechanik (die Mathematik, die beschreibt) ein Stück Physik.
welches Lehrbuch?

Antworten (2)

Die Determinante ist ziemlich einfach zu berechnen. Sie kennen im Wesentlichen bereits die Eigenwerte der Steifigkeitsmatrix; Genauer gesagt kennen Sie die Eigenwerte der Matrix M 1 k , weil das ω ich sind Nullstellen der Gleichung

0 = det ( M 1 k ω 2 ) .
(Die Ästhetischeren würden ersetzen M 1 k mit M 1 / 2 k M 1 / 2 um eine hermitesche Matrix zu erhalten, aber egal.) Wenn Sie die zweite Determinante in der entsprechenden Eigenbasis ausdrücken, erhalten Sie
det ( k M ω 2 ) = det ( M ) det ( M 1 k ω 2 ) = M 3 2 det ( ω 1 2 ω 2 0 0 0 ω 2 2 ω 2 0 0 0 ω 3 2 ω 2 ) ,
die Ihrem Lehrbuch Ausdruck verleiht. Allgemeiner ausgedrückt ist dies ein Ausdruck des Prinzips, dass die Determinante einer Matrix das Produkt ihrer Eigenwerte ist.


Das Adjunkt hingegen erfüllt (meines Wissens) keine so schöne Beziehung; auf jeden Fall ist es ein übles Biest, mit dem man sich auseinandersetzen muss, und ich denke, dass nur wenige Leute vernünftigerweise die Definition ersetzen k = M ω 2 2 / 2 von ω 2 anstatt k .

Ich vermute, Sie sind von der Fülle von Variablen überwältigt, die die grundlegende zyklometrische Symmetrie des Problems verschleiern. Sie können alle bis auf einen aus der symmetrischen Matrix M skalieren , von deren (symmetrischer) Umkehrung Sie tatsächlich einen Teil suchen, indem Sie sie neu definieren

M k ω 2 M = k [ 2 ( X 1 ) 1 0 1 2 ( X 1 ) 1 0 1 X 1 ] ,
Wo X ω 2 M 2 k .

Die zyklometrische (Mercedes-Benz) Symmetrie ihrer Eigenwerte ist dann aus ihrer einfachen Determinante ersichtlich, wenn man darin die Dreifachwinkelformel für den Sinus erkennt, det M = k 3 ( X 1 ) ( 4 ( X 1 ) 2 3 ) = k 3 ( X 1 ) ( X 1 3 2 ) ( X 1 + 3 2 ) , die die Sinus der 3 Einheitswurzeln anzeigt. Dieses Ergebnis ist nur nützlich, um das kubische Polynom aus der Inversen von M herauszufaktorisieren , von dem Sie wissen, dass es möglich ist, da das Adjugat adj M = det M M ist 1 ist ein Polynom in x und kein Verhältnis von Polynomen.

Trotzdem ist es einfacher, direkt die 3. Spalte der Umkehrung zu finden (was alles ist, was Sie für Ihre Lösung benötigen), indem Sie nach den drei trivialen Bedingungen nur der 3. Spalte von lösen M   M 1 = ICH . Sie wissen aus der obigen Diskussion, dass Ihre Antwort einen Umkehrfaktor von haben muss ( X 1 ) ( X 1 3 2 ) ( X 1 + 3 2 ) darin die Resonanzen des Systems zusammenfassend, auf die es geprüft wird; und die Transponierung dieser 3. Spalte des Adjugats ist einfach k 2 ( 1 , 2 ( 1 X ) , 4 ( 1 X ) 2 1 ) , was nur der Ausdruck Ihres Lehrbuchs ist, und natürlich auch Ihrer, sobald die Wurzeln der Determinante, die Nullen, eingesteckt sind. Der letzte Eintrag ist k 2 ( ( X 1 3 2 ) ( X 1 + 3 2 ) + 2 ) , Natürlich.

Auch wenn dies, wie Kyle Kanos in dem Kommentar argumentiert, streng genommen ein Problem der linearen Algebra ist, sind die verwendeten Symmetrietechniken und die verwendete Methodik das A und O der Physik, und man könnte durchaus argumentieren, dass Physiker normalerweise schneller, wenn nicht sogar besser damit umgehen können diese.

Akademischer Overkill PS : Wenn Sie über den Rahmen Ihres Problems hinaus an der vollständigen Umkehrung interessiert waren, die hier nutzlos ist, müssen Sie nur die obige trigonometrische Beobachtung ausnutzen und die Variablen ein letztes Mal ändern. X 1 Sünde ϕ , die im Wesentlichen von der Zyklometrie des Problems bestimmt wird. Dann ist es einfach zu beobachten, dass das Adjugat eigentlich elegant ist,

ICH = [ 2 Sünde ϕ 1 0 1 2 Sünde ϕ 1 0 1 Sünde ϕ ] [ cos ( 2 ϕ ) Sünde ϕ 1 Sünde ϕ 2 Sünde 2 ϕ 2 Sünde ϕ 1 2 Sünde ϕ 1 2 cos ( 2 ϕ ) ]   / Sünde ( 3 ϕ ) ,
die natürlich die obige Antwort enthält.