!aus dem Forum von MurinSplint
!Werkzeugheizung

finish   $/clear
*abbr,murinsplint,/input,murinsplint,txt

!Definieren des Berechnungsmodells
!----------------------------------------------------------------------------------------

breit=0.4
hoch=0.4
lang=1  
tumg=20

/Prep7
   et,1,solid70
et,1,solid226,10
   mp,dens,1,7800				! Dichte kg/m^3
   mp,kxx,1,48					! Wärmeleitkoeffizient W/(m^2*K)
   mp,c,1,452					! spez. Wärmekapazität J/(kg*K)
   
!----Geometrie----------
   Block,0,breit,0,hoch,0,lang
/view,,1,1,1

!----Vernetzung----------
   esize,0.05
   vmesh,all

!----Randbedingungen------
nsel,s,ext
nsel,u,loc,z,0
sf,all,conv,20,tumg
sf,all,conv,30,tumg
sf,all,conv,40,tumg

!-----Lasten--------
   nsel,s,loc,z,0
   sf,all,hflux,18750				! Wärmequelle an t=0, Stirnfläche
 
  
allsel  
!-------Lösung---------------------------------------------------------------------------------
/solu
!antype,trans
solve
!-------Auswertung------
/post1
/plopts,minm,0
plns,temp

/eof

!----Blick in Innere------
nsel,s,loc,x,0,breit/2
esln,s,1
plns,temp



! Berechnungsschleife Sensor1,Sensor2,Waermequelle
!----------------------------------------------------------------------------------------
 
*DO,Vergleich,1,1000,1
   
  
   TempMAX=90
   *get,Sensor1
   *get,Sensor2
   *get,TempMAX
 
   nsel,s,loc,z
   d,all,hflux,18750				! Wärmequelle an t=0, Stirnfläche

   				
  *IF,Sensor1,LT,TempMAX,Heizen,THEN		! Heizzyklus an

   nsel,s,loc,z
   sf,all,hflux,18750

  *ELSE,Sensor1,GT,TempMAX,nichtHeizen		! Heizzyklus aus

   nsel,s,loc,z
   sf,all,hflux,0

  *ELSEIF,Sensor1,EQ,Sensor2,90			! Berechnungsstop Temperatur erreicht

   STOP
  
  *ENDIF
  *ENDDO

!----------------------------------------------------------------------------------------