---

Hallo Leute,

ich möchte gerne eine nichtlineare Festigkeitsberechnung eines Hakens durchführen. Es sollen dabei kleine Dehnungen angenommen werden, die Nichtlinearität des Materials werden einer Datenbank entnommen. Weiterhin will ich eine Abfolge von 4 Lastschritten simulieren, die Kraftangriffsfläche befindet sich an der Innenseite.
Soweit habe ich eigentlich kein Problem, jedoch stürzt mein Ansys ab, sobald ich im /solu solve eingebe.

Mein Makro möchte ich auch gerne mal anfügen:

Code:

---

fini
/cle
l=40
l1=l/2
b1=26
b2=24
r=50
rf=8
h1=15
h2=10
h3=20
a=60
b=80

F1=10000
F2=100000
F=0

/prep7

et,1,solid185
et,10,mesh200,6

!MP, ex, 1, 210000
!MP, nuxy, 1, 0.3

tb,mkin,1
tbtemp,,strain
tbdata,1,0.002,0.006,0.01,0.02,0.05
tbtemp,
tbdata,1,420,540,600,660,710
!tbplot,mkin,1

/pnum,kp,1
/pnum,line,1
/pnum,area,1

csys,1
k,1,r,180
k,2,r,-90
k,3,r,0
k,4,r,a
csys,0
k,5,0,2*r,0
k,6,0,3*r,0
k,7,

l,6,5
larc,1,2,7,r
larc,2,3,7,r
larc,3,4,7,r
l2tan,4,-1
lsel,s,line,,1
ldele,all,,,1
ksel,s,kp,,7
kdele,all
alls
lcomb,all
numcmp,kp

k,3,-r+l1,0,b1/2
k,4,-r+l1,0,-b1/2
k,5,-r-l1,0,-b1/2
k,6,-r-l1,0,b1/2

l,3,4
l,4,5
l,5,6
l,3,6

lfillt,1,3,rf
lfillt,4,3,rf
lfillt,4,5,rf
lfillt,1,5,rf
lsel,s,line,,4
lsel,a,line,,7
lsel,a,line,,8
lcomb,all
lsel,s,line,,1
lsel,a,line,,6
lsel,a,line,,9
lcomb,all

alls
lesize,1,,,5
lesize,3,,,16

type, 10
mshkey,2
al,1,3,4,5
amesh,all
nummrg,all

type,1
mat,1
esize,,40
vdrag,1,,,,,,2,

esize,,5
vext,1,,,,h3,,0.8,,0.8
esize,,5
vext,6,,,,h1,,0.6,,0.6
esize,,2
vext,12,,,,h2

fini

/solu

nsel,s,loc,y,2*r+h1+h2
d,all,all,0

nlgeom,1

!Kraft 1
time,0
csys,1
nsel,s,loc,x,r
nsel,r,loc,y,-90+b/2,-90-b/2
f,all,fy,-f1
solve

!Kraft 2
time,2
fdele
solve

!Kraft 3
time,4
nsel,s,loc,x,r
nsel,r,loc,y,-90+b/2,-90-b/2
f,all,fy,-f2
solve

!Kraft 4
time,6
fdele
solve

fini

/post1

plnsol,s,eqv
plnsol,eppl,eqv

---

Funktioniert das bei euch? Habe ich gravierende Fehler im Code?

Ich hoffe, ihr könnt mir weiterhelfen.

Gruß
Wildy

Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP

---

Auf den ersten Blick fällt das fehlende allsel vor dem solve auf. Glaube allerdings nicht, dass das die Hauptursache für den Absturz ist.
Was sagt denn das error-file?

Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP

---

Bei mir verabschiedet sich Ansys auch ohne Kommentar. Es gibt also kein errorlog.
Durch das Einfügen von alls funktioniert es. und

!MP, ex, 1, 210000
!MP, nuxy, 1, 0.3

sollte ohne Kommentar da stehen.

------------------
Martin Feierabend
Fachgebiet Getriebetechnik
Technische Universität Ilmenau

[Diese Nachricht wurde von medin am 08. Jan. 2010 editiert.]

Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP

---

Das error-file ist im Arbeitsverzeichnis als 'filnemame'.err abgespeichert.
Gruß, A.

Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP

---

Das ist mir schon klar, wenn sich Ansys aber einfach selbst beendet und das ohne Fehlermeldung, dann steht da nichts drin.

------------------
Martin Feierabend
Fachgebiet Getriebetechnik
Technische Universität Ilmenau

Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP

---

Erst einmal danke euch beiden.

Das Einfügen von alls hat in der Tat geholfen. Ansys stürzt tatsächlich nicht mehr ab. Dass ich die mp Befehle auskommentiert habe, ist mir schleierhaft. Dass ich das die ganze Zeit nicht bemerkt habe, ebenso!

Wenn ich bis zum 1. solve Befehl bearbeiten lasse, dann kommt ein Popup mit "A check of your load data produced 1 warnings. SHOULD THE SOLV COMMAND BE EXECUTED?"

Weiterhin kann ich bei der Ergebnisauswertung weder plnsol,s,eqv noch plnsol,eppl,eqv ausführen, weil der Befehl ignoriert wird.

Bin da irgendwie ratolos.

Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP

---

deine Lösung konvergiert schon im ersten lastschritt nicht. Eventuell schaffen da die befehle deltime und nequit abhilfe. falls nicht, ist die aufgebrachte kraft einfach zu groß.

------------------
Martin Feierabend
Fachgebiet Getriebetechnik
Technische Universität Ilmenau

Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP

---

Hallo medin,

leider konvergiert kein einziger Lastfall, der Grund dafür ist "DOF limit exceeded". Ist vielleicht die Auswahl des solid185 nicht geeignet? Vielleicht deshalb die Fehlermeldung mit dem Freiheitsgrad DOF.
Auch die Berechnung dauert etwa 30 Minuten, das kann so nicht sein. Die verwendeten Kräfte sind aus einer FEM Aufgabenstellung entnommen und können nicht Grund dafür sein, dass die Simulation fehlschlägt.
Einen kleinen Fehler habe ich noch gefunden: Die Knotenauswahl für die Kraftangriffspunkte muss jeweils

Code:

---

nsel,s,loc,x,r-l1

---

heißen, jedoch hat dies keinen Einfluss auf die Lösung.

Mit den beiden Befehlen nequit und deltime kann ich leider nichts anfangen, weil die in der Hilfe nicht vorhanden sind.

Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP

---

"nequit" ist unter 'NEQIT' in der Hilfe, DELTIM so wie geschrieben. Sieh mal in der Hilfe von ANSYS 11 nach, da ist manches (meiner Meinung nach) leichter zu finden.
Versuche die Lastschritte zu verkleinern (mit neqit, deltim), möglicherweise konvergieren die Lastschritte dann besser.

Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP

---

Hallo Wildy!

unter der gegebenen Last fließt der Haken einfach weg. mit dem effekt das deine Elemente extrem verzerrt werden. Deshalb wird die Rechnung abgebrochen und der Lastschritt konvergiert nicht. Wenn das Material zu fließen beginnt muss die Last in kleinen Schritten aufgebracht werden.
Du solltst die Last per vorgabe nsubst,10 (oder besser noch mehr) langsam aufbringen. Bei mir läuft die rechnung noch und ist akutell bei einer Last von 0.33*Endlast. Dabei fließts anscheinend schon ordentlich.

Gruß
Jens

------------------
Jens Friedrich
TU-Dresden

Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP

---

Haken.png

ERROR IN ELEMENT FORMULATION!!!
Ich hab gerade gesehen dass die Last nicht verteilt wurde. Auf jeden Knoten hast du die gesamte Last aufgebracht. damit ergibt sich logischerweise eine riesige Last welche unrealistische Ergebnisse erzeugt.

so läufts
:
fini
/cle
l=40
l1=l/2
b1=26
b2=24
r=50
rf=8
h1=15
h2=10
h3=20
a=60
b=80
F1=10000
F2=100000
F=0

/prep7

et,1,solid185
et,10,mesh200,6

MP, ex, 1, 210000
MP, nuxy, 1, 0.3

tb,mkin,1
tbtemp,,strain
tbdata,1,0.002,0.006,0.01,0.02,0.05
tbtemp,
tbdata,1,420,540,600,660,710
!tbplot,mkin,1

/pnum,kp,1
/pnum,line,1
/pnum,area,1

csys,1
k,1,r,180
k,2,r,-90
k,3,r,0
k,4,r,a
csys,0
k,5,0,2*r,0
k,6,0,3*r,0
k,7,

l,6,5
larc,1,2,7,r
larc,2,3,7,r
larc,3,4,7,r
l2tan,4,-1
lsel,s,line,,1
ldele,all,,,1
ksel,s,kp,,7
kdele,all
alls
lcomb,all
numcmp,kp

k,3,-r+l1,0,b1/2
k,4,-r+l1,0,-b1/2
k,5,-r-l1,0,-b1/2
k,6,-r-l1,0,b1/2

l,3,4
l,4,5
l,5,6
l,3,6

lfillt,1,3,rf
lfillt,4,3,rf
lfillt,4,5,rf
lfillt,1,5,rf
lsel,s,line,,4
lsel,a,line,,7
lsel,a,line,,8
lcomb,all
lsel,s,line,,1
lsel,a,line,,6
lsel,a,line,,9
lcomb,all

alls
lesize,1,,,5
lesize,3,,,16

type, 10
mshkey,2
al,1,3,4,5
amesh,all
nummrg,all

type,1
mat,1
esize,,40
vdrag,1,,,,,,2,

esize,,5
vext,1,,,,h3,,0.8,,0.8
esize,,5
vext,6,,,,h1,,0.6,,0.6
esize,,2
vext,12,,,,h2

fini

/solu

nsel,s,loc,y,2*r+h1+h2
d,all,all,0

nlgeom,1

!Kraft 1
time,0
csys,1
nsel,s,loc,x,r
nsel,r,loc,y,-90+b/2,-90-b/2
*get,NANZ,node,,count
f,all,fy,-f1/nanz
alls
nsubst,10
solve

!Kraft 2
time,2
fdele
solve

!Kraft 3
time,4
nsel,s,loc,x,r
nsel,r,loc,y,-90+b/2,-90-b/2
f,all,fy,-f2
solve

!Kraft 4
time,6
fdele
solve

fini

/post1

plnsol,s,eqv

------------------
Jens Friedrich
TU-Dresden

Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP

---

Hallo Jens,

ach herrje, wieder so ein dummer Fehler, den ich übersehen hab. Klar, die Last muss verteilt werden.
Ich werde ebenso mal deine Tips aus dem vorherigen Post berücksichtigen!
Mal schauen, ob ich es dann hinbekomme.
Großes Dankeschön schon mal.

Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP

---

Ist das eigentlich immer noch der Abschlussbeleg für die FEM-Praktika an der TU in Ilmenau? den hatte ich doch auch schon vor Jahren.

grüße Martin

------------------
Martin Feierabend
Fachgebiet Getriebetechnik
Technische Universität Ilmenau

Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP

---

Ja, ist immer noch der Fall.
Da ich aber den Beleg fertig machen wollte, bevor ich alle dafür notwendigen Befehle hatte, hat natürlich einiges nicht auf Anhieb geklappt. Aber mit eurer Hilfe gings dann natürlich.
Danke nochmal.

Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP