---

hallo gemeinde,

ich bin relativ neu und arbeite mich in Abaqus noch ein. Vielleicht kann mir jemand helfen, ich werde einfach nicht schlau aus der Steiffigkeitsmatrix als *.mtx (*MATRIX GENERATE, STIFFNESS)

Ich habe ein paar einfache Beispiele erstellt mit je einem Element:

• Einfacher Zugstab mit 2 Knoten (Truss Element):  1------2
  steifigkeitsmatrix:
  

  Code:

  ---

  1,1,1,1,1.000000000000000e+036
  2,1,1,1,-2.100000000000000e+004
  1,2,1,2,1.000000000000000e+036
  2,1,2,1,2.100000000000000e+004
  2,2,2,2,0.000000000000000e+000
  

  ---

• 2D Platte mit 4 Knoten (2D Planar Shell)
  steifigkeitsmatrix:
  

  Code:

  ---

  1,1,1,1,1.000000000000000e+036
  1,2,1,1,-3.750000000000000e+004
  2,1,1,1,3.744951923076924e+004
  2,2,1,1,-2.884615384615387e+003
  3,1,1,1,-7.783413461538462e+004
  3,2,1,1,3.750000000000000e+004
  4,1,1,1,-3.755048076923077e+004
  4,2,1,1,2.884615384615387e+003
  ....usw.  36 zeilen lang
  

  ---

• 3D Würfel mit 8 Knoten
  die steifigkeitsmatrix hat 300 zeilen wobei der 2. und 4. wert zwischen 1,2 und 3 wechselt.
  

Um mal beim ersten Beispiel zu bleiben. In jeder Zeile habe ich 5 attributes. Welche Informationen sind in den Werten? Das einzige was ich ansatzweise darüber in der Documentation gefunden habe war unter 5.1.2 und das hilft mir nicht wirklich weiter weil die globale Steifigkeitsmatrix nicht erwähnt wird.

Würde mich um jeden Tipp freuen

Gruß

Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP

---

Ich denke das ist gar nicht so schwer. Die Hilfe verrät ja eigentlich alles.
Die Steifigkeitsmatrix eines Stabes mit nur einem DOF pro Knoten ist ja bekannt.
K=EA/l*[1 -1,-1 1]
Deine erste Matrix hat nun 5 Zeilen und pro Zeile 5 Informationen.
Zunächst sei geklärt warum 5 Zeilen:
ABAQUS nutzt offenbar 2 DOF pro Knoten, damit ist deine Matrix von der Dimension 4x4. Alle Einträge müssen aber nicht abgespeichert werden, da Steifigkeitsmatrizen (meist) symmetrisch sind. Demzufolge wird nur die obere Dreiecksmatrix (alle Einträge die gleich Null sind und nicht auf der Hauptdiagonalen liegen werden NICHT gespeichert) abgespeichert.
Nun zu den 5 Informationen pro Zeile:
Das steht eigentlioch sehr gut in der Hilfe. An deinem Knoten 1 liegt offenbar eine feste Einspannung vor, daher Steifigkeit 1e36. Der negative Eintrag ist die Koppelsteifigkeit.

Hoffe es ist nun etwas klarer.

Beste Grüße

Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP

---

danke für deine Hife, leider ist mir noch nicht Klar wie ich aus den 4 ersten informationen herauslesen kann an welche stelle der Wert in der steifigkeitsmatrix steht.
die erste und dritte zahl steht für die knoten, die zweite und vierte für die freiheitsgrade?
Ich habe das erste Beispiel etwas geändert, damit es anschaulicher ist: der Stab ist jetzt diagonal ausgerichtet und die Dreiecksmatrix ist damit voll besetzt:

Code:

---

1,1,1,1,1.000000000000000e+036
1,2,1,1,1.484924240491750e+004
2,1,1,1,-1.484924240491750e+004
2,2,1,1,-1.484924240491750e+004
1,2,1,2,1.000000000000000e+036
2,1,1,2,-1.484924240491750e+004
2,2,1,2,-1.484924240491750e+004
2,1,2,1,1.484924240491750e+004
2,2,2,1,1.484924240491750e+004
2,2,2,2,1.484924240491750e+004

---

in dieser form brauche ich die ausgabe:
| a11 a12 a13 a14 |
| xxx a22 a23 a24 |
| xxx xxx a33 a34 |
| xxx xxx xxx a44 |

Gruß

Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP

---

Ok ich geb mir Mühe es dir näher zu bringen.
Dein Stab liegt nun also unter einem Winkel alpha in der Ebene. Somit ist die Matrix (wie du bereits sagtest) und von der Dimension 4x4 (in der Ebene). Somit müssen 10 Einträge abgespeichert werden. Di du ja auch vorliegen hast.
Du musst dir nun die Steifigkeitsmatrix unterteilt in Submatrizen vorstellen:
|K11 K12|
|K21 K22|
Diese Indizes werden durch den ersten und dritten Eintrag (pro Zeile) in deinem .mtx file definiert. Beachte das K11,K12,K21 und K22 ebenfalls Matrizen sind. Der zweite und vierte Eintrag im .mtx file geben nun die Position innerhalb der Submatrizen an.
Hoffe es ist verständlich.
Wenn nicht, nochmal fragen

Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP