Die Rayleigh-Ritz-Methode ist eine numerische Methode zum Auffinden von Näherungen an Eigenwertgleichungen, die sich nur schwer analytisch lösen lassen, insbesondere im Zusammenhang mit der Lösung physikalischer Grenzwertprobleme, die sich als Matrixdifferentialgleichungen ausdrücken lassen. Es wird im Maschinenbau verwendet, um die Eigenmoden eines physikalischen Systems zu approximieren, beispielsweise um die Resonanzfrequenzen einer Struktur zu ermitteln, um eine angemessene Dämpfung zu steuern. Der Name ist eine verbreitete Fehlbezeichnung zur Beschreibung der Methode, die besser als Ritz-Methode oder Galerkin-Methode bezeichnet wird. [ Zitat erforderlich Diese Methode wurde 1909 von Walther Ritz erfunden. aber es hat eine gewisse Ähnlichkeit mit dem Rayleigh-Quotienten, und so bleibt die Fehlbezeichnung bestehen. [ Benötigtes Zitat
Beschreibung der Methode [
Die Rayleigh-Ritz-Methode ermöglicht die Berechnung von Ritz-Paaren
die die Lösungen für das Eigenwertproblem annähern [1]
-

wobei
.
Das Verfahren ist wie folgt: [2]
- Berechne eine orthonormale Basis
der sich dem -Eigenraum nähert, der [19459006entspricht] m Eigenvektoren
- Berechnen
- Berechne die Eigenwerte der R-Lösung
- Bilden Sie die Ritz-Paare
Die Genauigkeit einer solchen Annäherung kann immer über
Die Methode in der Variationsrechnung [ edit
In dieser Technik approximieren wir das Variationsproblem und erhalten ein endliches Dimensionsproblem. Beginnen wir also mit dem Problem der Suche nach einer Funktion
![y (x) [19659020] das ein Integral extremisiert <span class=](https://wikimedia.org/api/rest_v1/media/math/render/svg/e871993bfd131a8b0c3591c26084cf8171a74dcd)
. Angenommen, wir können y (x) durch eine lineare Kombination bestimmter linear unabhängiger Funktionen des Typs approximieren:
wobei
sind Konstanten, die durch ein Variationsverfahren wie eines, das nachstehend beschrieben wird, zu bestimmen sind.
Die Auswahl der Approximationsfunktionen
Die Verwendung ist mit Ausnahme der folgenden Überlegungen willkürlich:
a) Wenn das Problem Randbedingungen wie feste Endpunkte aufweist, dann
wird ausgewählt, um die Randbedingungen des Problems zu erfüllen, und alle anderen
verschwinden an der Grenze.
b) Wenn die Form der Lösung bekannt ist, dann
kann so gewählt werden, dass
wird diese Form haben.
Die Erweiterung von
in etwa Funktionen ersetzen das Variationsproblem der Extremisierung des Funktionsintegrals
auf ein Problem beim Finden einer Menge von Konstanten
das
. Wir können dies jetzt lösen, indem wir die partiellen Ableitungen auf Null setzen. Für jeden Wert von i
Das Verfahren besteht darin, zuerst eine anfängliche Schätzung von
durch die Approximation
. Als nächstes die Näherung
wird verwendet (wobei
erneut bestimmt wird) . Der Prozess wird mit
als dritte Annäherung und so weiter. In jeder Phase sind die folgenden zwei Punkte wahr:
- Im i-ten Stadium sind die Begriffe
wird neu bestimmt [19659049] Die Annäherung an das
Stadium
wird nicht schlechter sein als der Annäherung an den
Stadium
Konvergenz der Prozedur bedeutet, dass, wenn i gegen Unendlich tendiert, die Approximation gegen die exakte Funktion tendiert
das ein Integral
.
In vielen Fällen verwendet man einen vollständigen Satz von Funktionen e. G. Polynome oder Sinus und Cosinus. Eine Reihe von Funktionen
heißt vollständig über [a, b]wenn für jede integrierbare Riemann-Funktion
es gibt einen Satz von Koeffizientenwerten
das
.
Die oben beschriebene Prozedur kann auf Fälle mit mehr als einer unabhängigen Variablen erweitert werden.
Anwendungen im Maschinenbau [
Die Rayleigh-Ritz-Methode wird im Maschinenbau häufig verwendet, um die ungefähren realen Resonanzfrequenzen von Systemen mit mehreren Freiheitsgraden, wie z. B. Federn, zu ermitteln Massensysteme oder Schwungräder auf einer Welle mit unterschiedlichem Querschnitt. Es ist eine Erweiterung von Rayleighs Methode. Es kann auch zum Auffinden von Knicklasten und Nachknickverhalten für Stützen verwendet werden.
Betrachten wir den Fall, in dem wir die Resonanzfrequenz der Schwingung eines Systems ermitteln wollen. Schreiben Sie zuerst die Schwingung in die Form,
mit einer unbekannten Modusform
![{ displaystyle Y (x)} [19659037]. Bestimmen Sie als Nächstes die Gesamtenergie des Systems, die aus einem kinetischen und einem potenziellen Energieterm besteht. Der kinetische Energieterm beinhaltet das Quadrat der Zeitableitung von <span class=](https://wikimedia.org/api/rest_v1/media/math/render/svg/fca9de77ea29879f48f3f08c469fb73b8945a094)
und erhält damit einen Faktor von
. So können wir die Gesamtenergie des Systems berechnen und in folgender Form ausdrücken:
Durch Energieerhaltung muss die mittlere kinetische Energie gleich der mittleren potentiellen Energie sein. Somit,
der auch als Rayleigh-Quotient bezeichnet wird. Wenn wir also die Modenform
![{ displaystyle Y (x)} [19659037]könnten wir berechnen <span class=](https://wikimedia.org/api/rest_v1/media/math/render/svg/fca9de77ea29879f48f3f08c469fb73b8945a094)
und
und erhalten die Eigenfrequenz. Die Modenform ist uns jedoch noch nicht bekannt. Um dies zu finden, können wir
als Kombination einiger Näherungsfunktionen
wobei
c 1 c 2 ⋯ c N { displaystyle c_ {1}, c_ {2}, cdots, c_ {N}}
sind Konstanten bestimmt werden. Wenn wir eine zufällige Menge von
c 1 c 2 ⋯ c [19659159] N { displaystyle c_ {1}, c_ {2}, cdots, c_ {N}}
wird eine Überlagerung der tatsächlichen Eigenmoden des Systems beschrieben. Wenn wir jedoch nach
c 1 c 2 ⋯ c N [19659068suchen] { displaystyle c_ {1}, c_ {2}, cdots, c_ {N}}
wie z dass die Eigenfrequenz
ω 2 { displaystyle omega ^ {2}}
minimiert wird, dann wird der durch diesen Satz von beschriebene Modus
c 1 c 2 ⋯ c N { displaystyle c_ {1}, c_ {2}, cdots, c_ {N}}
liegen nahe am niedrigsten möglicher tatsächlicher Eigenmodus des Systems. Somit findet dies die niedrigste Eigenfrequenz. Wenn wir Eigenmoden finden, die orthogonal zu dieser angenäherten niedrigsten Eigenmode sind, können wir auch die nächsten Eigenfrequenzen finden.
Im Allgemeinen können wir
A [ Y ( x ) { displaystyle A [Y(x)] ausdrücken. }
und
B [ Y ( x ) { displaystyle B [Y(x)]}
als Sammlung von Begriffen mit quadratischen Koeffizienten
c i { displaystyle c_ {i}}
:
B Y ( x ] 19659080] = ∑ i ∑ j c i c j K 19659055] j = c T K c { displaystyle B [Y(x)] = sum _ {i} sum _ {j} c_ {i} c_ {j} K_ {ij} = { bf {c ^ {T} Kc}}
A Y ( x ] 19659080] = ∑ i ∑ j c i c j i 19659055] j = c T M c { displaystyle A [Y(x)] = sum _ {i} sum _ {j} c_ {i} c_ {j} M_ {ij} = { bf {c ^ {T} Mc}}
Die Minimierung von
ω 2 { displaystyle omega ^ {2}}
wird:
∂ ω 2 ∂ 19659273] c i = ∂ ∂ c i c T K c c T M c = 0 { displaystyle { partial omega ^ {2} over partial c_ {i}} = { partial over partial c_ {i}} { frac { bf {c ^ {T} Kc}} { bf {c ^ {T} Mc}} = 0}
Lösen Sie dies,
c T M c 19659830] ∂ c T K c ∂ c – c T K c ∂ c T M c ∂ c = 0 { displaystyle { bf {{ c ^ {T} Mc} { partiell { bf {c ^ {T} Kc}} über partiell c} – { bf {{c ^ {T} Kc} { partiell { bf {c ^ {T} Mc}} over partial c} = 0}}}}
K c – c T K c c T M c M c = 0 { displaystyle { bf { {Kc} – { frac { bf {c ^ {T} Kc}} { bf {c ^ {T} Mc}} { bf {{Mc} = 0}}}}
K c – ω [19659874] 2 M c = 0 { displaystyle { bf {{Kc} – omega ^ {2} { bf {{Mc} = 0}} }}}
Für eine nicht triviale Lösung von c muss die Determinante des Matrixkoeffizienten von c Null sein.
det ( K – ω 2 [19659888] M ) = 0 { displaystyle det ({ bf {{K} – omega ^ {2} { bf {M}}}) = 0}
Dies ergibt eine Lösung für die ersten N Eigenfrequenzen und Eigenmoden des Systems, wobei N die Anzahl der Approximationsfunktionen ist.
Einfacher Fall eines Systems mit doppelter Federmasse
In der folgenden Diskussion wird der einfachste Fall verwendet, in dem das System zwei konzentrierte Federn und zwei konzentrierte Massen und nur zwei Modi aufweist Formen werden angenommen. Daher M = m 1 m 2 und K = k [19659898] 1 k 2 .
Für das System wird eine Modenform mit zwei Begriffen angenommen, von denen einer mit einem Faktor B, z. Y = [1, 1] + B [1, −1].
Einfache harmonische Bewegungstheorie besagt, dass die Geschwindigkeit zum Zeitpunkt der Ablenkung Null ist, die Winkelfrequenz
ω { displaystyle omega}
multipliziert mit der Ablenkung (y) zum Zeitpunkt von maximale Auslenkung. In diesem Beispiel beträgt die kinetische Energie (KE) für jede Masse
1 2 ω 2 Y 1 2 m 1 { displaystyle { frac {1} {2}} omega ^ {2} Y_ {1} ^ {2} m_ {1}}
usw., und die potentielle Energie (PE) für jede Feder beträgt
1 2 [19659922] k 1 Y 1 2 { displaystyle { frac {1} {2}} k_ {1} Y_ {1} ^ {2}} [19659928] { frac {1} {2}} k_ {1} Y_ {1} ^ {2} “/> etc.
Wir wissen auch, dass das maximale KE ohne Dämpfung gleich dem maximalen PE ist. Somit,
-
∑ i = 1 2 ( 1 2 ω 2 Y i 2 M i ) = ∑ i = 1 2 ( 1 2 K i Y i 2 ) { displaystyle sum _ {i = 1} ^ {2 } left ({ frac {1} {2}} omega ^ {2} Y_ {i} ^ {2} M_ {i} right) = sum _ {i = 1} ^ {2} left ({ frac {1} {2}} K_ {i} Y_ {i} ^ {2} right)}

Beachten Sie, dass die Gesamtamplitude der Modenform immer von jeder Seite ausgeglichen wird. Das heißt, die tatsächliche Größe der angenommenen Auslenkung spielt keine Rolle, nur die Form der Mode .
Durch mathematische Manipulationen erhält man dann einen Ausdruck für
ω { displaystyle omega}
nach B, der sich nach B unterscheiden lässt Finden Sie das Minimum, dh wenn
d ω / d B = 0 { displaystyle d omega / dB = 0 }
. Dies ergibt den Wert von B, für den
ω { displaystyle omega}
am niedrigsten ist. Dies ist eine Lösung mit Obergrenze für
ω [196599056] { displaystyle omega}
if
ω { displaystyle omega} [19659906] omega “/> soll die vorhergesagte Grundfrequenz des Systems sein, da die Modenform angenommen wird aber wir haben den niedrigsten Wert dieser oberen Schranke unter unseren Annahmen gefunden, weil B verwendet wird um die optimale 'Mischung' der beiden angenommenen Modenformfunktionen zu finden.
Es gibt viele Tricks bei dieser Methode. Das Wichtigste ist, realistische angenommene Modenformen auszuwählen. Beispielsweise ist es im Fall von Strahlablenkungsproblemen ratsam, eine deformierte Form zu verwenden, die der erwarteten Lösung analytisch ähnlich ist. Ein Viertel kann zu den meisten einfachen Problemen einfach verbundener Balken passen, selbst wenn die Reihenfolge der deformierten Lösung niedriger sein kann. Die Federn und Massen müssen nicht diskret sein, sie können kontinuierlich (oder gemischt) sein, und diese Methode kann leicht in einer Kalkulationstabelle verwendet werden, um die Eigenfrequenzen recht komplexer verteilter Systeme zu ermitteln, wenn Sie das verteilte KE und beschreiben können PE-Begriffe leicht oder brechen die durchgehenden Elemente in diskrete Teile auf.
Diese Methode kann iterativ verwendet werden, indem der vorherigen besten Lösung zusätzliche Modusformen hinzugefügt werden, oder Sie können einen langen Ausdruck mit vielen Bs und vielen Modusformen erstellen und diese dann teilweise differenzieren.
Siehe auch Bearbeiten
Notizen und Referenzen Bearbeiten
Externe Links Bearbeiten ]