finish
/clear,start

/Title,Einzelentladungskrater bei der Mikrofunkenerosion

/prep7

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

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


Qq=0.12*we/(3.1416*(d/2)*(d/2)*te1)	!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!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

MPTEMP,1,20,90,200,315,425,500
mp,kxx,1,7.2,7.5,8.5,9.8,11.0,12.1
mp,c,1,560,570,600,650,700,740
mp,alpx,1,8.8e-6,9.0e-6,9.3e-6,9.6e-6,9.8e-6,10e-6
mp,ex,1,113e+9,107.35e+9,101.7e+9,92.66e+9,84.75e+9,70.06e+9

mp,nuxy,1,0.361
mp,dens,1,4450

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

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

TBTEMP,20
TBDATA,1,900e+6,110e+9
TBTEMP,90
TBDATA,1,810e+6,102e+9
TBTEMP,200
TBDATA,1,630e+6,94e+9
TBTEMP,315
TBDATA,1,580e+6,90e+9
TBTEMP,425
TBDATA,1,540e+6,83e+9
TBTEMP,500
TBDATA,1,500e+6,78e+9

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