finish
/clear,start

/Title,Einzelentladungskrater bei der Mikrofunkenerosion

/prep7

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!Parameter

we=60e-6 			![J]
d=18e-6				![m]
te=0.4e-6				![s]
te2=100

Qq=0.12*we/(3.1416*(d/2)*(d/2)*te)	!Wärmestromdichte [W/m²]

Tu=20					!Umgebungstemperatur
n=50

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!Definition des Elemtetyps

ET,1,PLANE13 

keyopt,1,1,4       !bezieht therm. und mech.Freiheitsgrade mit ein

!!!!temp.abhängige Materialwerte!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

mp,nuxy,1,0.361
mp,dens,1,7750

MPTEMP,1,20,100,200,300,400,500
mp,kxx,1,15.3,,,,,23
mp,c,1,460
mp,alpx,1,9.54e-6,10.8e-6,10.8e-6,11.3e-6,11.7e-6,12.0e-6
mp,ex,1,196e+9,191e+9,185e+9,179e+9


!!!!!plastische Verformungsoption!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

tb,BKIN,1,6,,1     ! Bilinear Kinematic Hardening für 6 Temp.

TBTEMP,20
TBDATA,1,1280e+6,6e+8
TBTEMP,100
TBDATA,1,1200e+6,6.5e+8
TBTEMP,200
TBDATA,1,1130e+6,7e+8
TBTEMP,300
TBDATA,1,1050e+6,7.5e+8
TBTEMP,400
TBDATA,1,1000e+6,6.5e+8
TBTEMP,500
TBDATA,1,880e+6,6e+8

TBLIST,BKIN,1
/XRANGE,0,0.02
TBPLOT,BKIN,1

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!Geometrie

rectng,-2*d,6*d,0,3*d
rectng,,4*d,0,d
AOVLAP,1,2

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!Meshing Parameter

LESIZE,5,,,4*n
LESIZE,8,,,n

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!Meshing

MSHAPE,0,2D 
MSHKEY,1
AMESH,2

MSHAPE,1,2D 
MSHKEY,0
AMESH,3

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!Lastschritte...

finish

/SOLU

ANTYPE,TRANS,NEW	!Berechnung wird transient durchgeführt
TINTP,,,,1		!Zeit-Integrationsparameter (Methode nach Euler Rückwärtseinsetzen)
outpr,all,all
OUTRES,all,all
nlgeom,on
autots,on

lsel,s,,,9
lsel,a,,,4
lsel,a,,,3
lsel,a,,,2
lsel,a,,,10
nsll,s,1
d,all,temp,TU
alls

lsel,s,,,3
nsll,s,1
d,all,uy
d,all,ux
alls

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!Lastschritt 1 (Zimmertemp)
! 1.Lastschritt-Raumtemperatur

TIMINT,OFF		!Ausstellen der transienten Efffekte

time,te		!Zeit am Ende dieses Lastfalls		
D,ALL,TEMP,TU		!Zimmertemp
alls

solve

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!Lastschritte: bewegte Wärmequelle
!2.Lastschritt-Aufwärmen

TIMINT,on
TIME,te+te1

ddel,all,temp
NSEL,S,LOC,Y,0  
NSEL,r,LOC,X,d+d/2,2*d+d/2
SF,all,hflux,Qq

neqit,20
deltim,te1/10,te1/20,te1/5
kbc,1
alls

Solve
   
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

!3.Lastschritt-Abkühlen

/solu


time,te+te1+te2

ESEL,S,LIVE        ! Selects all active elements
NSEL,S,LOC,Y,0  
NSEL,r,LOC,X,d+d/2,2*d+d/2
SFDEL,all,HFLUX
alls

NSEL,S,LOC,Y,0  
NSEL,r,LOC,X,d+d/2,2*d+d/2
SF,all,conv,1500,20
alls

neqit,200
deltime,te2/10,/te2/20,te2/5
kbc,1
alls
solve

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

finish