!aus dem Forum von franky85
!Änderungen Schätzing, 25.6.2013

finish  $/clear

*abbr,laserpolieren2,/input,laserpolieren2,txt



abma=20e-3		!Abmessung Aussen
rand=4e-3		!Rand

/prep7
et,1,mesh200,6		!Hilfselement für die Vernetzung
et,2,solid70

!----Materialvereinbarung-----
mp,kxx,1,54
mp,dens,1,7850
mp,c,1,461

!-----Erstellung der Geometrie------
rectng,rand,abma-rand,abma/2,abma-rand

type,1
esize,abma/120
amesh,all


!-----Netz um 0.005 in die Tiefe extrudieren mit 5 Elementen----

abmz=1e-3		!Abmessung in die Z-Richtung
esize,,5
type,2
mat,1
vext,all,,,,,-abmz
aclear,all			!Hilfselemente löschen

r=(4.07e9*0.12)		!Wärmestromdichte mit der Absorption multipliziert

vv=7.5e-5		!Vorschubweg    0.075mm
spu=1.5e-4		!Spurabstand

finish
/config,nres,30000

/solu
antype,trans
kbc,1
tunif,293
time,1e-12
outres,all,all
solve
eplot


Spur=1

!n: aktuelle Spur
!i: aktueller Puls

!*do,n,0,Spur
n=0

!----Schleife zum Einschalten des Laserpuls----

Pulse=127
pulse=27

*do,i,0,Pulse

time,(0.0005+(i*0.01)+(n*2.56))
!time,(0.0005+(i*0.01)+(i*2.56))
sfdele,all,all,all
local,11,1,(0.005+(i*vv)),(0.015-(2*n*spu)),0,0,0,0
wpcsys,,11
nsel,s,loc,x,0,0.00025
nsel,r,loc,z,0
sf,all,hflux,r
nsub,1
allsel
outres,all,all
solve

!-----Ausschalten des Laserpuls-----

time,(0.001+(i*0.01)+(n*2.56))
!time,(0.001+(i*0.01)+(i*2.56))
sfdele,all,all,all
nsub,9
nsub,1
allsel
outres,all,all
solve

*enddo  !end i



!/eof

/post1
/cont,,9,293,,1300
/plopts,minm,0
plns,temp

antime,pulse*2,0.1,,1,1