Autor
|
Thema: Abaqus - lokales Koordinatensystem (1583 / mal gelesen)
|
VR46M1 Mitglied M.Sc.
Beiträge: 5 Registriert: 16.05.2017
|
erstellt am: 16. Mai. 2017 16:13 <-- editieren / zitieren --> Unities abgeben:
Hallo zusammen, ich arbeite gerade an einem Materialmodell deren Verhalten transversal isotrop ist. jetzt ist die Faserorientierung jedoch nicht gleichbleibend sondern verhält sich von Rand zu Rand nach der Laplace-Gleichung (Interpolation). Jetzt war meine Idee erstmal für den vereinfachten Fall (Würfel) eine thermische Simulation durchzuführen(steady-state)und mir dann die Temperaturen ausgeben zu lassen. Diese Verhalten sich ja dementsprechend nach der Wärmeleitungsgleichung in stationären Fall und ohne Einwirkung äußerer Einflüsse (div=0) gleich der Laplace-Gleichung. Deshalb habe ich nun einfach am Würfel auf beiden Seiten die jeweilige Temperatur angelegt und mir dann die Temperaturen an den Nodes (NT11) in ein File schreiben lassen (ausgelesen aus dem odb-file). Nun möchte ich diese in die subroutine Orient einbinden. Jedoch habe ich ja die Knoten-Werte der Temperatur ausgelesen und ich weiß aktuell nicht genau auf was die subroutine genau ein lokales Koordinatensystem erzeugt(Im Manual steht Materialpoint...)... Geschieht das elementweiße, sprich ich muss die jeweiligen Knoten-Werte des jeweiligen Elementes(linear) irgendwie mitteln um dann dem KS den richtigen (gemittelten) Winkel für die Rotationsmatrix T zu übergeben? Ich bin für Vorschläge dankbar. Ps. Später möchte ich die Laplace-Gleichung für die 6 notwendigen Richtungen der Rotationsmatrix berechnen und somit eine vollständige Interpolation für die ersten beiden Richtungen durchführen. (Die dritte Raumrichtung wird ja laut Manual automatisch orthogonalisiert). Beste Grüße VR46M1 Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP |
Mustaine Ehrenmitglied V.I.P. h.c.
Beiträge: 3554 Registriert: 04.08.2005 Abaqus
|
erstellt am: 18. Mai. 2017 11:40 <-- editieren / zitieren --> Unities abgeben: Nur für VR46M1
Ich habe die Routine noch nie benutzt, aber ich denke das Manual ist recht eindeutig. Es wird jeder Material Point von jedem Element durch die Routine geschickt und bekommt seine Orientierung. Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP |
VR46M1 Mitglied M.Sc.
Beiträge: 5 Registriert: 16.05.2017
|
erstellt am: 18. Mai. 2017 13:34 <-- editieren / zitieren --> Unities abgeben:
|
Mustaine Ehrenmitglied V.I.P. h.c.
Beiträge: 3554 Registriert: 04.08.2005 Abaqus
|
erstellt am: 18. Mai. 2017 13:55 <-- editieren / zitieren --> Unities abgeben: Nur für VR46M1
|
VR46M1 Mitglied M.Sc.
Beiträge: 5 Registriert: 16.05.2017
|
erstellt am: 22. Mai. 2017 14:03 <-- editieren / zitieren --> Unities abgeben:
Ok. Dadurch entsteht natürlich das Problem, dass ich die berechneten Temperaturen nicht einfach lokal auf das Koordinatensystem aufbringen kann (Temperatur an Knoten <-> lokales Koordinatensystem an Integrationspunkten). Gibt es hierfür ein geeignetes, Verfahren um diese Lücke zu schließen? Danke! Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP |
Mustaine Ehrenmitglied V.I.P. h.c.
Beiträge: 3554 Registriert: 04.08.2005 Abaqus
|
erstellt am: 22. Mai. 2017 14:54 <-- editieren / zitieren --> Unities abgeben: Nur für VR46M1
|
VR46M1 Mitglied M.Sc.
Beiträge: 5 Registriert: 16.05.2017
|
erstellt am: 29. Mai. 2017 13:30 <-- editieren / zitieren --> Unities abgeben:
Hallo nochmal, bei der Definition des lokalen Koordinatensystems bin ich auf ein weiteres Problem gestoßen. Ich möchte lokal am Knoten den Normalenvektor bezogen auf die angrenzenden Flächen berechnen. Dazu würde ich gerne alle "benachbarten" Knoten herausfiltern (Dann das jeweilige Kreuzprodukt bilden und über alle Normalvektoren mitteln). Ich habe bereits ein Nodeset erstellt, das alle relevanten Knoten enthält und mir dazu die jeweiligen Koordinaten und Nodenumber ausgeben lassen. Ich müsste jetzt nur noch die Verbindung der einzelnen Knoten untereinander herstellen (Bsp. Node1={2,12,14,53,64}). Kann mir hier jemand behilflich sein? Ich bin wirklich dankbar, falls Ihr zur Lösung der Problematik beitragen könntet. Beste Grüße und vielen Dank! Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP |
Mustaine Ehrenmitglied V.I.P. h.c.
Beiträge: 3554 Registriert: 04.08.2005 Abaqus
|
erstellt am: 29. Mai. 2017 23:13 <-- editieren / zitieren --> Unities abgeben: Nur für VR46M1
Du müsstest ermitteln, welche Elemente deinen Knoten verwenden, von den Elementen dann alle Knoten ermitteln und dann die verwenden, die eine Schnittmenge mit deinem Knotenset bilden. Nicht einfach... Ich glaube bei Schalen werden die Knotennormalen in die .dat geschrieben, wenn man bei *Preprint eine extra Option aktiviert (evtl. Model=yes). Falls du Solids verwendest, kannst du ggf. eine dummy-Skin aus Schalen auf die Oberfläche legen.
Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP |