momentan versuche ich eine nichtlineare Berechnung über die Konsole zu starten. Leider bricht mit Berechnung immer mit einem Fehler ab, auch wenn sie im linearen Fall funktionierte. Die Suche nach einem Fehler hatte leider keinen Erfolg.
Die Verschiebungen scheinen alle 0 zu sein, genau so wie die Spannungen. Können sie mir bitte weiterhelfen? Ich bedanke mich bereits im Vorraus.
nach einigen Tests und Vergleichen von der Z88Aurora-Berechnung und der Berechnung des Solvers z88nl.exe mit den selbsterstellten Dateien ist mir ein Unterschied aufgefallen. Die z88mat.txt hat bei den selbsterstellten Dateien keinen 4. Wert, der für das Materialgesetzt stehen soll. Leider habe ich keine weiteren Angaben zu dieser Zahl bzw. zu dem Materialgesetz gefunden.
Welchen Wert muss die 4. Zahl in der z88mat.txt bei einer materiellen-nichtlinearen Berechnung annehmen?
vielen Dank für die Erklärung der Flags.
Die Vermischung wurde automatisch von Aurora erstellt. Dabei sollten alle Elemente eine Materiale Nichtliearität (Von Mises) aufweisen.
In Ihrer z88mat.txt sind nicht nur verschiedene Materialgesetze (1 bzw. sondern auch verschiedene Materialien (nl_1 bzw. nl_53) hinterlegt. Anscheinend haben Sie verschiedene Materialsets definiert und dafür dann verschiedene Gesetze ausgewählt. Es wird nicht abgeprüft, ob eine Vermischung der Materialgesetze eingestellt wurde. Die Vermischung von Materialgesetzen ist allerdings nicht rechenbar, verschiedene Materialen stellen sonst kein Problem dar. Falls sowieso das gesamte Bauteil aus einem Material sein soll, löschen Sie die aktuellen Sets und verwenden eines für sämtliche Elemente.
die Materialgesetzte sind nun überall auf "8" festgelegt. Die zwei verschiedenen Materialien sind erwünscht und notwendig. Nun kommt die Berechnung bei irgendeinem Punkt zum Stehen und bricht ab.
Der Inhalt der z88nl.log ist wie folgt:
### Fehler bei der Berechnung. Programmabbruch.
### Fehler bei der Berechnung von Element 102395
### Integrationspunkt: 1
### Unerlaubter Deformationszustand
Verschiebungen der Knoten dieses Elementes:
Knoten x-Verschiebung y-Verschiebung z-Verschiebung
117053 +2.6692586E-02 -5.2785424E-02 +2.1540567E-02
117052 +2.8062662E-02 -5.8304219E-02 +2.7528670E-02
117057 +1.0324979E-01 -7.0695523E-02 -3.4820742E-02
117058 +2.6916263E-02 -5.4686007E-02 +2.0715609E-02
116923 +2.7132667E-02 -5.4110150E-02 +2.1183406E-02
116922 +5.0884966E-02 -6.5924805E-02 +2.1233012E-02
116927 +1.3013712E-01 -9.2389544E-02 -4.6984453E-02
116928 +2.8263415E-02 -5.7127601E-02 +2.0042459E-02
Technische Dehnungen epsilon (gesamt) dieses Elements:
IP eps11 eps22 eps33 eps12 eps23 eps31
1 +3.2682951E-01 -9.3970072E-02 -2.2162035E-01 +1.6914101E-01 -2.6480572E-01 +4.2846280E-02
Plastische Vergleichsdehnung:
IP divpl
1 +3.7718788E-01
Plastisches Dehnungsinkrement:
IP delta_gam
1 +1.9864902E-03
zum einen könnte das Problem an der Integrationsordnung liegen. Was steht denn in Ihrer z88int.txt?
Zum anderen könnte eine zu starke Verformung des Netzes die Ursache sein. Könnten Sie mal ein Bild des Netzes hochladen?
Eine Darstellung der Verformung ist leider nicht möglich, da die Berechnung bis jetzt nicht erfolgreich war. Das Element scheint direkt an der Kraftrandbedingung zu liegen, wie es die Bilder Zeigen.
die z88int.txt ist so in Ordnung und rechenbar. Sollte die zu starke Verformung des Elements Schuld an dem Problem sein, könnten Sie mal die Kraft reduzieren und einen neuen Berechnungsversuch starten.
Wahrscheinlicher ist aber, dass es immer noch an der Materialzuweisung liegt. Mir ist aufgefallen, dass Sie 8100 Material-Sets nutzen. Ist das richtig? Falls ja wird das Problemelement, dasjenige sein, bei dem die Materialgesetze "ausgegangen" sind. Im NL ist nur eine begrenzte Anzahl (<8100) an Materialgesetzen erlaubt. Durch eine ungünstige Elementnummerierung kann es aber sein, dass diese Anzahl überschritten wird. Versuchen Sie doch mal ein Material zu verwenden und erneut zu rechnen. Alternativ können Sie uns das Projekt auch schicken und wir schauen uns das Ganze mal detailliert an.
leider sind die zwei verschiedenen Materialien wichtig. Die hohe Anzahl der Materialgesetzte in z88mat.txt ist der, wie du bereits erwähnt hast, nicht auf geringe Anzahl optimierte Nummerierung verschuldet. Nach der Verwendung nur eines Materials kommt es aber leider auch zu keinem Ergebnis.
wir haben die zugesendeten Dateien durchforstet und haben das Problem identifizieren können:
Und zwar lag der Fehler in der mat_NL_E295.txt bzw. der mat_NL_CuSn12Ni.txt.
Es muss die Anzahl der einzelnen Werte angegeben werden, nicht die Anzahl der Wertepaare.
vielen Dank für das Finden meines Problems - das hat mich schon lange verfolgt
Leider gab es nicht viel Dokumentation zu dem manuellen Erstellen von Materiealdateiel, vorallem nicht für den nichtlinearen Fall. Desweiteren ist mir aufgefallen, dass wenn man in z88 Aurora für alle Elemente ein Material ausgewählt, und dann in der Spalte für Materialgesetz das erwünschte auswählt, es nicht für alle Elemente übernommen wird, sondern nur für eine Zeile (siehe Beitrag).
Eine Frage habe ich noch: Stellt die große Anzahl der Materialgesetze in z88mat.txt (>24000) ein Problem dar? Es werden in Zukunft wahrscheinlich noch größere Geometrien mit mehr Materialgesetzen folgen.
den von Ihnen entdeckten Fehler konnten wir nachvollziehen. Ich habe das an unsere Entwickler weitergeleitet. Danke für den Hinweis.
Die Anzahl der Materialgesetzte ist momentan auf 10000 beschränkt. Wir könnten allerdings einen Solver mit mehr Möglichkeiten zur Verfügung stellen. Können Sie abschätzen wie viele Materialgesetze maximal bei Ihnen auftreten?
bei der aktuellen, erfolgreichen Rechnung sind es 56700 Materialgesetzte, das liegt vor allem daran, dass die Geometrie zweigeteilt ist. Die Nabe besteht aus Stahl und das Schneckenrad aus Bronze (siehe Abbildung).
Leider kommt es bei einer anderen, aber ähnlichen Geometrie wieder zu dem bekannten Fehler :
### Fehler bei der Berechnung. Programmabbruch.
### Fehler bei der Berechnung von Element 305525
### Integrationspunkt: 1
### Unerlaubter Deformationszustand
Verschiebungen der Knoten dieses Elementes:
Knoten x-Verschiebung y-Verschiebung z-Verschiebung
338982 +1.3525834E-01 -9.8597480E-02 -1.5686811E-02
338981 +2.4137926E-01 -1.5837802E-01 -4.2598553E-02
338986 +7.7653981E-01 -1.1964611E-01 +3.2032337E-02
338987 +1.4341233E-01 -1.5827713E-01 -3.0039969E-02
338852 +1.6197930E-01 -1.0141423E-01 -1.3529474E-02
338851 +2.0551211E-01 -9.7027283E-02 -2.3893277E-02
338856 +1.4950909E-01 -5.8756184E-02 -1.7050079E-02
338857 +1.5680287E-01 -7.8816289E-02 -1.5426050E-02
nach erfolgreicher Durchsicht der Daten konnten wir den Fehler finden.
In der Materialdatei mat_NL_E295 sind die Spannungswerte nur bis zu einer Dehnung von 1.0 angegeben. Dieser Wert wird allerdings überschritten. Sie müssten entsprechend die Wertepaare erweitern. Dann sollte alles wieder funktionieren