fini
/CLEAR,START

!#######################################################################
!#                                                                     #
!#    Tiefziehen                                                       #
!#                                                                     #
!#    Kontaktanalyse - rigid-deformable ohne Reibung                   #
!#                                                                     #
!#    Elastisch-plastisches Material                                   #
!#    (kinematische Verfestigung, plastische Anisotropie nach Hill)    #
!#                                                                     #
!#######################################################################

/filnam,nlfem8-tz
/title, Tiefziehen
!Preprozessing
/prep7

!Definition der Elementtypen

et,1,solid45
r,1 
et,2,conta173
et,3,targe170
et,4,shell143

!Realkonstanten für Zuordnung der Kontaktpartner

r,2 !Kontakt Werkstück - Matrize
r,3 !Kontakt Werkstück - Stempel 
r,4 !Kontakt Werkstück - Niederhalter

!Reibwerte für die Kontaktpaarungen

mp,mu,2,0.1
mp,mu,3,0.1


!Materialdaten

!elastisch

mp,ex,1,210000
mp,prxy,1,0.3

!plastisch

!bilineares Modell mit kinematischer Verfestigung

tb,bkin,1
tbdat,1,300,1000                     !Flie�grenze,Tangentenmodul

!Anisotropie nach Hill

tb,hill,1
tbdata,1,1.0,1.2,0.85,0.95,0.8,1.1   !Faktoren f�r Flie�spannungen

!Modellgenerierung

!Parameter

l1=30
l2=30
r1=6
r2=3
r3=5
d1=30
d2=35
d3=60
d4=70
t=3
h=1

!Rotationsachse
k,6
k,7,,l2

!Matrize
k,1,d2,0
k,2,d2,l2-r2
k,3,d2+r2,l2
k,4,d4,l2
k,5,d2+r2,l2-r2
larc,3,2,5,r2     				!Kreisbogen
cm,rundung1,line
l,1,2
l,3,4
cm,lines_matrize,line
arotat,lines_matrize,,,,,,6,7,360,4     	!Fläche durch Rotation von Linien
cm,matrize,area
asel,u,area,,all
lsel,u,line,,all

!Stempel
k,101,,l2+t+h
k,102,d1-r1,l2+t+h
k,103,d1,l2+t+r1+h
k,104,d1,l2+t+l1+h
k,105,d1-r1,l2+t+r1+h
larc,103,102,105,r1
cm,rundung2,line
l,101,102
l,103,104
cm,lines_stempel,line
arotat,lines_stempel,,,,,,6,7,360,4
cm,stempel,area
asel,u,area,,all
lsel,u,line,,all

!Niederhalter
k,201,d2,l2+t+r3+h
k,202,d2+r3,l2+t+h
k,203,d4,l2+t+h
k,204,d2+r3,l2+t+r3+h
larc,202,201,204,r3
cm,rundung3,line
l,202,203
cm,lines_niederhalter,line
arotat,lines_niederhalter,,,,,,6,7,360,4
cm,niederhalter,area
asel,u,area,,all
lsel,u,line,,all

!Ronde
k,301,,l2
k,302,d3,l2
l,301,302
cm,lines_ronde,line
arotat,lines_ronde,,,,,,6,7,360,4
cm,areas_ronde,area
vext,areas_ronde,,,,3     !Extrudieren von Flaechen zu einem Volumen
cm,ronde,volu
allsel

cmgrp,rundungen,rundung1,rundung2,rundung3

!Vernetzung Target-Flächen

n,10000,0,70   ! Pilot-Knoten Stempel
n,10001,0,80   ! Pilot-Knoten Niederhalter

mshkey,1    !mapped meshing
esize,5
lesize,rundungen,1

!Matrize

type,3
real,2
mat,2
amesh,matrize
allsel

!Stempel

type,3
real,3
mat,3
amesh,stempel

TSHAP,PILO  		! Create a pilot node
E,10000
allsel

! Niederhalter

type,3
real,4
mat,3
amesh,niederhalter

TSHAP,PILO  		! Create a pilot node
E,10001
allsel

!Vernetzung Ronde

mshkey,1
esize,5
type,4
real,1
mat,1
amesh,areas_ronde
allsel

esize,,3
type,1
mat,1
vsweep,ronde
!Löschen shell-Elemente
aclear,areas_ronde
!Vernetzung Kontaktelemente

!untere Seite der Ronde

type,2
real,2
mat,2
esel,s,type,,1
nsle,s
nsel,r,loc,y,l2-0.0001,l2+(t/6)
esln,s
esurf
allsel

!obere Seite der Ronde

type,2
real,3
mat,3
esel,s,type,,1
nsle,s
nsel,r,loc,y,l2+((5/6)*t),l2+t+0.0001
esln,s
esurf
real,4
esurf
allsel

!Drehen der falsch orientierten Elementkoordinatensysteme

asel,s,area,,matrize
esla,s
esurf,,reverse
allsel

fini

/solu
RESCONTROL,,all,last

!1. Lastschritt - Niederhalter anlegen

!Randbedingungen

!Fixierung der Mitte der Ronde & der Pilotknoten in x- und z-Richtung

nsel,s,loc,x,0
nsel,r,loc,z,0
d,all,ux
d,all,uz
allsel

!Fixierung des Randes der Ronde in y-Richtung

csys,5
esel,s,type,,1
nsle,s
nsel,r,loc,x,d3
nsel,r,loc,z,l2
d,all,uy
csys,0
allsel

!Bewegung Niederhalter

d,10001,uy,-(h-0.0001)

!Lösungsoptionen
antype,static,new
nlgeom,on
autots,on
time,5
nsubst,1,5,1
solve

!2. Lastschritt � Stempel anlegen

!Randbedingungen

!Stempelbewegung

d,10000,uy,-(h+0.0001)

!Lösungsoptionen

nlgeom,on
autots,on
time,9
nsubst,1,5,1
solve

!3. Fixierung des Randes der Ronde in y-Richtung  aufheben

!Randbedingungen

csys,5
esel,s,type,,1
nsle,s
nsel,r,loc,x,d3
nsel,r,loc,z,l2
ddele,all
csys,0
allsel

!Lösungsoptionen

nlgeom,on
autots,on
time,10
nsubst,1,5,1
solve

!4. Lastschritt � Stempelbewegung � Teil I

!Randbedingungen

!Stempelbewegung

d,10000,uy,-(15+h)


!Lösungsoptionen

nlgeom,on
autots,on
time,310
nsubst,750,2000,50
outres,all,all
solve

!5. Lastschritt � Stempelbewegung � Teil II

!Randbedingungen

!Stempelbewegung

d,10000,uy,-(25+h)

!Lösungsoptionen

nlgeom,on
autots,on
time,510
nsubst,2500,5000,50
outres,all,all
solve