Autor
|
Thema: lineare Tets (2575 mal gelesen)
|
vw-student Mitglied Berechnungsingenieur
Beiträge: 543 Registriert: 07.06.2005 IDEAS 12m4 ABAQUS/CAE 6.10-1 ADAMS MD R3 Matlab R2007a Mathcad 14
|
erstellt am: 12. Nov. 2008 13:35 <-- editieren / zitieren --> Unities abgeben:
Hallo. Ich versuche grad die Simulationsergebnisse einer dünnwandigen Struktur bezüglich der verwendeten Elemente zu vergleichen.Als erstes habe ich 8-Knoten Schalenelemente verwendet.Ergebnisse erscheinen plausibel.Danach habe ich 20-Knoten Hexa genommen.Bei geniügend feiner Diskretisierung weichen die Ergenisse von denen der Schalen ab.Nun soll ich aber versuchen diese Struktur mit linearen Tets zu vernetzen.Und genau hier liegt der Hund begraben.Egal wie fein ich vernetze,die Ergebnisse sind ehrlich gesagt ziemlich kacke.Und wenn ich zu fein werde bricht ABAQUS ab (insufficient memory). Hat evtl jemand Ideen, wie ich auch mit den linearen Tets zumindest mal in die Nähe der Schalenerbenisse komme?Und meint ihr,dass es hoffnungslos ist? Wäre für Kritik und Anregungen sehr dankbar. Grüße aus Berlin 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: 12. Nov. 2008 21:27 <-- editieren / zitieren --> Unities abgeben: Nur für vw-student
Die linearen Tets sind zu steif in ihrer Definition. Man muss schon recht viele nehmen um gute Ergebnisse zu bekommen. Deswegen werden sie auch kaum benutzt. An dem Problem kann man nichts ändern. Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP |
vw-student Mitglied Berechnungsingenieur
Beiträge: 543 Registriert: 07.06.2005 IDEAS 12m4 ABAQUS/CAE 6.10-1 ADAMS MD R3 Matlab R2007a Mathcad 14
|
erstellt am: 13. Nov. 2008 09:47 <-- editieren / zitieren --> Unities abgeben:
Recht viele ist gut Wenn ich die 8-Knoten-Shells nehme brauche ich 96 Elemente für ein konvergentes Netz.Bei den linearen Tets habe ich mittlerweile über 300.000 und die Ergebnisse sind immernoch deutlich zu steif! Naja...das ist ja auch eine Erkenntnis bzw ein Ergebnis.Jetzt muß ich wohl dran setzen und versuchen es mathematisch zu beweisen.Habe nur das dumme Gefühl, dass dies ein langer und harter Weg ist. Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP |
Mechanix Mitglied Mechanix
Beiträge: 29 Registriert: 20.12.2006
|
erstellt am: 14. Nov. 2008 09:18 <-- editieren / zitieren --> Unities abgeben: Nur für vw-student
hallo! dass die linearen tetraeder und linearen dreieckselemente viel zu steif und damit unbrauchbar sind, ist ein alter hut! in JEDEM besseren lehrbuch über fem wird dies angesprochen. stw. constant strain element lineare ansatzfunktionen approximieren das verschiebunsfeld linear dehnungsfeld bzw dehnungstensor enthält 1.partielle ableitungen des verschiebungsfelds. (lineare elasti-theo. ) wenn du eine lineare funktion ableitest erhältst du eine konstante (mathe gymmi-stoff) -> constant strain elements oder hier nachschauen: bathe, belytschko, zienkiewicz, hughes, skriptum, vorlesung, bibliothek, etc.... den chef fragen LG mechanix
Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP |
vw-student Mitglied Berechnungsingenieur
Beiträge: 543 Registriert: 07.06.2005 IDEAS 12m4 ABAQUS/CAE 6.10-1 ADAMS MD R3 Matlab R2007a Mathcad 14
|
erstellt am: 14. Nov. 2008 10:15 <-- editieren / zitieren --> Unities abgeben:
Ok.Bevor du mir indirekt vorwirfst,dass ich nichts von der Theorie verstehe,solltest du evtl erst nach den Hintergründen fragen.Es hat nämlich einen Sinn,dass ich diese Elemente verwenden möchte!! Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP |
Mechanix Mitglied Mechanix
Beiträge: 29 Registriert: 20.12.2006
|
erstellt am: 14. Nov. 2008 19:49 <-- editieren / zitieren --> Unities abgeben: Nur für vw-student
hallo! 1. ich werfe dir gar nichts vor. Zitat: Wäre für Kritik und Anregungen sehr dankbar
.... Zitat: Nun soll ich aber versuchen diese Struktur mit linearen Tets zu vernetzen
2. war das deine idee? 3. welche krümmungen hat dein bauteil, wie sieht es aus? etwa ein blechteil? belastungen? ... beschreib mal das problem nicht beleidigt sein.
LG Mechanix Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP |
vw-student Mitglied Berechnungsingenieur
Beiträge: 543 Registriert: 07.06.2005 IDEAS 12m4 ABAQUS/CAE 6.10-1 ADAMS MD R3 Matlab R2007a Mathcad 14
|
erstellt am: 18. Nov. 2008 10:34 <-- editieren / zitieren --> Unities abgeben:
Es handelt sich um eine ganz simple ebene Platte mit den Abmessungen 1m x 1,5m x 0,02m. Werkstoff ist Stahl.Als Belastung wirkt eine konstante Flächenlast q=50kPa. Und als Lagerung habe ich die beiden Kanten der Länge 1m gelenkig gelagert. Wie gesagt mit Shells oder Hex funktioniert alles wunderbar und eine Konvergenz wird rasch erreicht. Und die linearen Tets...naja...wie bereits erwähnt. Habe mich nun darauf beschränkt zu sagen, dass diese Elemente keineswegs geeignet sind und wir (meine Betreuer und ich)uns auf die Realisierung unseres Vorhabens mittels anderer Elemente konzentrieren müssen. Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP |
KubaG Mitglied
Beiträge: 165 Registriert: 18.11.2002
|
erstellt am: 19. Nov. 2008 08:51 <-- editieren / zitieren --> Unities abgeben: Nur für vw-student
Hallo vw-student, mit linearen Tetraedern solltest Du mit ca. 30 000 000 Elementen sehr brauchbare Ergebnisse hinsichtlich der Verschiebung erhalten. Deine Längen-Breiten-Höhen Verhältnisse sind für 3D-Rechnungen etwas ungünstig. Ich würde bei 2D-Rechnung bleiben. Wie hast Du die gelenkige Lagerung der Blechkanten realisiert ? Gruß KubaG Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP |
vw-student Mitglied Berechnungsingenieur
Beiträge: 543 Registriert: 07.06.2005 IDEAS 12m4 ABAQUS/CAE 6.10-1 ADAMS MD R3 Matlab R2007a Mathcad 14
|
erstellt am: 19. Nov. 2008 12:22 <-- editieren / zitieren --> Unities abgeben:
30 000 000 Elemente? Hmmm...mein Rechner gibt bereits bei 1 600 000 Elemente auf. ABAQUS bricht dann die Rechnung ab. Zu den Ranbedingungen: Da ich bei der Schale ohne Offset gearbeitet habe, habe ich das Solid in der Mitte partitioniert und die RB an den entsprechenden Mittelflächenkanten angebracht. Zur Frage 3D vs 2D: Mir ist schon klar, dass Schalen für solch dünnwandige Strukturen die sinnvollere effizientere Methode darstellen. Aber die Verwendung von Konti-Elementen hat seine Berechtigung, da wir hier am Institut bestimmte Dinge mit den linearen Tets vorhaben. Allerdings sieht es momentan so aus, als müssten wir auf die 10 Knoten Tets umsteigen, da eine solch immense Elementanzahl nicht hinnehmbar ist. Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP |
vw-student Mitglied Berechnungsingenieur
Beiträge: 543 Registriert: 07.06.2005 IDEAS 12m4 ABAQUS/CAE 6.10-1 ADAMS MD R3 Matlab R2007a Mathcad 14
|
erstellt am: 20. Nov. 2008 10:53 <-- editieren / zitieren --> Unities abgeben:
Ich habe jetzt mal 4 Elemente über die Dicke genommen und "in-plane" so vernetzt, dass das Seitenverhältnis der Elementkanten nicht größer als 2 ist. Nur leider bricht ABAQUS mir den Job mit folgender Meldung ab: Error in job linear_Tet_0005: The executable C:\ABAQUS\6.7-4\exec\pre.exe aborted with system error code 5. Please check the .dat, .msg, and .sta files for error messages if the files exist. If there are no error messages and you cannot resolve the problem, please run the command "abaqus job=support information=support" to report and save your system information. Use the same command to run Abaqus that you used when the problem occurred. Please contact your local Abaqus support office and send them the input file, the file support.log which you just created, the executable name, and the error code. Job linear_Tet_0005 aborted due to errors. Weiß evtl jemand Rat??? Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP |
vw-student Mitglied Berechnungsingenieur
Beiträge: 543 Registriert: 07.06.2005 IDEAS 12m4 ABAQUS/CAE 6.10-1 ADAMS MD R3 Matlab R2007a Mathcad 14
|
erstellt am: 21. Nov. 2008 15:15 <-- editieren / zitieren --> Unities abgeben:
@KubaG: Hast du das mal mit 30 000 000 Elementen durchgerechnet? Oder woher weißt du, dass ich mit dieser Elementanzahl mit brauchbaren Ergebnissen rechnen kann? Ich werd das gleich mal bei uns auf den Großrechner schmeißen. Mal schauen, was der dann sagt. Allerdings vernetzt ABAQUS schon seit 12 Minuten Mal schauen, ob dat heut noch was wird. Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP |
Mechanix Mitglied Mechanix
Beiträge: 29 Registriert: 20.12.2006
|
erstellt am: 23. Nov. 2008 18:18 <-- editieren / zitieren --> Unities abgeben: Nur für vw-student
hallo! darf ich mal zusammenfassen: Zitat: .... 30 000 000 Elemente
Zitat: .... Großrechner
Zitat: Es handelt sich um eine ganz simple ebene Platte mit den Abmessungen 1m x 1,5m x 0,02m. Werkstoff ist Stahl.Als Belastung wirkt eine konstante Flächenlast q=50kPa.
hoffe Du meinst das nicht ernst. Zitat: Aber die Verwendung von Konti-Elementen hat seine Berechtigung, da wir hier am Institut bestimmte Dinge mit den linearen Tets vorhaben.
willst du schalenartige strukturen untersuchen? möglichweise kannst du diese "bestimmten dinge" auch mit schalenelementen umsetzen. wenn's nicht standardmäßig funktioniert besteht die möglichkeit ein user-element zu programmieren. hast du das schon mal überlegt? p.s. "bestimmte dinge" (geheimnisvoll, ich tippe auf ein schädigungsmodell) LG Mechanix
Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP |
carsten1983 Mitglied Doktorand
Beiträge: 125 Registriert: 11.10.2007
|
erstellt am: 24. Nov. 2008 17:38 <-- editieren / zitieren --> Unities abgeben: Nur für vw-student
Mich würde ebenfalls der Hintergrund interessieren, warum es unbedingt Tetraeder sein sollen. Lässt sich das nicht umgehen? Bekommst du irgendwelchen Input eines anderen Programms nur als Tetraedernetz? Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP |
KubaG Mitglied
Beiträge: 165 Registriert: 18.11.2002
|
erstellt am: 15. Dez. 2008 14:28 <-- editieren / zitieren --> Unities abgeben: Nur für vw-student
Hallo VW-Student, das mit den 30 Mio. Elementen ist grob überschlagen und entspricht einer luxuriösen Lösung mit regulären Hexaedern (15Mio 1x1x2 HEX8 halbiert 30Mio Tet4). Mit etwas Mühe wirst Du noch etwas abspecken können, Du wirst aber m.E. im Mio.-Bereich hinsichtlich der Elementzahl bleiben, falls Du ein konvergentes Ergebnis anstrebst. Gruß KubaG [Diese Nachricht wurde von KubaG am 15. Dez. 2008 editiert.] Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP |