| | | Gut zu wissen: Hilfreiche Tipps und Tricks aus der Praxis prägnant, und auf den Punkt gebracht für Ansys |
Autor
|
Thema: Berechnung einfache Membran (5500 mal gelesen)
|
heunigreenfreak Mitglied wiss. Mitarbeiter
Beiträge: 27 Registriert: 22.12.2009 Win XP, Inventor Professionel 2010 und 2012, ANSYS v15, LS-DYNA
|
erstellt am: 22. Dez. 2009 11:06 <-- editieren / zitieren --> Unities abgeben:
Hallo Zusammen, eine Frage, kennt sich jemand mit der Nutzung von Membranelementen shell41 aus? Ein ganz einfaches Beispiel funktioniert schon nicht, aus unerklärlichen Gründen??? Rechteckmembran, an vier Ecken eingespannt, mit Druck belastet. Es soll einfach nur die Deformation erhalten werden. Habe auch schon mit Vorspannung (Verschiebungsvorgabe an Eckpunkt ->siehe Beispiel unten) im ersten Lösungsschritt gearbeitet. Weiss jemand an was es liegt, was bitte mache ich falsch? Vielen Dank!! ps.: Es gibt dazu auch ein Beispiel auf CAD-FEM, das im Prinzip nichts anders beeinhaltet (Verschiebungsvorgabe für Vorspannung+ anschließend Druck-RB)
Code:
fini /clear /PREP7 et,1,shell41, r,1,1.5 mp,ex,1,2.1e5 mp,prxy,1,.3 rect,0,100,0,100 esize,4 MOPT,AMESH,MAIN MSHAPE,1,2D amesh,all save /solu nlgeom,on ! 1.Berechung d,1,all d,2,all d,52,all d,27,ux,.05 d,27,uy,.05 d,27,uz, alls solve ! 2.Berechung sf,all,pres,0.5 alls solve fini /post1 PLDISP,2 PRNSOL,U,COMP
Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP |
simulmi Mitglied Berechnungsingenieur
Beiträge: 295 Registriert: 27.08.2008 ANSYS Workbench 13 ANSYS Classic 13 MAXWELL 2D/3D V14 Pro/E Wildfire<P>Windows XP 64 bit
|
erstellt am: 22. Dez. 2009 13:30 <-- editieren / zitieren --> Unities abgeben: Nur für heunigreenfreak
Hallo, dazu bitte im Solution-File nachsehen mit welcher Meldung der Solver abbricht und ggf. hier posten. Gruß simulmi ------------------ Matthias Ulmer Institut für Feinwerktechnik Universität Stuttgart Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP |
heunigreenfreak Mitglied wiss. Mitarbeiter
Beiträge: 27 Registriert: 22.12.2009 Win XP, Inventor Professionel 2010 und 2012, ANSYS v15, LS-DYNA
|
erstellt am: 22. Dez. 2009 14:31 <-- editieren / zitieren --> Unities abgeben:
Die Lösung konvergiert nicht (Kraftresiduum), Verschiebung ja, aber nur extrem langsam. Bei diesem einfachen Beispiel sind eigentlich keine derart schwierigen Nichtlinearitäten vorhanden. Für einen Kreis bekommt man ja sogar die Lösung noch relativ einfach per Hand(Poisson-PDGL). Es muss also an etwas anderem liegen??? Error-File-Ausgabe:
Code: *** WARNING *** SUPPRESSED MESSAGE CP = 188.172 TIME= 14:24:14 There are 3 small equation solver pivot terms. *** ERROR *** SUPPRESSED MESSAGE CP = 188.172 TIME= 14:24:14 Solution not converged at time 1.01 (load step 2 substep 1). Run terminated. *** WARNING *** SUPPRESSED MESSAGE CP = 188.172 TIME= 14:24:14 The unconverged solution (identified as time 2 substep 999999) is output for analysis debug purposes. Results should not be used for any other purpose.
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: 31. Dez. 2009 18:15 <-- editieren / zitieren --> Unities abgeben: Nur für heunigreenfreak
|
heunigreenfreak Mitglied wiss. Mitarbeiter
Beiträge: 27 Registriert: 22.12.2009 Win XP, Inventor Professionel 2010 und 2012, ANSYS v15, LS-DYNA
|
erstellt am: 04. Jan. 2010 09:18 <-- editieren / zitieren --> Unities abgeben:
Gesundes Neues und Dank für die bisherigen Vorschläge, auch wenn noch nicht zielführend... Aber nein, es geht so nicht! Denn: 1.) Ich meinte nur den Membranspannungszustand. Der Biegeanteil spielt nur bei einem Plattenelement eine Rolle, bei einem reinen Membranelement nicht (verhält sich quasi wie ein 2-dimensionales Seil), die Lösung im "wosch-Beispiel" basiert aber nur darauf. Folglich ist beim shell63 keyopt,1,1,1 zu setzen oder gleich shell41 zu nehmen. Und dann kommt eben keine Verschiebung raus bei kleinen Def. oder konvergiert nicht bei großen Def. 2.) Große Deformationen sind zu betrachten! Hat denn niemand eine Idee an was es liegen könnte? Danke schonmal für die Hilfe. Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP |
ChristophN Mitglied Berechnungsingenieur
Beiträge: 773 Registriert: 13.12.2008
|
erstellt am: 04. Jan. 2010 09:47 <-- editieren / zitieren --> Unities abgeben: Nur für heunigreenfreak
|
heunigreenfreak Mitglied wiss. Mitarbeiter
Beiträge: 27 Registriert: 22.12.2009 Win XP, Inventor Professionel 2010 und 2012, ANSYS v15, LS-DYNA
|
erstellt am: 04. Jan. 2010 13:29 <-- editieren / zitieren --> Unities abgeben:
Hi Christoph...Gruss an die TUC Mechaniker nsubst hilft auch nicht weiter. Wenn ich nachstehendes ganz einfaches Beispiel verwende, treten zwar die Konvergenzprobleme nicht auf (wie bei der Vorgehensweise Vorspannung+Druckbelastung), aber die Verschiebungen aller Knoten sind im Ergebnis = null. Wird die Belastung irgendwie nicht korrekt aufgebracht? Der Fehler mit den pivot-Elementen tritt bei diesem Bsp auch nicht auf. Code: fini /clear /PREP7 et,1,shell41, r,1,1.5 mp,ex,1,2.1e5 mp,prxy,1,.3 rect,0,100,0,100 esize,4 MOPT,AMESH,MAIN MSHAPE,1,2D amesh,all save /solu nlgeom,on sf,all,pres,0.5 !Druck sftran d,1,all ! alle Ecken festhalten d,2,all d,52,all d,27,all alls solcontrol,off nsubst,15 solve fini /post1 PLDISP,2 PRNSOL,U,COMP
Error-File-Ausgabe (wiederholende Fehler dieser Art):
Code: *** WARNING *** CP = 6.984 TIME= 13:15:38 Structural elements without mid nodes usually produce much more accurate results in quad or brick shape. *** WARNING *** CP = 7.797 TIME= 13:15:41 The calculated reference FORCE CONVERGENCE VALUE = 0 is less than the default minimum of 1 (MINREF field of CNVTOL command). Check results carefully.
Gruss Christian ps.:Eigentlich ist das Bsp. total simpel, irgendwo ist da noch ein kleiner Fehler. (hab aber auch nicht Zeit mir h damit zu vertun) ------------------ "Die Welt ist nichtlinear" Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP |
heunigreenfreak Mitglied wiss. Mitarbeiter
Beiträge: 27 Registriert: 22.12.2009 Win XP, Inventor Professionel 2010 und 2012, ANSYS v15, LS-DYNA
|
erstellt am: 07. Jan. 2010 13:13 <-- editieren / zitieren --> Unities abgeben:
Niemand eine Idee? Es muss doch in ANSYS für eine einfache druckbelastete Membran nicht vorgespannt/vorgespannt eine Lösung geben! Oder, was ist anders zum Beispiel von CAD-FEM (siehe Erster Beitrag). Danke schonmal! Grüße ------------------ "Die Welt ist nichtlinear" Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP |
Vic09 Mitglied Student
Beiträge: 5 Registriert: 11.01.2010
|
erstellt am: 11. Jan. 2010 09:59 <-- editieren / zitieren --> Unities abgeben: Nur für heunigreenfreak
Hallo, ich hatte genau das gleiche Probelm wie du, habs jetzt aber (zumindest halbwegs) hingekriegt. Im Grunde hab ich es genau wie du gemacht, nur mit anderen BCs. Ich hab nämlich ein Line Pressure aufgebracht zur Vorspannung und desweiteren noch die Bewegung in Normalenrichtung an den Rändern zu Null gesetzt. Dann hat es bei mir funktioniert, jedoch nur bei Element Shell41. Bei Shell63 mit keyopt,1,1,1 klappt es aus irgendwelchen Gründen nicht!! Ich hoffe es hilft dir weiter... Viele Grüße, Victor Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP |
heunigreenfreak Mitglied wiss. Mitarbeiter
Beiträge: 27 Registriert: 22.12.2009 Win XP, Inventor Professionel 2010 und 2012, ANSYS v15, LS-DYNA
|
erstellt am: 11. Jan. 2010 13:05 <-- editieren / zitieren --> Unities abgeben:
Danke für den Hinweis. So richtig verständlich ist mir das aber nicht. Macht denn ein Line Pressure auf festgehaltene Ränder Sinn? Erwartungsgemäß kommt da doch für alle Knoten null heraus. Oder wie meintest du das? Dein skript? ------------------ "Die Welt ist nichtlinear" Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP |
Vic09 Mitglied Student
Beiträge: 5 Registriert: 11.01.2010
|
erstellt am: 11. Jan. 2010 13:37 <-- editieren / zitieren --> Unities abgeben: Nur für heunigreenfreak
Line Pressure (oder alternativ Kräfte an den Knoten) macht dann Sinn, wenn ich in x- und z- richtung vorspanne und die Verschiebung in y-richtung verhindere. Blöde Frage: Wie erstelle ich ein Script in Ansys? Viele Grüße, Victor Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP |
heunigreenfreak Mitglied wiss. Mitarbeiter
Beiträge: 27 Registriert: 22.12.2009 Win XP, Inventor Professionel 2010 und 2012, ANSYS v15, LS-DYNA
|
erstellt am: 12. Jan. 2010 08:29 <-- editieren / zitieren --> Unities abgeben:
Ok, nur mit Line Pressure muss aber zwingend mindestens ein weiterer Knoten in x und z festgehalten werden, das Kräftegleichgewicht ist in der FEM ja nie genau erfüllt, die Membran würde sich ins unendliche verschieben. So richtig verständlich ist das nicht. Hättest du das Eingabe-Skript? In was arbeitest du denn, Ansys Classic oder Workbench? Hast du dazu kein Skript per Handeingabe? (erstellen geht in Classic über List->File->Logfile oder in Workbench über Extras->Schreibe Ansys-Input-File (oder so ähnlich), dazu erst auf die Lösung klicken) Grüße ------------------ "Die Welt ist nichtlinear" Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP |
Vic09 Mitglied Student
Beiträge: 5 Registriert: 11.01.2010
|
erstellt am: 12. Jan. 2010 09:36 <-- editieren / zitieren --> Unities abgeben: Nur für heunigreenfreak
Hallo, ich hab jetzt mal ein log file erstellt. Ich hoffe es hilft weiter. Grüße, Victor /BATCH ! /COM,ANSYS RELEASE 12.0.1 UP20090415 09:30:43 01/12/2010 /input,menust,tmp,'',,,,,,,,,,,,,,,,1 ! /GRA,POWER ! /GST,ON ! /PLO,INFO,3 ! /GRO,CURL,ON ! /CPLANE,1 ! /REPLOT,RESIZE WPSTYLE,,,,,,,,0 !* /NOPR /PMETH,OFF,0 KEYW,PR_SET,1 KEYW,PR_STRUC,1 KEYW,PR_THERM,0 KEYW,PR_FLUID,0 KEYW,PR_ELMAG,0 KEYW,MAGNOD,0 KEYW,MAGEDG,0 KEYW,MAGHFE,0 KEYW,MAGELC,0 KEYW,PR_MULTI,0 KEYW,PR_CFD,0 /GO !* ! /COM, ! /COM,Preferences for GUI filtering have been set to display: ! /COM, Structural !* /prep7 et,1,43 r,1,0.001 mp,ex,1,2.1e11 mp,prxy,1,.3 mp,dens,1,7850 rectng,0,0.1,0,0.1 esize,0.01 MOPT,AMESH,MAIN MSHAPE,1,2D amesh,all FLST,2,4,4,ORDE,2 FITEM,2,1 FITEM,2,-4 !* /GO DL,P51X, ,UZ,0 FLST,2,11,1,ORDE,4 FITEM,2,1 FITEM,2,22 FITEM,2,32 FITEM,2,-40 !* /GO F,P51X,FX,-1000 FLST,2,11,1,ORDE,3 FITEM,2,2 FITEM,2,12 FITEM,2,-21 !* /GO F,P51X,FX,1000 FLST,2,11,1,ORDE,3 FITEM,2,12 FITEM,2,22 FITEM,2,-31 !* /GO F,P51X,FY,1000 FLST,2,11,1,ORDE,2 FITEM,2,1 FITEM,2,-11 !* /GO F,P51X,FY,-1000 sf,all,pres,0.5 FINISH /SOL nlgeom,on solve FINISH /POST1 ! /SHOW,WIN32C SET,FIRST ! /PLOPTS,INFO,3 ! /CONTOUR,ALL,18 ! /PNUM,MAT,1 ! /NUMBER,1 ! /REPLOT,RESIZE ! PLNSOL,U,Z ! /UIS,ABORT,1 ! /SHOW,WIN32 ! /REPLOT,RESIZE ! LGWRITE,'membran','lgw','D:\Victor\ANSYSW~1\',COMMENT
Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP |
heunigreenfreak Mitglied wiss. Mitarbeiter
Beiträge: 27 Registriert: 22.12.2009 Win XP, Inventor Professionel 2010 und 2012, ANSYS v15, LS-DYNA
|
erstellt am: 13. Jan. 2010 23:50 <-- editieren / zitieren --> Unities abgeben:
Also habe es versucht, bei mir kommt da keine Lösung heraus. Bitte auch beachten! Dein Element ist das shell43! Das ist NICHT das Membran-Element. Grüße [Diese Nachricht wurde von heunigreenfreak am 14. Jan. 2010 editiert.] Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP |
Vic09 Mitglied Student
Beiträge: 5 Registriert: 11.01.2010
|
erstellt am: 14. Jan. 2010 16:40 <-- editieren / zitieren --> Unities abgeben: Nur für heunigreenfreak
Hallo, ja da hast du recht! Ich hab mich wohl einfach vertippt. Nichts desto trotz hab ich die Membran mit shell63 keyopt 1,1,1 erstellt. Ich werd in den nächsten Tagen nochmal ein vernünftiges Skript schreiben... Viele Grüße, Victor Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP |
| Anzeige.:
Anzeige: (Infos zum Werbeplatz >>)
|