Die Mathe-Redaktion [Home] - Neues auf einen Blick 31.07.2021 01:34 [Neues] [Forum] [Suche] - Registrieren/Login
 
Matroids Matheplanet Forum Index
Moderiert von Berufspenner Ueli rlk MontyPythagoras
Physik » Ingenieurwesen » Spektrale Finite Elemente		Druckversion
Druckversion		Antworten
Antworten
Autor
Universität/Hochschule Spektrale Finite Elemente
civilengineer
Wenig Aktiv Letzter Besuch: vor mehr als 3 Monaten
Dabei seit: 25.05.2011
Mitteilungen: 1762
Wohnort: Bayern, München
Zum letzten BeitragZum nächsten BeitragZum vorigen BeitragZum erstem Beitrag  Themenstart: 2017-05-30


Hallo,

bei der Lösung der Wellengleichung bin ich in der Literatur auf eine Matrizengleichung gestossen, die folgendermassen ausschaut:

fed-Code einblenden

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



Eine Notiz zu diese Forumbeitrag schreiben Notiz   Profil  Quote  Link auf diesen Beitrag Link
civilengineer
Wenig Aktiv Letzter Besuch: vor mehr als 3 Monaten
Dabei seit: 25.05.2011
Mitteilungen: 1762
Wohnort: Bayern, München
Zum letzten BeitragZum nächsten BeitragZum vorigen BeitragZum erstem Beitrag  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...

Eine Notiz zu diese Forumbeitrag schreiben Notiz   Profil  Quote  Link auf diesen Beitrag Link
Spock
Senior Letzter Besuch: in der letzten Woche
Dabei seit: 25.04.2002
Mitteilungen: 8199
Wohnort: Schi'Kahr/Vulkan
Zum letzten BeitragZum nächsten BeitragZum vorigen BeitragZum erstem Beitrag  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

Eine Notiz zu diese Forumbeitrag schreiben Notiz   Profil  Quote  Link auf diesen Beitrag Link
civilengineer
Wenig Aktiv Letzter Besuch: vor mehr als 3 Monaten
Dabei seit: 25.05.2011
Mitteilungen: 1762
Wohnort: Bayern, München
Zum letzten BeitragZum nächsten BeitragZum vorigen BeitragZum erstem Beitrag  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 <math> (\kappa) </math> und den Eigenvektoren gelöst. Ich weiss mittlerweile, dass man für die Verschiebungen den Ansatz <math> u(x) = A \cdot e^{i \kappa x} </math> 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 <math> \omega </math>) 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 <math> e^{i \kappa x} </math> 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 <math> (\kappa) </math> komplexe Zahlen sind. Die Verschiebung somit auch?!?! ... Die Verschiebung sollte aber eine reelle Grösse sein ... ?

Viele Grüße

CE

Eine Notiz zu diese Forumbeitrag schreiben Notiz   Profil  Quote  Link auf diesen Beitrag Link
civilengineer
Wenig Aktiv Letzter Besuch: vor mehr als 3 Monaten
Dabei seit: 25.05.2011
Mitteilungen: 1762
Wohnort: Bayern, München
Zum letzten BeitragZum nächsten BeitragZum vorigen BeitragZum erstem Beitrag  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? ...

Eine Notiz zu diese Forumbeitrag schreiben Notiz   Profil  Quote  Link auf diesen Beitrag Link
civilengineer
Wenig Aktiv Letzter Besuch: vor mehr als 3 Monaten
Dabei seit: 25.05.2011
Mitteilungen: 1762
Wohnort: Bayern, München
Zum letzten BeitragZum nächsten BeitragZum vorigen BeitragZum erstem Beitrag  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?!...

Eine Notiz zu diese Forumbeitrag schreiben Notiz   Profil  Quote  Link auf diesen Beitrag Link
civilengineer
Wenig Aktiv Letzter Besuch: vor mehr als 3 Monaten
Dabei seit: 25.05.2011
Mitteilungen: 1762
Wohnort: Bayern, München
Zum letzten BeitragZum nächsten BeitragZum vorigen BeitragZum erstem Beitrag  Beitrag No.6, vom Themenstarter, eingetragen 2017-06-09


Alle machen Informatik jetzt, niemand will Physik machen

Eine Notiz zu diese Forumbeitrag schreiben Notiz   Profil  Quote  Link auf diesen Beitrag Link
civilengineer
Wenig Aktiv Letzter Besuch: vor mehr als 3 Monaten
Dabei seit: 25.05.2011
Mitteilungen: 1762
Wohnort: Bayern, München
Zum letzten BeitragZum nächsten BeitragZum vorigen BeitragZum erstem Beitrag  Beitrag No.7, vom Themenstarter, eingetragen 2017-06-10


Jemand eine Idee?

Eine Notiz zu diese Forumbeitrag schreiben Notiz   Profil  Quote  Link auf diesen Beitrag Link
civilengineer
Wenig Aktiv Letzter Besuch: vor mehr als 3 Monaten
Dabei seit: 25.05.2011
Mitteilungen: 1762
Wohnort: Bayern, München
Zum letzten BeitragZum nächsten BeitragZum vorigen BeitragZum erstem Beitrag  Beitrag No.8, vom Themenstarter, eingetragen 2017-06-10


push

Eine Notiz zu diese Forumbeitrag schreiben Notiz   Profil  Quote  Link auf diesen Beitrag Link
civilengineer
Wenig Aktiv Letzter Besuch: vor mehr als 3 Monaten
Dabei seit: 25.05.2011
Mitteilungen: 1762
Wohnort: Bayern, München
Zum letzten BeitragZum nächsten BeitragZum vorigen BeitragZum erstem Beitrag  Beitrag No.9, vom Themenstarter, eingetragen 2017-06-13


Push

Eine Notiz zu diese Forumbeitrag schreiben Notiz   Profil  Quote  Link auf diesen Beitrag Link
civilengineer
Wenig Aktiv Letzter Besuch: vor mehr als 3 Monaten
Dabei seit: 25.05.2011
Mitteilungen: 1762
Wohnort: Bayern, München
Zum letzten BeitragZum nächsten BeitragZum vorigen BeitragZum erstem Beitrag  Beitrag No.10, vom Themenstarter, eingetragen 2017-06-13


push

Eine Notiz zu diese Forumbeitrag schreiben Notiz   Profil  Quote  Link auf diesen Beitrag Link
Spock
Senior Letzter Besuch: in der letzten Woche
Dabei seit: 25.04.2002
Mitteilungen: 8199
Wohnort: Schi'Kahr/Vulkan
Zum letzten BeitragZum nächsten BeitragZum vorigen BeitragZum erstem Beitrag  Beitrag No.11, eingetragen 2017-06-14


Hallo!
fed-Code einblenden
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


Eine Notiz zu diese Forumbeitrag schreiben Notiz   Profil  Quote  Link auf diesen Beitrag Link
civilengineer
Wenig Aktiv Letzter Besuch: vor mehr als 3 Monaten
Dabei seit: 25.05.2011
Mitteilungen: 1762
Wohnort: Bayern, München
Zum letzten BeitragZum nächsten BeitragZum vorigen BeitragZum erstem Beitrag  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

Eine Notiz zu diese Forumbeitrag schreiben Notiz   Profil  Quote  Link auf diesen Beitrag Link
civilengineer
Wenig Aktiv Letzter Besuch: vor mehr als 3 Monaten
Dabei seit: 25.05.2011
Mitteilungen: 1762
Wohnort: Bayern, München
Zum letzten BeitragZum nächsten BeitragZum vorigen BeitragZum erstem Beitrag  Beitrag No.13, vom Themenstarter, eingetragen 2017-06-19


In einem Paper mit dem folgenden Namen:

'A waveguide finite element and boundary element approach to calculating the sound radiated by railway and tram rails' vom 27.10.2008 habe ich genau diese Gleichung gefunden (wenn man F = 0 setzt):



fed-Code einblenden
...

Viele Grüße

CE

Eine Notiz zu diese Forumbeitrag schreiben Notiz   Profil  Quote  Link auf diesen Beitrag Link
civilengineer
Wenig Aktiv Letzter Besuch: vor mehr als 3 Monaten
Dabei seit: 25.05.2011
Mitteilungen: 1762
Wohnort: Bayern, München
Zum letzten BeitragZum nächsten BeitragZum vorigen BeitragZum erstem Beitrag  Beitrag No.14, vom Themenstarter, eingetragen 2017-06-19


Push

Eine Notiz zu diese Forumbeitrag schreiben Notiz   Profil  Quote  Link auf diesen Beitrag Link
Spock
Senior Letzter Besuch: in der letzten Woche
Dabei seit: 25.04.2002
Mitteilungen: 8199
Wohnort: Schi'Kahr/Vulkan
Zum letzten BeitragZum nächsten BeitragZum vorigen BeitragZum erstem Beitrag  Beitrag No.15, eingetragen 2017-06-20


Hallo!

2017-06-19 16:19 - civilengineer in Beitrag No. 13 schreibt:
...
fed-Code einblenden
...

Ja!

Gruß
Juergen

Eine Notiz zu diese Forumbeitrag schreiben Notiz   Profil  Quote  Link auf diesen Beitrag Link
civilengineer
Wenig Aktiv Letzter Besuch: vor mehr als 3 Monaten
Dabei seit: 25.05.2011
Mitteilungen: 1762
Wohnort: Bayern, München
Zum letzten BeitragZum nächsten BeitragZum vorigen BeitragZum erstem Beitrag  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 <math> \kappa_{i} </math> und der zugehörige Eigenvektor <math> U^{~}_{i}\left(\kappa_{i} \right) </math> sind immer konstant und größtenteils komplex...

Also wenn der Eigenvektor <math> U^{~}_{i}\left(\kappa_{i} \right) </math> konstant ist, dann kann ich ihn aus der Integration herausnehmen: <math> U^{~}_{i}\left(x \right) = U^{~}_{i}\left(\kappa_{i} \right) \cdot \frac{1}{2 \cdot \pi} \cdot \int_{-\infty}^{\infty} e^{-i \cdot \kappa_{i} \cdot x} d \kappa </math>. Korrekt?

Dann muss ich nur die Exponentialfunktion von <math> -\infty </math> bis <math> +\infty </math> integrieren: <math> U^{~}_{i}\left(x \right) = U^{~}_{i}\left(\kappa_{i}\right) \cdot \frac{1}{2 \cdot \pi} \cdot \frac{1}{-i \cdot x} \cdot \left(e^{-i \cdot \kappa_{i} \cdot x} \right) \limits_{- \infty}^{\infty} </math> Korrekt? ... Ich weiss nicht, wie das weitergehen soll ... cauchy residue theorem?!?! ...

Die Integration über <math> \kappa </math> von <math> -\infty </math> bis <math> +\infty </math> 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 <math> \kappa_{i} </math> und Eigenvektoren <math> U^{~}_{i}\left(\kappa_{i} \right) </math> 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

Eine Notiz zu diese Forumbeitrag schreiben Notiz   Profil  Quote  Link auf diesen Beitrag Link
civilengineer
Wenig Aktiv Letzter Besuch: vor mehr als 3 Monaten
Dabei seit: 25.05.2011
Mitteilungen: 1762
Wohnort: Bayern, München
Zum letzten BeitragZum nächsten BeitragZum vorigen BeitragZum erstem Beitrag  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

Eine Notiz zu diese Forumbeitrag schreiben Notiz   Profil  Quote  Link auf diesen Beitrag Link
civilengineer
Wenig Aktiv Letzter Besuch: vor mehr als 3 Monaten
Dabei seit: 25.05.2011
Mitteilungen: 1762
Wohnort: Bayern, München
Zum letzten BeitragZum nächsten BeitragZum vorigen BeitragZum erstem Beitrag  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

Eine Notiz zu diese Forumbeitrag schreiben Notiz   Profil  Quote  Link auf diesen Beitrag Link
civilengineer
Wenig Aktiv Letzter Besuch: vor mehr als 3 Monaten
Dabei seit: 25.05.2011
Mitteilungen: 1762
Wohnort: Bayern, München
Zum letzten BeitragZum nächsten BeitragZum vorigen BeitragZum erstem Beitrag  Beitrag No.19, vom Themenstarter, eingetragen 2017-06-22


Push

Eine Notiz zu diese Forumbeitrag schreiben Notiz   Profil  Quote  Link auf diesen Beitrag Link
civilengineer
Wenig Aktiv Letzter Besuch: vor mehr als 3 Monaten
Dabei seit: 25.05.2011
Mitteilungen: 1762
Wohnort: Bayern, München
Zum letzten BeitragZum nächsten BeitragZum vorigen BeitragZum erstem Beitrag  Beitrag No.20, vom Themenstarter, eingetragen 2017-06-25


Push

Eine Notiz zu diese Forumbeitrag schreiben Notiz   Profil  Quote  Link auf diesen Beitrag Link
civilengineer
Wenig Aktiv Letzter Besuch: vor mehr als 3 Monaten
Dabei seit: 25.05.2011
Mitteilungen: 1762
Wohnort: Bayern, München
Zum letzten BeitragZum nächsten BeitragZum vorigen BeitragZum erstem Beitrag  Beitrag No.21, vom Themenstarter, eingetragen 2017-07-05


Push

Eine Notiz zu diese Forumbeitrag schreiben Notiz   Profil  Quote  Link auf diesen Beitrag Link
civilengineer hat die Antworten auf ihre/seine Frage gesehen.
Neues Thema [Neues Thema] Antworten [Antworten]    Druckversion [Druckversion]

 


Wechsel in ein anderes Forum:
 Suchen    
 
All logos and trademarks in this site are property of their respective owner. The comments are property of their posters, all the rest © 2001-2021 by Matroids Matheplanet
This web site was originally made with PHP-Nuke, a former web portal system written in PHP that seems no longer to be maintained nor supported. PHP-Nuke is Free Software released under the GNU/GPL license.
Ich distanziere mich von rechtswidrigen oder anstößigen Inhalten, die sich trotz aufmerksamer Prüfung hinter hier verwendeten Links verbergen mögen.
Lesen Sie die Nutzungsbedingungen, die Distanzierung, die Datenschutzerklärung und das Impressum.
[Seitenanfang]