Autor |
Spektrale Finite Elemente |
|
civilengineer
Ehemals Aktiv  Dabei seit: 25.05.2011 Mitteilungen: 1762
Herkunft: Bayern, München
 |
Hallo,
bei der Lösung der Wellengleichung bin ich in der Literatur auf eine Matrizengleichung gestossen, die folgendermassen ausschaut:
 
(\kappa^2 * [K_2] + \kappa * [K_1] + ([K_0] - \omega^2 * [M])) * \v^> = 0^>
k ist die Wellenzahl. K0, K1 und K2 sind Matrizen, die mit der Steifigkeit des Systems zu tun haben. Omega ist eine Frequenz, die man vorgibt, um nach k (Eigenwerte) und den Eigenvektoren aufzulösen. M ist die Massenmatrix.
Ich würde gerne wissen, woher genau diese Gleichung kommt?! Wie hängt diese Matrizengleichung mit der Wellendifferentialgleichung zusammen? Und wie rechnet man die Schwingformen (die freien Vibrationen) mithilfe von den Eigenwerten und Eigenvektoren dieser obigen Gleichung aus?
Ich würde mich freuen, falls jemand, der sich mit der Lösung der Wellengleichung mithilfe der Fouriertransformation auskennt, mir einen Ansatz dazu machen könnte.
Vielen Dank für eure Hilfe im Voraus!
Edit:
Nebenfrage
Edit:
Nebenfrage
|
Notiz Profil
Quote
Link |
civilengineer
Ehemals Aktiv  Dabei seit: 25.05.2011 Mitteilungen: 1762
Herkunft: Bayern, München
 |     Beitrag No.1, vom Themenstarter, eingetragen 2017-06-03
|
Im Prinzip geht es um die Konstruktion einer vollständigen Lösung von irgendeiner Wellengleichung im Frequenzbereich und die Rücktransformation dieser Lösung in den Zeitbereich...
|
Notiz Profil
Quote
Link |
Spock
Senior  Dabei seit: 25.04.2002 Mitteilungen: 8144
Herkunft: Schi'Kahr/Vulkan
 |     Beitrag No.2, eingetragen 2017-06-03
|
Hallo!
2017-05-30 10:08 - civilengineer im Themenstart schreibt:
...
bei der Lösung der Wellengleichung bin ich in der Literatur auf eine Matrizengleichung gestossen
...
Deine Frage wird wahrscheinlich klarer, wenn Du die "Literatur" dazu angibst, Wellengleichungen gibt es wie Sand am Meer, :-)
Gruß
Juergen
|
Notiz Profil
Quote
Link |
civilengineer
Ehemals Aktiv  Dabei seit: 25.05.2011 Mitteilungen: 1762
Herkunft: Bayern, München
 |     Beitrag No.3, vom Themenstarter, eingetragen 2017-06-03
|
Hi Juergen,
danke für deine Antwort. Hier ist die Gleichung, die ich gemeint habe, mit Ansatz in dieser PDF-Datei: Chapter 8: Waveguide Finite Element Method
Ich habe genauso ein Gleichungssystem (8.27) vor mir liegen und ich hab nach den Eigenwerten und den Eigenvektoren gelöst. Ich weiss mittlerweile, dass man für die Verschiebungen den Ansatz macht ... A sind die Amplituden der Wellen ... Ich kann mit der Lösung gerade trotzdem nichts anfangen ... Ich dachte die Lösung (des Eigenwertproblems für eine gegebene Frequenz ) hängt nicht nur von x sondern auch von der Zeit t ab ...
Im Prinzip sind die Eigenvektoren, die man aus diesem Gleichungssystem rausbekommt genau diese Amplituden, richtig? Wenn man also die Eigenvektoren aus der Lösung mit multipliziert, dann hat man schon eine Lösung, korrekt? ...
Das grösste Problem dabei ist, dass sowohl die Eigenvektoren (A) als auch die zugehörigen Eigenwerte komplexe Zahlen sind. Die Verschiebung somit auch?!?! ... Die Verschiebung sollte aber eine reelle Grösse sein ... ?
Viele Grüße
CE
|
Notiz Profil
Quote
Link |
civilengineer
Ehemals Aktiv  Dabei seit: 25.05.2011 Mitteilungen: 1762
Herkunft: Bayern, München
 |     Beitrag No.4, vom Themenstarter, eingetragen 2017-06-07
|
Soll ich einfach den Realteil der Verformung v(x) nehmen? ... Was ist mit der Zeit t? ...
|
Notiz Profil
Quote
Link |
civilengineer
Ehemals Aktiv  Dabei seit: 25.05.2011 Mitteilungen: 1762
Herkunft: Bayern, München
 |     Beitrag No.5, vom Themenstarter, eingetragen 2017-06-07
|
Der Lösungsansatz in Gleichung (8.29)... Der hängt nur von x ab... Wo bleibt die Abhängigkeit von der Zeit t? Ausserdem sind sowohl die Amplituden als auch die Wellenzahlen (die Eigenwerte) komplexe Zahlen... Wo bleibt die reale Verschiebung?!...
|
Notiz Profil
Quote
Link |
civilengineer
Ehemals Aktiv  Dabei seit: 25.05.2011 Mitteilungen: 1762
Herkunft: Bayern, München
 |     Beitrag No.6, vom Themenstarter, eingetragen 2017-06-09
|
Alle machen Informatik jetzt, niemand will Physik machen
|
Notiz Profil
Quote
Link |
civilengineer
Ehemals Aktiv  Dabei seit: 25.05.2011 Mitteilungen: 1762
Herkunft: Bayern, München
 |     Beitrag No.7, vom Themenstarter, eingetragen 2017-06-10
|
Notiz Profil
Quote
Link |
civilengineer
Ehemals Aktiv  Dabei seit: 25.05.2011 Mitteilungen: 1762
Herkunft: Bayern, München
 |     Beitrag No.8, vom Themenstarter, eingetragen 2017-06-10
|
Notiz Profil
Quote
Link |
civilengineer
Ehemals Aktiv  Dabei seit: 25.05.2011 Mitteilungen: 1762
Herkunft: Bayern, München
 |     Beitrag No.9, vom Themenstarter, eingetragen 2017-06-13
|
Notiz Profil
Quote
Link |
civilengineer
Ehemals Aktiv  Dabei seit: 25.05.2011 Mitteilungen: 1762
Herkunft: Bayern, München
 |     Beitrag No.10, vom Themenstarter, eingetragen 2017-06-13
|
Notiz Profil
Quote
Link |
Spock
Senior  Dabei seit: 25.04.2002 Mitteilungen: 8144
Herkunft: Schi'Kahr/Vulkan
 |     Beitrag No.11, eingetragen 2017-06-14
|
Hallo!
 
Normalerweise wird man die zweite Zeitableitung (und damit die Zeitabhängigkeit) in der Wellengleichung los, wenn man vom Zeitbereich in den Frequenzbereich fouriertransformiert. Unter Umständen hat der Verfasser auch einen stationären Lösungsansatz der Form v^>(x^>,t)=v^>^~(x^>) exp(-\ii\.\w t) gewählt. Beides würde das \w^2 in den Gleichungen erklären.
Aber so richtig verstanden habe ich die Gleichung (8.24) und (8.27) in Deinem Beitrag No.3 noch nicht, zumal auch Dein Link zur gesamten Datei bei mir nicht funktioniert.
Ist Dein Problem, daß Du den Artikel verstehen mußt, oder wie lautet Deine Wellengleichung bevor Du bei der Gleichung landest, die Du im Themenstart hingeschrieben hast?
Gruß
Juergen
|
Notiz Profil
Quote
Link |
civilengineer
Ehemals Aktiv  Dabei seit: 25.05.2011 Mitteilungen: 1762
Herkunft: Bayern, München
 |     Beitrag No.12, vom Themenstarter, eingetragen 2017-06-15
|
Hallo Jürgen,
danke für deine Antwort!
Wie man auf dieses Eigenwertproblem kommt, ist mir auch noch ein Rätsel.
Ich weiss nur, dass es dieses Problem gibt und dass man für die Transformation der Gleichung (8.27) den Lösungsansatz in Gleichung (8.29) wählt, um auf ein quadratisches Eigenwertproblem zu kommen.
Die Lösung dieses Eigenwertproblems liefert sowohl die Eigenwerte (Wellenzahlen, kappa) als auch die (Rechts-)Eigenvektoren (Amplituden der Wellenschwingungen).
Sowohl die Eigenwerte als auch die Eigenvektoren können komplex sein. Die Eigenwerte liegen genau doppeltsymmetrisch zu den Achsen der komplexen Zahlenebene.
Wie kann ich aus diesen (angeblichen) Lösungen eine Wellenschwingung genau konstruieren? Ich würde gerne die Wellen animieren. Ich habe mir dazu ein Mathebuch (Spectral Theory of Guided Waves von A.S Silbergleit,Y Kopilevich) von der Bib bestellt und hoffe, dass die Erklärung dort zu finden ist?!
Viele Grüße
CE
|
Notiz Profil
Quote
Link |
civilengineer
Ehemals Aktiv  Dabei seit: 25.05.2011 Mitteilungen: 1762
Herkunft: Bayern, München
 |     Beitrag No.13, vom Themenstarter, eingetragen 2017-06-19
|
Notiz Profil
Quote
Link |
civilengineer
Ehemals Aktiv  Dabei seit: 25.05.2011 Mitteilungen: 1762
Herkunft: Bayern, München
 |     Beitrag No.14, vom Themenstarter, eingetragen 2017-06-19
|
Notiz Profil
Quote
Link |
Spock
Senior  Dabei seit: 25.04.2002 Mitteilungen: 8144
Herkunft: Schi'Kahr/Vulkan
 |     Beitrag No.15, eingetragen 2017-06-20
|
Hallo!
2017-06-19 16:19 - civilengineer in Beitrag No. 13 schreibt:
...
 
Ich habe die Lösungen für Gleichung 15 rausgefunden. Das sind dann die U^~(\kappa) Lösungen? Muss ich jetzt die Formel U^~(x) = 1/(2 \pi) int(U^~(\kappa) e^(-i*\kappa*x),x,-\inf ,+\inf) anwenden, um die Lösungen in Abhängigkeit von x rausfinden zu können?
...
Ja!
Gruß
Juergen
|
Notiz Profil
Quote
Link |
civilengineer
Ehemals Aktiv  Dabei seit: 25.05.2011 Mitteilungen: 1762
Herkunft: Bayern, München
 |     Beitrag No.16, vom Themenstarter, eingetragen 2017-06-20
|
Hallo Jürgen,
ok, da habe ich meine Probleme, diese Formel anzuwenden bzw. überhaupt zu verstehen ...
Und zwar ist es so: Ein Eigenwert und der zugehörige Eigenvektor sind immer konstant und größtenteils komplex...
Also wenn der Eigenvektor konstant ist, dann kann ich ihn aus der Integration herausnehmen: . Korrekt?
Dann muss ich nur die Exponentialfunktion von bis integrieren: Korrekt? ... Ich weiss nicht, wie das weitergehen soll ... cauchy residue theorem?!?! ...
Die Integration über von bis macht gerade nicht viel Sinn für mich ... Es gibt reelle, imaginäre und komplexe Eigenwerte. Sie liegen doppeltsymmetrisch zu den Achsen der komplexen Zahlenebene...
Ist diese Formel so zu interpretieren, dass man einfach über alle vorhandenen Eigenwerte und Eigenvektoren aufsummiert? ...
Kann man anhand eines Eigenwertes mit seinem zugehörigen Eigenvektor keine Formel im Zeitbereich herleiten, die die Verschiebungen der Welle (Wellenschwingungen/Vibrationen/...) in Abhängigkeit von x und t beschreibt? ...
Viele Grüße
CE
|
Notiz Profil
Quote
Link |
civilengineer
Ehemals Aktiv  Dabei seit: 25.05.2011 Mitteilungen: 1762
Herkunft: Bayern, München
 |     Beitrag No.17, vom Themenstarter, eingetragen 2017-06-20
|
In einem wissenschaftlichen Artikel (Guided wave dispersion curves for a bar with an arbitrary cross-section, a rod and rail example) habe ich das Folgende gefunden:
Bei mir ist N(x,y) = 1 erstmal.
Viele Grüße
CE
|
Notiz Profil
Quote
Link |
civilengineer
Ehemals Aktiv  Dabei seit: 25.05.2011 Mitteilungen: 1762
Herkunft: Bayern, München
 |     Beitrag No.18, vom Themenstarter, eingetragen 2017-06-21
|
In einem Paper (On the study of vibration of a supported railway rail using the semi-analytical finite element method - Wenxu Li. Richard A. Dwight. Tieling Zhang.) habe ich das Folgende gefunden:
...
Viele Grüße
CE
|
Notiz Profil
Quote
Link |
civilengineer
Ehemals Aktiv  Dabei seit: 25.05.2011 Mitteilungen: 1762
Herkunft: Bayern, München
 |     Beitrag No.19, vom Themenstarter, eingetragen 2017-06-22
|
Notiz Profil
Quote
Link |
civilengineer
Ehemals Aktiv  Dabei seit: 25.05.2011 Mitteilungen: 1762
Herkunft: Bayern, München
 |     Beitrag No.20, vom Themenstarter, eingetragen 2017-06-25
|
Notiz Profil
Quote
Link |
civilengineer
Ehemals Aktiv  Dabei seit: 25.05.2011 Mitteilungen: 1762
Herkunft: Bayern, München
 |     Beitrag No.21, vom Themenstarter, eingetragen 2017-07-05
|
Notiz Profil
Quote
Link |