---

1.jpg

Hallo,
ich bin ein absoluter Ansysneulin und bin grad dabei mich in die materie hineinzuarbeiten. (Classic)
Folgendes Problem:
Ich möchte einen Kontakt zwischen einem Zylinder und einer Platte simulieren, habe das Modell aufgebaut(inkl. symmetrie) und es vernetzt (zugegeben ist das Netz noch nicht so doll und nicht das Entstadium).
Weiterhin ist der Kontakt mit dem kontaktmanager erstellt worden.

Definiere ich nun die Knoten der oberen waagerechte Fläche mit einer verschiebung in y-Richtung, so erhalte ich ein Ergebnis.
Quelltext:

/BATCH 
/COM,ANSYS RELEASE 13.0    UP20101012       08:16:25    05/04/2011
/input,menust,tmp,'',,,,,,,,,,,,,,,,1  
/GRA,POWER
/GST,ON
/PLO,INFO,3
/GRO,CURL,ON
/CPLANE,1  
/REPLOT,RESIZE 
WPSTYLE,,,,,,,,0
/REPLOT,RESIZE 

/prep7 

!Element definieren
!==================

ET,1,SOLID65
 
KEYOPT,1,1,0
KEYOPT,1,3,0
KEYOPT,1,5,2
KEYOPT,1,6,0
KEYOPT,1,7,0
KEYOPT,1,8,0

 
!Körper erstellen, Materialdaten
!===============================
 
mp,ex,1,200e9  
mp,prxy,1,0.3
 
cylind,,1,-0.5,0.5,,-90
block,0,2,-1,-2,-0.5,0.5

!Sicht ändern
!============
/VIEW,1,1,2,3  
/ANG,1 
/REP,FAST
 
!Vernetzen
!=========

lesi,16,,,10,0.2
lesi,11,,,10,0.2
lesi,20,,,5
lesi,8,,,5
lesi,7,,,5
lesi,9,,,5
lesi,3,,,10,5
lesi,4,,,10,5

MSHAPE,0,3D
MSHKEY,1
!* 
FLST,5,2,6,ORDE,2  
FITEM,5,1  
FITEM,5,-2 
CM,_Y,VOLU 
VSEL, , , ,P51X
CM,_Y1,VOLU
CHKMSH,'VOLU'  
CMSEL,S,_Y 
!* 
VMESH,_Y1  
!* 
CMDELE,_Y  
CMDELE,_Y1 
CMDELE,_Y2

!Kontaktbedingungen
!==================

CM,_TARGET,AREA
!* 
CM,_CONTACT,AREA
!* 
/COM, CONTACT PAIR CREATION - START
CM,_NODECM,NODE
CM,_ELEMCM,ELEM
CM,_KPCM,KP
CM,_LINECM,LINE
CM,_AREACM,AREA
CM,_VOLUCM,VOLU
/GSAV,cwz,gsav,,temp
MP,MU,1,1000
MAT,1  
R,3
REAL,3 
ET,2,170
ET,3,174
KEYOPT,3,9,0
KEYOPT,3,10,2  
R,3,
RMORE, 
RMORE,,0
RMORE,0
! Generate the target surface  
ASEL,S,,,6 
ASEL,A,,,7 
ASEL,A,,,8 
ASEL,A,,,9 
ASEL,A,,,10
ASEL,A,,,11
CM,_TARGET,AREA
TYPE,2 
NSLA,S,1
ESLN,S,0
ESLL,U 
ESEL,U,ENAME,,188,189  
NSLE,A,CT2 
ESURF  
CMSEL,S,_ELEMCM
! Generate the contact surface 
ASEL,S,,,1 
ASEL,A,,,2 
ASEL,A,,,3 
ASEL,A,,,4 
ASEL,A,,,5 
CM,_CONTACT,AREA
TYPE,3 
NSLA,S,1
ESLN,S,0
NSLE,A,CT2 ! CZMESH patch (fsk qt-40109 8/2008)
ESURF  
ALLSEL 
ESEL,ALL
ESEL,S,TYPE,,2 
ESEL,A,TYPE,,3 
ESEL,R,REAL,,3 
/PSYMB,ESYS,1  
/PNUM,TYPE,1
/NUM,1 
EPLOT  
ESEL,ALL
ESEL,S,TYPE,,2 
ESEL,A,TYPE,,3 
ESEL,R,REAL,,3 
CMSEL,A,_NODECM
CMDEL,_NODECM  
CMSEL,A,_ELEMCM
CMDEL,_ELEMCM  
CMSEL,S,_KPCM  
CMDEL,_KPCM
CMSEL,S,_LINECM
CMDEL,_LINECM  
CMSEL,S,_AREACM
CMDEL,_AREACM  
CMSEL,S,_VOLUCM
CMDEL,_VOLUCM  
/GRES,cwz,gsav 
CMDEL,_TARGET  
CMDEL,_CONTACT 
/COM, CONTACT PAIR CREATION - END  
/MREP,EPLOT 

!Aufbringen der Randbedingungen
!==============================

FINISH 
/SOL

asel,s,area,,11
nsla,s,1
d,all,ux

asel,s,area,,8
nsla,s,1
d,all,uy
d,all,uz

!
asel,s,area,,5
nsla,s,1
d,all,uy,-0.2
!

asel,s,area,,10
nsla,s,1
d,all,ux

asel,s,area,,4
nsla,s,1
d,all,ux
alls

nlgeom,on
auto,on
nsub,20,100
pred,on
solve

Wenn ich nun jedoch auf die linke obere Kante eine Last aufbringe, so:
-Konvergiert ANSYS nicht wenn die Kraft klein ist
-Ist die maximale Verschiebung eines Knotens in y-Richtung erreicht, was doch bedeutet, das der Körper abhaut!?
Quelltext:

/BATCH 
/COM,ANSYS RELEASE 13.0    UP20101012       08:16:25    05/04/2011
/input,menust,tmp,'',,,,,,,,,,,,,,,,1  
/GRA,POWER
/GST,ON
/PLO,INFO,3
/GRO,CURL,ON
/CPLANE,1  
/REPLOT,RESIZE 
WPSTYLE,,,,,,,,0
/REPLOT,RESIZE 

/prep7 

!Element definieren
!==================

ET,1,SOLID65
 
KEYOPT,1,1,0
KEYOPT,1,3,0
KEYOPT,1,5,2
KEYOPT,1,6,0
KEYOPT,1,7,0
KEYOPT,1,8,0

 
!Körper erstellen, Materialdaten
!===============================

 
mp,ex,1,200e9  
mp,prxy,1,0.3
 
cylind,,1,-0.5,0.5,,-90
block,0,2,-1,-2,-0.5,0.5

!Sicht ändern
!============
/VIEW,1,1,2,3  
/ANG,1 
/REP,FAST
 
!Vernetzen
!=========

lesi,16,,,10,0.2
lesi,11,,,10,0.2
lesi,20,,,5
lesi,8,,,5
lesi,7,,,5
lesi,9,,,5
lesi,3,,,10,5
lesi,4,,,10,5

MSHAPE,0,3D
MSHKEY,1
!* 
FLST,5,2,6,ORDE,2  
FITEM,5,1  
FITEM,5,-2 
CM,_Y,VOLU 
VSEL, , , ,P51X
CM,_Y1,VOLU
CHKMSH,'VOLU'  
CMSEL,S,_Y 
!* 
VMESH,_Y1  
!* 
CMDELE,_Y  
CMDELE,_Y1 
CMDELE,_Y2

!Kontaktbedingungen
!==================

CM,_TARGET,AREA
!* 
CM,_CONTACT,AREA
!* 
/COM, CONTACT PAIR CREATION - START
CM,_NODECM,NODE
CM,_ELEMCM,ELEM
CM,_KPCM,KP
CM,_LINECM,LINE
CM,_AREACM,AREA
CM,_VOLUCM,VOLU
/GSAV,cwz,gsav,,temp
MP,MU,1,1000
MAT,1  
R,3
REAL,3 
ET,2,170
ET,3,174
KEYOPT,3,9,0
KEYOPT,3,10,2  
R,3,
RMORE, 
RMORE,,0
RMORE,0
! Generate the target surface  
ASEL,S,,,6 
ASEL,A,,,7 
ASEL,A,,,8 
ASEL,A,,,9 
ASEL,A,,,10
ASEL,A,,,11
CM,_TARGET,AREA
TYPE,2 
NSLA,S,1
ESLN,S,0
ESLL,U 
ESEL,U,ENAME,,188,189  
NSLE,A,CT2 
ESURF  
CMSEL,S,_ELEMCM
! Generate the contact surface 
ASEL,S,,,1 
ASEL,A,,,2 
ASEL,A,,,3 
ASEL,A,,,4 
ASEL,A,,,5 
CM,_CONTACT,AREA
TYPE,3 
NSLA,S,1
ESLN,S,0
NSLE,A,CT2 ! CZMESH patch (fsk qt-40109 8/2008)
ESURF  
ALLSEL 
ESEL,ALL
ESEL,S,TYPE,,2 
ESEL,A,TYPE,,3 
ESEL,R,REAL,,3 
/PSYMB,ESYS,1  
/PNUM,TYPE,1
/NUM,1 
EPLOT  
ESEL,ALL
ESEL,S,TYPE,,2 
ESEL,A,TYPE,,3 
ESEL,R,REAL,,3 
CMSEL,A,_NODECM
CMDEL,_NODECM  
CMSEL,A,_ELEMCM
CMDEL,_ELEMCM  
CMSEL,S,_KPCM  
CMDEL,_KPCM
CMSEL,S,_LINECM
CMDEL,_LINECM  
CMSEL,S,_AREACM
CMDEL,_AREACM  
CMSEL,S,_VOLUCM
CMDEL,_VOLUCM  
/GRES,cwz,gsav 
CMDEL,_TARGET  
CMDEL,_CONTACT 
/COM, CONTACT PAIR CREATION - END  
/MREP,EPLOT 

!Aufbringen der Randbedingungen
!==============================

FINISH 
/SOL

asel,s,area,,11
nsla,s,1
d,all,ux

asel,s,area,,8
nsla,s,1
d,all,uy

asel,s,area,,2
nsla,s,1
d,all,uz

!
asel,s,area,,7
nsla,s,1
d,all,uz
!

asel,s,area,,10
nsla,s,1
d,all,ux

asel,s,area,,4
nsla,s,1
d,all,ux
alls

lsel,s,line,,9
nsll,s,1
f,all,fy,-1000
alls

nlgeom,on
auto,on
nsub,20,100

pred,on

solve

Meine Frage: Was mache ich hier falsch?

Danke für Antworten!!

[Diese Nachricht wurde von Neupi am 04. Mai. 2011 editiert.]

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

---

Prinzipiel ist es wenig ratsam Modell die nur durch kontakt in eine richtung gehalten werden kraftgesteuert zu berechnen, da man erst den Kontakt etablieren muss bevor dieser Das Bauteil halten kann. Bei Kraftsteuerung und vorallem deiner relativ grßen Schrittweite fliegt dir das Bauteil direkt weg.

ich würde wie folgt vorgehen:

1. Loadstep 1: Bauteil weggesteuert auf den kontaktpartner zu bewegen wir der Kontakt definitiv gefunden ist.
2: Loadstep 2: Reaktionskrauft aus der wegsteuerung auslesen und direkt ( stepped nicht ramped) als Kraft aufbringen. Der  ANSYS Befehl dazu ist mir gerande entfallen, aber es gibt ihn
3. Loadstep 3: Kraft bis zum gewünschten wert erhöhen.

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

---

Ok danke, ich werde mir die Hinweise in den nächsten Tagen erarbeiten. Wäre denn eine kleinere Schrittweite möglicherweise auch zielführend?

mfg

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

---

Hallo,

ich hätte folgendes anzumerken:

- Würde auch eine Wegsteuerung empfehlen, die ist im Konvergenzverhalten deutlich robuster.
- Kleinere Schrittweite halte ich auch für empfehlenswert
- Das Netz ist, gelinde gesagt, nicht so der Hit. Gerade im Kontaktbereich muss das Netz deutlich feiner werden. Dann klappts auch mit der Kovergenz.
- Dann wäre da noch die nächste Stellschraube, die Kontaktsteifigkeit. Probiere mal einen Wert von 0,1. (Nach kovergenter Lösung ist die Durchdringung auszuwerten, diese sollte in "einem vernünftig kleinen Bereich" liegen.)
- Und stelle bitte die Aktualisierung der Kontaktsteifigkeit so ein, dass diese bei jedem Iterationsschritt aktualisiert wird.

Und, klappts?

------------------
Gruß 

Markus

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

---

OK danke für die Tips,

ich werde erst im laufe der nächsten Tage dazu kommen alles auszuprobieren.
Ja wie ich sagte ist das Netz noch ziemlich mies. Es ging mir ja erstmal um eine grobe antastung an die materie.

Wie verfeinere ich das Netz am besten im Bereich des Kontaktes? Ich hatte es ja mit dem Befehl lesize,,, probiert aber das geht sicher irgendwie noch besser und einfacher...

bei der wegsteuerung, die ja funktionierte, habe ich das netz mal insgesamt verfeinert, also über das ganze Bauteil (ich weiß ist unpraktikabel) und wollte mal sehen, was dann passiert. Leider hat die Rechnung dann ewig gedauert und nicht konvergiert...

Habt ihr einen Tip, wie der Zylinder gut vernetzt werden könnte?

mfg

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

---

Ja danke, der Tip mit der weggesteuerten annäherung war genau richtig.

Nochmal die Frage, wie verfeinere ich (mapped meshing) am besten das Netz im Bereich des Kontaktes? Gibt es da noch andere Methoden als lesize?

mfg

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

---

ich würde bei lesize bleiben aber an der Stelle für SPACE ( sieh hilfe) einen wert setzen.

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

---

Wenn die Kontakte schon erstellt sind, gibt es (in der WB) im Mesh-Rider die Option "Kontaktelementgröße". Damit werden Kontakt- und Targetseite kontaktoptimiert vernetzt.

------------------
Gruß, A.

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