finish
/clear,start

/title,Enthalpietest

/prep7

radanf=1e-6   !Radius Anfang
strend=40     !Strom-Ende
pulsd=2e-6    !Pulsdauer
radwa=80       !Radiuswachstum [µm/µs]
strwa=200e+6   !Stromwachstum [A/s]
radmi=25e-6   !Radius Mittlerer
radma=50e-6   !Radius Maximum
fupusp=14     !Fußpunktspannungsabfall
te=1e-12


et,1,plane55,,,1 
mp,kxx,1,84
mp,dens,1,7860

!nichtlineare Enthalpie

mptemp,1,300,1800,1810,2500,2800,3000
mpdata,enth,1,1,9.8e+8,6.5e+9,8.7e+9,1.2e+10,6.2e+10,6.3e+10
mpplot,enth,1

!vereinfachte Enthalpie

!mptemp,1,300,2800,5000
!mpdata,enth,1,1,9.8e+8,6.2e+10,6.21e+10
!mpplot,enth,1

!Geometrie

pcirc,0,radmi,0,-90
pcirc,0,radma,0,-90
aovlap,1,2    !A1=Innen A3=Außen
aglue,all

!Vernetzung innerer Teil

esize,1e-6
mshape,0,2d
mshkey,1
amesh,1

!Vernetzung des äußeren Teils

csys,1
lsel,s,loc,x,radmi,radma
lsel,u,loc,x,radmi
lsel,u,loc,x,radma
lesize,all,,,6,6
esize
amesh,3

TUNIF,300

/solu
antype,trans
kbc,1
time,te
nsubst,1
solve


radend=strend/strwa*radwa
sw=50e-9          !Schrittweite
schritte=pulsd/sw
cnvtol,heat,1000,0.1

*do,i,1,schritte
zeit=sw*i

!Strom
strom=strwa*zeit
*if,strom,gt,strend,then
strom=strend
*endif

!Radius
rad=radanf+radwa*zeit
*if,rad,gt,radend,then
rad=radend
*endif

!Fläche
fla=3.1415*rad*rad

nsel,s,loc,y,0
nsel,r,loc,x,0,rad
sf,all,hflux,strom*fupusp/fla*1
nsel,all

time,zeit
autots,on
neqit,50
nsubst,10,200,2
solve
*enddo

finish

/post1
/edge,1,1
/plots,minm,0
/cont,,,300,200,2900

set,first
*do,i,1,schritte
plns,temp
set,next
*enddo