DDBanner02

Zur Startseite

TM-aktuell
Aufgaben zur Kinematik und Kinetik und den Prinzipien der Mechanik

Kontrollrechnung mit Matlab-Femset

Gegeben: EI = 3000 Nm2 ; ρA = 3 kg/m l = 1 m  .

Der Finite-Elemente-”Baukasten” FEMSET dient eigentlich als Basis für den Aufbau individueller FEM-Programme. Weil in der MATLAB-Interface-Version aber ein Baustein “Eigenschwingungen” und das Element “Biegeträger” als Standard-Element vorhanden sind, können die Eigenfrequenzen für die Aufgabe 32-7 damit auf einfache Weise nachgerechnet werden. Zunächst wird die oben zu sehende sehr grobe Einteilung in nur 2 Elemente (3 Knoten) realisiert. Um das nebenstehende Matlab-Scripts laufen zu lassen, wird nur die DLL femeig.dll benötigt.

Das nebenstehende Matlab-Script zeigt den Aufbau der 4 Matrizen, die das Berechnungsmodell definieren:

  • xy enthält die Koordinaten der drei Knoten (bezogen auf den willkürlich in den Punkt 1 gelegten Ursprung),
  • km enthält die zu den beiden Elementen gehörenden Knotennummern (Koinzidenzmatrix),
  • ep enthält die (jeweils 2) Elementparameter der beiden Elemente: Biegesteifigkeit und Massebelegung,
  • kr enthält die Informationen über die verhinderten Verformungen (1 - Verformung verhindert, 0 - Verformung möglich), für jeden Knoten je ein Indikator für die Vertikalverschiebung und den Biegewinkel.

Mit dem Aufruf von femeig_m wird die Eigenkreisfrequenzberechnung gestartet, der erste Parameter gibt die Anzahl der gewünschten (kleinsten) Eigenkreisfrequenzen vor, ausgegeben wird in das Command Window die Eigenfrequenz (bei der äußerst groben Elementeinteilung ist es nicht sinnvoll, mehr als eine Frequenz berechnen zu lassen).

Das Ergebnis ist recht genau, um allerdings auch höhere Frequenzen zu berechnen, muss eine feinere Diskretisierung gewählt werden. Das folgende Matlab-Script realisiert eine beliebig feine (gleichmäßige) Elementeinteilung:

Die Anzahl der Elemente wird in Zeile 9 festgelegt. Schon bei einer Einteilung in nur 20 Elemente sind die ersten drei Eigenfrequenzen, die in das Command Window ausgegeben werden, praktisch exakt:

Zum Download verfügbar sind

Zur Übersicht der Aufgaben zur Kinematik und Kinetik und den Prinzipien der Mechanik

www.DankertDankert.de

TM1-Aufgaben TM2-Aufgaben
TM3-Aufgaben TM3-Aufgaben