| | | Gut zu wissen: Hilfreiche Tipps und Tricks aus der Praxis prägnant, und auf den Punkt gebracht für Ansys |
Autor
|
Thema: NUMMRG -> LS-DYNA Konvertierung (2308 mal gelesen)
|
Deffensor Mitglied Student
Beiträge: 36 Registriert: 15.04.2016
|
erstellt am: 07. Mai. 2016 11:03 <-- editieren / zitieren --> Unities abgeben:
Hallo liebe Ansys Freunde, bei meiner Belegarbeit bin auf ein kleines Problem gestoßen. Ich habe einen Scherprüfkörper (siehe Anhang) in Ansys modelliert. Dabei habe ich 5 Teilvolumen erstellt, welche ich jeweils vernetzt habe. Das FE-Netz in Teilvolumen 1 und 2 habe ich nach der Vernetzung durch nummrg,all zusammengefügt. Volumen 3 und 4 habe ich auch selektiert und durch den Befehl nummrg zusammengeführt. Das Volumen 5 soll ein anderes Material als den anderen Volumen zugeordnet werden. Es ist quasi das Kohäsivelement zwischen den Schichten. Jetzt habe ich folgende Frage: Sollte ich nach dem ich das Volumen 5 erstellt habe, noch ein Mal nummrg,all durchführen? Dann könnte ich mir jedoch den nummrg Befehl bei den anderen Teilvolumen sparen, da ich ja sowieso alles verbinde. Sehe ich das richtig? An den Stirnkanten des Prüfkörpers soll eine Zugkraft angreifen, wodurch dementsprechend der Prüfkörper an der Kohäsivschicht abscheren wird. Demenstprechend müsste ich doch wenn ich die gleiche Elementgröße gewählt habe genauso viele Knoten an beiden Stirnkanten haben, oder etwa nicht? Nach der Konvertierung in das Programm LS-DYNA habe ich leider eine Diskripanz von 50 Knoten Vielen Dank im Voraus
Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP |
wosch Ehrenmitglied V.I.P. h.c. Elektrotechniker im Ruhestand
Beiträge: 2607 Registriert: 16.12.2004 Rechne zuerst ein Problem nach, für das Du eine analytische Lösung kennst.
|
erstellt am: 08. Mai. 2016 08:05 <-- editieren / zitieren --> Unities abgeben: Nur für Deffensor
|
Deffensor Mitglied Student
Beiträge: 36 Registriert: 15.04.2016
|
erstellt am: 08. Mai. 2016 10:15 <-- editieren / zitieren --> Unities abgeben:
Wahrscheinlich ja, doch dann entsteht das Problem, dass ich eine Parametrisierung des Winkels vornehme und die Vernetzung an der Kerbe nicht funktioniert. Mit dem Befehl Nummrg sieht das resultierende FE-Netz sehr gut aus. Das möchte ich schon beibehalten. Nummrg sagt ja aus, dass doppelte Knoten verschmolzen werden, was ja das Ziel ist. Nur bin ich mir nicht sicher, ob bei einem "kontinuierlichen" Bauteil die FE-Netze aller Volumen verbinden sollte. Wahrscheinlich ja, oder? Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP |
Rainer Schulze Ehrenmitglied V.I.P. h.c. Dipl.-Ing. im Ruhestand
Beiträge: 4419 Registriert: 24.09.2012
|
erstellt am: 08. Mai. 2016 11:17 <-- editieren / zitieren --> Unities abgeben: Nur für Deffensor
>>Nur bin ich mir nicht sicher, ob bei einem "kontinuierlichen" Bauteil die FE-Netze aller Volumen verbinden sollte Nach meinem Verständnis handelt es sich für die Simulation nicht um EIN Bauteil sondern um DREI. >>Das Volumen 5 soll ein anderes Material als den anderen Volumen zugeordnet werden. Denkst Du, es würde Dir gelingen, EINEM Bauteil mehrere Materialien zuzuweisen? >>dass ich eine Parametrisierung des Winkels vornehme Welches Winkels? Habe ich da in Deiner Skizze etwas übersehen? Bist Du ausdrücklich aufgefordert, mit APDL zu arbeiten? Mit der Workbench und dem Design Modeler wäre die Arbeit schnell erledigt. ------------------ Rainer Schulze Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP |
Deffensor Mitglied Student
Beiträge: 36 Registriert: 15.04.2016
|
erstellt am: 08. Mai. 2016 12:44 <-- editieren / zitieren --> Unities abgeben:
Hallo, ja eigentlich handelt es sich ja um drei Teile. Volumen 1+2 sind Bauteil 1. Volumen 3+4 sind Bauteil 2. Volumen 5 ist Bauteil 3. Da habe ich mich falsch ausgedrückt. Jetzt kann man unterschieldliche Materialien zuweisen Ich meine aber, dass Bauteil 1 und 2 durch Bauteil 3 miteinander verbunden sind. Dadurch kann ich ja quasi nummrg auf alle 5 Volumen anwenden, um die Verbindung der FE-Netze zu gewährleisten? Also nachdem ich jedes Bauteil einzeln vernetzt habe: nummrg,all anwenden müsste ja dann reichen, oder? Ja ich bin aufgefordert mit dem APDL-Skript zu arbeiten. Ich muss in einer Parameterstudie den Einfluss des Kerbwinkel auf die Scherfestigkeit untersuchen. Das habe ich nochmal als 2D Skizze angehängt. Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP |
Rainer Schulze Ehrenmitglied V.I.P. h.c. Dipl.-Ing. im Ruhestand
Beiträge: 4419 Registriert: 24.09.2012
|
erstellt am: 08. Mai. 2016 15:47 <-- editieren / zitieren --> Unities abgeben: Nur für Deffensor
>>Also nachdem ich jedes Bauteil einzeln vernetzt habe: nummrg,all anwenden müsste ja dann reichen, oder? Da sag ich doch glatt: Probieren geht über studieren. Das Ergebnis kannst Du in wenigen Minuten ohne die Hilfe des Forums sehen. ------------------ Rainer Schulze Da es mir an Übung in APDL mangelt zeige ich Dir anhand der Workbench, was geschehen wird, wenn die Netze auf beiden Seiten eines Kontaktes nicht deckungsgleich sind. [Diese Nachricht wurde von Rainer Schulze am 08. Mai. 2016 editiert.] Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP |
Deffensor Mitglied Student
Beiträge: 36 Registriert: 15.04.2016
|
erstellt am: 08. Mai. 2016 21:40 <-- editieren / zitieren --> Unities abgeben:
|
wosch Ehrenmitglied V.I.P. h.c. Elektrotechniker im Ruhestand
Beiträge: 2607 Registriert: 16.12.2004 Rechne zuerst ein Problem nach, für das Du eine analytische Lösung kennst.
|
erstellt am: 09. Mai. 2016 07:33 <-- editieren / zitieren --> Unities abgeben: Nur für Deffensor
Grundsätzlich: erst alle Volumen zusammenkleben und erst danach vernetzen! Wenn ich den Quelltext bekomme, sehe ich mir ihn an. ------------------ Viel Erfolg wünscht Wolfgang Schätzing [Diese Nachricht wurde von wosch am 09. Mai. 2016 editiert.] Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP |
Deffensor Mitglied Student
Beiträge: 36 Registriert: 15.04.2016
|
erstellt am: 10. Mai. 2016 08:44 <-- editieren / zitieren --> Unities abgeben:
Also ich füge meinen Code mal ein. Leider erhalte ich aber nicht den "Merge-Effekt" wie ihn Rainer Schulze per Screenshot gezeigt hat. (siehe Screenshot) An welcher Stelle sollte ich am Besten Vglue mit einfügen? Ich habe es probiert, jedoch kam immer die Fehlermeldung, dass die Volumen nicht vernetzt werden können. Das was ich ja in Ansys erreichen möchte, ist nur die Geometrie + FE-Netz, da die Rechnung dann im anderen Programm ablaufen wird. Die Knoten müssten dazu ja aber verschmelzt werden. !---------------------------------------- C*** globale Parametereingabe !----------------------------------------
! Eingabegrößen L=200 ! Probenlänge [mm] B=25 ! Probenbreite [mm] D2=4 ! Probendicke [mm] ! Normbedingt konstant AK=12.5 ! Abstand der Kerben [mm] LK=2 ! Länge der Kerben [mm] KF=0.45 ! Faktor der Kerbe [-] ! zwischen 0.45 ... 0.48 pi=4*atan(1) ! Pi ! Resultierende Größen D1=KF*D2 ! Probendicke - Kerbtiefe [mm] LangX=0.5*L+0.5*AK ! Lange Seite X [mm] KurzX=0.5*L-0.5*AK-LK ! Kurze Seite X [mm] ! Kerbwinkel Alpha=0 ! Kerbwinkel [°] BM=((Alpha*pi)/(180)) ! Winkel im Bogenmaß
! Elementgrößen elementgroesse1=2.0 !Elemtgröße Grob -> Außen elementgroesse2=1.0 !Elemtgröße Fein -> Kerbe elementgroesse3=0.25 !Elemtgröße Geometrie 3 -> Kohäsivelemente ! Elementanzahl elementanzahl1=25 ! Elementanzahl 1 -> Grob - Außen elementanzahl2=50 ! Elementanzahl 2 -> Fein - Kerbe elementanzahl3=20 ! Elementanzahl 3 -> Kohäsivelemente /prep7 ! Starten Preprocess
!----------------------------------------- C*** Zuweisungen für Konvertierung !-----------------------------------------
! Für Deckschichten - Real Constant = 50 TYPE,1 ! Auswahl des 1. Elementtyps SECNUM,1 ! Auswahl der 1. Section MAT,1 ! Auswahl des 1. Materials ! Für Kohäsivelemente - Real Constant = 100 TYPE,2 ! Auswahl des 2. Elementtyps SECNUM,2 ! Auswahl der 2. Section MAT,2 ! Auswahl des 2. Materials !----------------------------------------- C*** Geometrie 1 !----------------------------------------- ! Elementzuweisung et,1,solid186 ! Elementtypauswahl !NEU k,1,0,0,0 ! Keypoints 1-8 definieren Volumen 1 k,2,KurzX,0,0 k,3,KurzX,0,D2 k,4,0,0,D2 k,5,0,B,0 k,6,KurzX,B,0 k,7,KurzX,B,D2 k,8,0,B,D2 v,1,2,3,4,5,6,7,8 k,9,KurzX,0,D1 ! Keypoints 9-16 definieren Volumen 2 k,10,LangX,0,D1 k,11,LangX-D1*tan(BM),0,0 k,12,KurzX,0,0 k,13,KurzX,B,D1 k,14,LangX,B,D1 k,15,LangX-D1*tan(BM),B,0 k,16,KurzX,B,0 v,9,10,11,12,13,14,15,16 !nummrg,all vsel,s,volu,,1,2 ! Zuweisung für Konvertierung vatt,1,50 esize,elementgroesse1 vmesh,1 esize,elementgroesse2 vmesh,2 !----------------------------------------- C*** Geometrie 2 !----------------------------------------- ! NEU k,101,LangX+LK,0,0 ! Keypoints 101-108 definieren Volumen 3 k,102,L,0,0 k,103,L,0,D2 k,104,LangX+LK,0,D2 k,105,LangX+LK,B,0 k,106,L,B,0 k,107,L,B,D2 k,108,LangX+LK,B,D2 v,101,102,103,104,105,106,107,108 k,109,KurzX+LK,0,D2-D1 ! Keypoints 109-116 definieren Volumen 4 k,110,LangX+LK,0,D2-D1 k,111,LangX+LK,0,D2 k,112,KurzX+LK+D1*tan(BM),0,D2 k,113,KurzX+LK,B,D2-D1 k,114,LangX+LK,B,D2-D1 k,115,LangX+LK,B,D2 k,116,KurzX+LK+D1*tan(BM),B,D2 v,109,110,111,112,113,114,115,116
!vsel,s,volu,,3,4 !nummrg,all vsel,s,volu,,3,4 ! Zuweisung für Konvertierung vatt,1,50 esize,elementgroesse1 vmesh,3 esize,elementgroesse2 vmesh,4 !----------------------------------------- C*** Geometrie 3 (Klebverbindung) !-------------------------------------------- !Keypointdefinition k,201,KurzX+LK,0,D1 ! Diese Keypoints definieren die Matrixschicht zwischen den Faserlagen k,202,LangX,0,D1 k,203,LangX,0,D2-D1 k,204,KurzX+LK,0,D2-D1 k,205,KurzX+LK,B,D1 k,206,LangX,B,D1 k,207,LangX,B,D2-D1 k,208,KurzX+LK,B,D2-D1 v,201,202,203,204,205,206,207,208 ! Vernetzung vsel,s,volu,,5 ! Zuweisung für Konvertierung vatt,2,100 esize,elementgroesse3 vmesh,5 allsel nummrg,all !----------------------------- C*** element parts (FÜR KONVERTIERUNG IN LS-DYNA) !----------------------------- !* element part id of first part !* element type of first part: 1...beam !* 2...shell ! 3...solid !DPATH='/mnt/lft/lft_home/studenten/s5591929/Grosser_Beleg/LS-DYNA' DPATH='/w/Studenten/s5591929' DIRECTORY='/0°'
el_pzahl=0 el_pzahl=el_pzahl+1 esel,none ESEL,s,REAL,,50,,,1 CM,elementpart_%el_pzahl%,elem elementpart_id_%el_pzahl% = 50 elementtype_%el_pzahl% = 3 !2-Shell, 3-Volumen
el_pzahl=el_pzahl+1 esel,none ESEL,s,REAL,,100,,,1 CM,elementpart_%el_pzahl%,elem elementpart_id_%el_pzahl% = 100 elementtype_%el_pzahl% = 3 C*** CONVERTING TO LS-DYNA *** C*** element parts elementpart_number = el_pzahl !* number of element parts
C*** NODE SETS *** allsel !*** Auswählen alle Mitglieder des Modells nsel,s,LOC,x,0 ! Knotensets cm,nodeset_1,node nodeset_id_1 = 1000
allsel nsel,s,LOC,x,200 cm,nodeset_2,node nodeset_id_2 = 2000 allsel nsel,all NSEL,s,ext cm,nodeset_3,node nodeset_id_3 = 3000 nodeset_number = 3 C*** converter esel,all eplot
/SYP,/w/bin/erstDYNAModellVerz.csh,'%DPATH%','%DIRECTORY%' *use,/w/bin/convert_to_dyna.enc *MSG,NOTE,DPATH,DIRECTORY Pfad: %C; Verzeichnis: %C /eof
[Diese Nachricht wurde von Deffensor am 10. Mai. 2016 editiert.] [Diese Nachricht wurde von Deffensor am 10. Mai. 2016 editiert.] Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP |
Deffensor Mitglied Student
Beiträge: 36 Registriert: 15.04.2016
|
erstellt am: 10. Mai. 2016 09:37 <-- editieren / zitieren --> Unities abgeben:
Mit folgendem Code habe ich in Längsrichtung eine Verschmelzung. Aber ob das richtig ich ist weiß ich leider auch nicht. Ich bin für Eure Hilfe dankbar. ! Eingabegrößen L=200 ! Probenlänge [mm] B=25 ! Probenbreite [mm] D2=4 ! Probendicke [mm]
! Normbedingt konstant AK=12.5 ! Abstand der Kerben [mm] LK=2 ! Länge der Kerben [mm] KF=0.45 ! Faktor der Kerbe [-] ! zwischen 0.45 ... 0.48 pi=4*atan(1) ! Pi ! Resultierende Größen D1=KF*D2 ! Probendicke - Kerbtiefe [mm] LangX=0.5*L+0.5*AK ! Lange Seite X [mm] KurzX=0.5*L-0.5*AK-LK ! Kurze Seite X [mm] ! Kerbwinkel Alpha=0 ! Kerbwinkel [°] BM=((Alpha*pi)/(180)) ! Winkel im Bogenmaß
! Elementgrößen elementgroesse1=2.0 !Elemtgröße Grob -> Außen elementgroesse2=1.0 !Elemtgröße Fein -> Kerbe elementgroesse3=0.25 !Elemtgröße Geometrie 3 -> Kohäsivelemente ! Elementanzahl elementanzahl1=25 ! Elementanzahl 1 -> Grob - Außen elementanzahl2=50 ! Elementanzahl 2 -> Fein - Kerbe elementanzahl3=20 ! Elementanzahl 3 -> Kohäsivelemente /prep7 ! Starten Preprocess
!----------------------------------------- C*** Zuweisungen für Konvertierung !-----------------------------------------
! Für Deckschichten - Real Constant = 50 TYPE,1 ! Auswahl des 1. Elementtyps SECNUM,1 ! Auswahl der 1. Section MAT,1 ! Auswahl des 1. Materials ! Für Kohäsivelemente - Real Constant = 100 TYPE,2 ! Auswahl des 2. Elementtyps SECNUM,2 ! Auswahl der 2. Section MAT,2 ! Auswahl des 2. Materials !----------------------------------------- C*** Geometrie 1 !----------------------------------------- ! Elementzuweisung et,1,solid186 ! Elementtypauswahl !NEU k,1,0,0,0 ! Keypoints 1-8 definieren Volumen 1 k,2,KurzX,0,0 k,3,KurzX,0,D2 k,4,0,0,D2 k,5,0,B,0 k,6,KurzX,B,0 k,7,KurzX,B,D2 k,8,0,B,D2 v,1,2,3,4,5,6,7,8 k,9,KurzX,0,D1 ! Keypoints 9-16 definieren Volumen 2 k,10,LangX,0,D1 k,11,LangX-D1*tan(BM),0,0 k,12,KurzX,0,0 k,13,KurzX,B,D1 k,14,LangX,B,D1 k,15,LangX-D1*tan(BM),B,0 k,16,KurzX,B,0 v,9,10,11,12,13,14,15,16 nummrg,all vsel,s,volu,,1,2 ! Zuweisung für Konvertierung vatt,1,50 esize,elementgroesse1 vmesh,1 esize,elementgroesse2 vmesh,2 !----------------------------------------- C*** Geometrie 2 !----------------------------------------- ! NEU k,101,LangX+LK,0,0 ! Keypoints 101-108 definieren Volumen 3 k,102,L,0,0 k,103,L,0,D2 k,104,LangX+LK,0,D2 k,105,LangX+LK,B,0 k,106,L,B,0 k,107,L,B,D2 k,108,LangX+LK,B,D2 v,101,102,103,104,105,106,107,108 k,109,KurzX+LK,0,D2-D1 ! Keypoints 109-116 definieren Volumen 4 k,110,LangX+LK,0,D2-D1 k,111,LangX+LK,0,D2 k,112,KurzX+LK+D1*tan(BM),0,D2 k,113,KurzX+LK,B,D2-D1 k,114,LangX+LK,B,D2-D1 k,115,LangX+LK,B,D2 k,116,KurzX+LK+D1*tan(BM),B,D2 v,109,110,111,112,113,114,115,116 nummrg,all vsel,s,volu,,3,4 ! Zuweisung für Konvertierung vatt,1,50 esize,elementgroesse1 vmesh,3 esize,elementgroesse2 vmesh,4 !----------------------------------------- C*** Geometrie 3 (Klebverbindung) !--------------------------------------------
!Keypointdefinition k,201,KurzX+LK,0,D1 ! Diese Keypoints definieren die Matrixschicht zwischen den Faserlagen k,202,LangX,0,D1 k,203,LangX,0,D2-D1 k,204,KurzX+LK,0,D2-D1 k,205,KurzX+LK,B,D1 k,206,LangX,B,D1 k,207,LangX,B,D2-D1 k,208,KurzX+LK,B,D2-D1 v,201,202,203,204,205,206,207,208 ! Vernetzung nummrg,all vsel,s,volu,,5 ! Zuweisung für Konvertierung vatt,2,100 esize,elementgroesse3 vmesh,5 /eof
Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP |
wosch Ehrenmitglied V.I.P. h.c. Elektrotechniker im Ruhestand
Beiträge: 2607 Registriert: 16.12.2004 Rechne zuerst ein Problem nach, für das Du eine analytische Lösung kennst.
|
erstellt am: 10. Mai. 2016 16:39 <-- editieren / zitieren --> Unities abgeben: Nur für Deffensor
Hier im Anhang eine Vernetzungsmöglichkeit. Natürlich kann man das auch anders machen, insbesondere wenn es nur quaderförmige Elemente sein sollen. ------------------ Viel Erfolg wünscht Wolfgang Schätzing Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP |
Deffensor Mitglied Student
Beiträge: 36 Registriert: 15.04.2016
|
erstellt am: 10. Mai. 2016 19:40 <-- editieren / zitieren --> Unities abgeben:
Vielen Dank für die Lösungsmöglichkeit Das Ergebnis werde ich mir erst morgen an der Uni anschauen können. Deshalb möchte ich mal nachfragen, ob denn damit quaderförmige Elemente entstehen? Oder was könnte ich noch diesbezüglich verändern? Vielen Dank im Voraus [Diese Nachricht wurde von Deffensor am 10. Mai. 2016 editiert.] Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP |
wosch Ehrenmitglied V.I.P. h.c. Elektrotechniker im Ruhestand
Beiträge: 2607 Registriert: 16.12.2004 Rechne zuerst ein Problem nach, für das Du eine analytische Lösung kennst.
|
erstellt am: 11. Mai. 2016 07:31 <-- editieren / zitieren --> Unities abgeben: Nur für Deffensor
|
Deffensor Mitglied Student
Beiträge: 36 Registriert: 15.04.2016
|
erstellt am: 11. Mai. 2016 08:29 <-- editieren / zitieren --> Unities abgeben:
Mein Betreuer an der Universität meint ja. Aber das wird ja wahrscheinlich nur schwer zu realisieren sein, kann ich mir vorstellen. Aufgrund, dass die Größe meiner Volumenelemente unterschiedlich ist, wird es doch wahrscheinlich nicht möglich sein nur mit Quaderelementen zu arbeiten. Diese würden ja auch verzerrt sein. [Diese Nachricht wurde von Deffensor am 11. Mai. 2016 editiert.] Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP |
wosch Ehrenmitglied V.I.P. h.c. Elektrotechniker im Ruhestand
Beiträge: 2607 Registriert: 16.12.2004 Rechne zuerst ein Problem nach, für das Du eine analytische Lösung kennst.
|
erstellt am: 11. Mai. 2016 13:14 <-- editieren / zitieren --> Unities abgeben: Nur für Deffensor
|
Deffensor Mitglied Student
Beiträge: 36 Registriert: 15.04.2016
|
erstellt am: 11. Mai. 2016 13:54 <-- editieren / zitieren --> Unities abgeben:
Hallo Wosch, vielen Dank für Deine super Hilfe Das hat mir sehr geholfen Ich denke jedoch, dass aufgrund der sehr hohen Elementenanzahl bei der Quadervariante wohl Dein erster Vorschlag der geeignetere ist. Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP |
Deffensor Mitglied Student
Beiträge: 36 Registriert: 15.04.2016
|
erstellt am: 15. Jun. 2016 12:10 <-- editieren / zitieren --> Unities abgeben:
Hallo nochmal, es hat sich im Verlauf meiner Arbeit herausgestellt, dass ich doch auf die Vernetzung mit quaderförmigen Volumenelementen zurückgreifen muss. Das ist dadurch bedingt, dass diese in LS-DYNA gebraucht werden. Weiterhin ist auch eine aufsteigende Knotennummerierung sehr wichtig. Soweit scheint es zu funktionieren, doch ich habe noch das Problem, dass Ansys mir die Geometrien bei verschiedenen Winkeln herausgibt. Kann mir da noch jemand einen Tipp geben wie ich meinen Code anpasse, dass ich ein quaderförmiges Volumenelementnetz bei verschiedenen Winkeln meiner Kerbe erhalte? Vielen Dank & beste Grüße /BATCH /COM,ANSYS RELEASE 14.0 UP20111024 11:04:25 04/17/2013 /input,menust,tmp, ,,,,,,,,,,,,,,,,1 /GRA,POWER /GST,ON /PLO,INFO,3 /GRO,CURL,ON *SET,_REMOTE_VIS_ID,'21' /CPLANE,1 /REPLOT,RESIZE WPSTYLE,,,,,,,,0 finish /clear /COM,ANSYS RELEASE 14.0 UP20111024 11:04:58 04/17/2013 /input,start140,ans,/mnt/appl/x86_64/ansys_inc/v140/ansys/apdl/ , fini /title,Quader ! Eingabegrößen L=200 ! Probenlänge [mm] B=25 ! Probenbreite [mm] D2=4 ! Probendicke [mm] ! Normbedingt konstant AK=12.5 ! Abstand der Kerben [mm] LK=2 ! Länge der Kerben [mm] KF=0.45 ! Faktor der Kerbe [-] ! zwischen 0.45 ... 0.48 pi=4*atan(1) ! Pi ! Resultierende Größen D1=KF*D2 ! Probendicke - Kerbtiefe [mm] LangX=0.5*L+0.5*AK ! Lange Seite X [mm] KurzX=0.5*L-0.5*AK-LK ! Kurze Seite X [mm] ! Kerbwinkel Alpha=45 ! Kerbwinkel [°] BM=((Alpha*pi)/(180)) ! Winkel im Bogenmaß
! Elementgrößen elementgroesse1=2 !Elemtgröße Grob -> Außen elementgroesse2=2 !Elemtgröße Fein -> Kerbe elementgroesse3=2 !Elemtgröße Geometrie 3 -> Kohäsivelemente eg1=2 eg2=2 eg3=2 ! Elementanzahl elementanzahl1=25 ! Elementanzahl 1 -> Grob - Außen elementanzahl2=50 ! Elementanzahl 2 -> Fein - Kerbe elementanzahl3=20 ! Elementanzahl 3 -> Kohäsivelemente /prep7 ! Starten Preprocess
!----------------------------------------- C*** Zuweisungen für Konvertierung !-----------------------------------------
! Für Deckschichten - Real Constant = 50 TYPE,1 ! Auswahl des 1. Elementtyps SECNUM,1 ! Auswahl der 1. Section MAT,1 ! Auswahl des 1. Materials ! Für Kohäsivelemente - Real Constant = 100 TYPE,2 ! Auswahl des 2. Elementtyps SECNUM,2 ! Auswahl der 2. Section MAT,2 ! Auswahl des 2. Materials ! Elementzuweisung
et,1,solid186 ! Elementtypauswahl !----------------------------------------- C*** Geometrie 1 !-----------------------------------------
k,1,0,0,0 ! Keypoints 1-8 definieren Volumen 1 k,2,KurzX,0,0 k,3,KurzX,B,0 k,4,0,B,0 k,5,0,0,D2 k,6,KurzX,0,D2 k,7,KurzX,B,D2 k,8,0,B,D2 v,1,2,3,4,5,6,7,8 k,9,0,0,0 ! Keypoints 9-16 definieren Volumen 2 k,10,LangX-D1*tan(BM),0,0 k,11,LangX-D1*tan(BM),B,0 k,12,0,B,0 k,13,0,0,D1 k,14,LangX,0,D1 k,15,LangX,B,D1 k,16,0,B,D1 v,9,10,11,12,13,14,15,16 !----------------------------------------- C*** Geometrie 2 (Kleber) !-----------------------------------------
k,101,KurzX+LK,0,0 k,102,LangX,0,0 k,103,LangX,B,0 k,104,KurzX+LK,B,0 k,105,KurzX+LK,0,D2 k,106,LangX,0,D2 k,107,LangX,B,D2 k,108,KurzX+LK,B,D2 v,101,102,103,104,105,106,107,108 !----------------------------------------- C*** Geometrie 3 !-----------------------------------------
k,201,LangX+LK,0,0 k,202,L,0,0 k,203,L,B,0 k,204,LangX+LK,B,0 k,205,LangX+LK,0,D2 k,206,L,0,D2 k,207,L,B,D2 k,208,LangX+LK,B,D2 v,201,202,203,204,205,206,207,208 k,209,KurzX+LK,0,D2-D1 k,210,L,0,D2-D1 k,211,L,B,D2-D1 k,212,KurzX+LK,B,D2-D1 k,213,KurzX+LK+D1*tan(BM),0,D2 k,214,L,0,D2 k,215,L,B,D2 k,216,KurzX+LK+D1*tan(BM),B,D2
v,209,210,211,212,213,214,215,216 !----------------------------------------- C*** Überlappung !----------------------------------------- vovlap,all numcmp,volu !----------------------------------------- C*** Vernetzung !-----------------------------------------
!-------------------------------------------- C*** DECKSCHICHT 1 !-------------------------------------------- Pr1=30000 ! Beginn der Objektnummerierung numstr,node,Pr1+1 ! Einschalten der Knotennummerierung vsel,s,volu,,1,3 ! Deckschicht 1 Zuweisung vsel,a,volu,,7 vatt,1,50 allsel esize,eg3 ! Geometrie 3 Vernetzung type,1 mat,1 vmesh,3 esize,eg3 ! Geometrie 1 und 2 Vernetzung type,1 mat,1 vmesh,1,2
esize,eg3 !Geometrie 7 Vernetzung vmesh,7
allsel nummrg,all !----------------------------------------- C*** KOHÄSIVSCHICHT !----------------------------------------- Pr0=60000 numstr,node,Pr0+1 ! Einschalten der Knotennummerierung vsel,s,volu,,8 !Kohäsivschicht Zuweisung vatt,3,100 allsel esize,eg3 !Kohäsivschicht (Geometrie 8) Vernetzung type,2 mat,2 vmesh,8 !-------------------------------------------- C*** DECKSCHICHT 2 !--------------------------------------------
Pr2=100000 ! Beginn der Objektnummerierung der Papierprobe numstr,node,Pr2+1 ! Einschalten der Knotennummerierung vsel,s,volu,,4,6 ! Deckschicht 2 Zuweisung vsel,a,volu,,9 vatt,2,60 allsel esize,eg3 ! Geometrie 4 Vernetzung type,1 mat,1 vmesh,4 esize,eg3 ! Geometrie 9 Vernetzung vmesh,9
esize,eg3 ! Geometrie 5 und 6 Vernetzung type,1 mat,1 vmesh,5,6 allsel nummrg,all /eof
[Diese Nachricht wurde von Deffensor am 15. Jun. 2016 editiert.] Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP |
Rainer Schulze Ehrenmitglied V.I.P. h.c. Dipl.-Ing. im Ruhestand
Beiträge: 4419 Registriert: 24.09.2012
|
erstellt am: 15. Jun. 2016 12:32 <-- editieren / zitieren --> Unities abgeben: Nur für Deffensor
>>...die Vernetzung mit quaderförmigen Volumenelementen ... bei verschiedenen Winkeln meiner Kerbe Wie passt das zu den Regeln der Geometrie? Quader haben immer Winkel von 90°. Oder brauchst Du eigentlich nur Oktaeder mit überwiegend parallelen Flächen (verzerrte Quader)? ------------------ Rainer Schulze Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP |
Deffensor Mitglied Student
Beiträge: 36 Registriert: 15.04.2016
|
erstellt am: 15. Jun. 2016 12:36 <-- editieren / zitieren --> Unities abgeben:
Soweit wie möglich möchte ich Quader (90°) anwenden. Natürlich werde ich an den Kerben entsprechend verzerrte Quader-Elemente verwenden müssen. Das wichtigste hierbei ist für mich das die Knoten der Deckschichten mit den Knoten der Klebschicht/Mittelschicht übereinstimmen. [Diese Nachricht wurde von Deffensor am 15. Jun. 2016 editiert.] Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP |
wosch Ehrenmitglied V.I.P. h.c. Elektrotechniker im Ruhestand
Beiträge: 2607 Registriert: 16.12.2004 Rechne zuerst ein Problem nach, für das Du eine analytische Lösung kennst.
|
erstellt am: 15. Jun. 2016 16:35 <-- editieren / zitieren --> Unities abgeben: Nur für Deffensor
Wenn man von quaderförmigen Elementen spricht, meint man damit Hexaeder. Für die Nutzung mit ANSYS ist keine Rechtwinkligkeit erforderlich. (Eigentlich versteht sich das von selbst.) ------------------ Viel Erfolg wünscht Wolfgang Schätzing Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP |
wosch Ehrenmitglied V.I.P. h.c. Elektrotechniker im Ruhestand
Beiträge: 2607 Registriert: 16.12.2004 Rechne zuerst ein Problem nach, für das Du eine analytische Lösung kennst.
|
erstellt am: 15. Jun. 2016 18:41 <-- editieren / zitieren --> Unities abgeben: Nur für Deffensor
Nach dem konstruktiven Beitrag von Rainer Schulze nun auch noch ein nicht lauffähiger Quelltext !!! (Wir sind von unserem Moderator angehalten, längere Quelltexte in den Anhang zu schieben.) Zum Problem: Wenn eine sinnvole Reihenfolge bei der Vernetzung der Volumen eingehalten wird funktioniert's. ------------------ Viel Erfolg wünscht Wolfgang Schätzing Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP |
Deffensor Mitglied Student
Beiträge: 36 Registriert: 15.04.2016
|
erstellt am: 15. Jun. 2016 18:55 <-- editieren / zitieren --> Unities abgeben:
Hallo, vielen Dank für den Hinweis. Heißt es, dass die Geometrie richtig erstellt wird, ich jedoch einfach die Reihenfolge austauschen muss? Ich dachte jedoch, dass der Fehler bei der Geometrierzeugung liegt, da ich beim Setzen des /eof (end-of-file)- Befehles nach vovlap bei Winkeln ungleich 0 Grad eine falsche Geometrie erhalte. Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP |
Deffensor Mitglied Student
Beiträge: 36 Registriert: 15.04.2016
|
erstellt am: 15. Jun. 2016 21:37 <-- editieren / zitieren --> Unities abgeben:
Eigentlich müsste der Quelltext jedoch voll funktionsfähig sein. Also für den Winkel von 0 Grad. Bei anderen Winkeln habe ich ja gerade das Problem, dass das Modell falsche Formen annimmt. Anders gesagt liegt der Fehler bei vovlap, dass eine falsche Kontur erzeugt wird. Siehe Screenshot. [Diese Nachricht wurde von Deffensor am 16. Jun. 2016 editiert.] Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP |
Deffensor Mitglied Student
Beiträge: 36 Registriert: 15.04.2016
|
erstellt am: 16. Jun. 2016 11:07 <-- editieren / zitieren --> Unities abgeben:
Ich kann natürlich anders vorgehen und versuchen jedes Teilvolumen selbst zu vernetzen dann passt die Geometrie, doch das Netz ist dann nicht "stetig". Damit meine ich, dass die Knoten der Teilvolumen nicht miteinander verschmelzt sind. Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP |
Deffensor Mitglied Student
Beiträge: 36 Registriert: 15.04.2016
|
erstellt am: 16. Jun. 2016 16:12 <-- editieren / zitieren --> Unities abgeben:
|
| Anzeige.:
Anzeige: (Infos zum Werbeplatz >>)
|