FINISH
/CLEAR,START
/ERASE

/COM,   
/COM,Preferences for GUI filtering have been set to display:
/COM,  Thermal  


*create,plot,mac
/POST1  
/EFACET,4
PLNSOL, TEMP,, 0
*end

*ABBR,PLOT,plot

*set,l,200
*set,h,50
*set,b,100
ANZAHL_STUNDEN=12
Anfangstemperatur_Kiste=10
Anfangstemperatur_Eis=-2
Anfangstemperatur_Viecher=25
Anfangstemperatur_Luft=15
Prozent_Viecher=40
Prozent_Eis=33
Prozent_Luft=100-(Prozent_Eis+Prozent_Viecher)

!Modelling
/PREP7  
 
!Luft
et,1,solid5 
KEYOPT,1,1,8
KEYOPT,1,3,0
KEYOPT,1,5,0
mp,kxx,1,400
mp,dens,1,2e3
mp,c,1,10
 
BLOCK,0,l,0,h,0,b, 

!VIEH
et,2,solid5 
KEYOPT,2,1,8
KEYOPT,2,3,0
KEYOPT,2,5,0
mp,kxx,2,0.9
mp,dens,2,2e3
mp,c,2,10

!Eis
et,3,solid5 
KEYOPT,3,1,8
KEYOPT,3,3,0
KEYOPT,3,5,0
mp,kxx,3,5
mp,dens,3,2e3
mp,c,3,10

allsel,all
   
LESIZE,all,10

TYPE,1   
MAT,1
REAL,   
ESYS,0   
SECNUM, 

MSHAPE,0,3D 
MSHKEY,1

VMESH,all  

/VIEW,1,1,1,1   
/ANG,1  
/REP,FAST   
 


!Verteilung der einzelnen Inhalte
*GET,anzahl_elemente,elem,,count				
Anzahl_Viecher=nint(Prozent_Viecher/100*anzahl_elemente)
Anzahl_Eis=nint(Prozent_Eis/100*anzahl_elemente)
Anzahl_Luft=anzahl_elemente-Anzahl_Eis-Anzahl_Viecher
*del,elementel
*dim,elemente,array,anzahl_elemente,,,,,
*do,i,1,anzahl_elemente,1
	elemente(i)=1
*enddo
*do,i,1,Anzahl_Viecher,1
	frage=1
	*DOWHILE, frage
		Elementnr_Viecher=nint(rand(1,anzahl_elemente))
		*if,elemente(Elementnr_Viecher),eq,1,then
			frage=0
		*endif
	*enddo	
	elemente(Elementnr_Viecher)=2
	ESEL,S, , ,Elementnr_Viecher
	emodif,all,mat,2
*enddo
*do,i,1,Anzahl_Eis,1
	frage=1
	*DOWHILE, frage
		Elementnr_Eis=nint(rand(1,anzahl_elemente))
		*if,elemente(Elementnr_Eis),eq,1,then
			frage=0
		*endif
	*enddo	
	elemente(Elementnr_Eis)=3
	ESEL,S, , ,Elementnr_Eis
	emodif,all,mat,3
*enddo

FINISH  


!Solution
/SOLU

ANTYPE,trans

!konstante Temperatur für alle am Anfang dauerhaft
allsel,all
tunif,Anfangstemperatur_Luft

allsel,all

!VIECHER
ESEL,S,MAT,,2   
NSLE,S  
ic,all,temp,Anfangstemperatur_Viecher
allsel,all

!EIS
ESEL,S,MAT,,3
NSLE,S  
ic,all,temp,Anfangstemperatur_Eis
allsel,all

allsel,all


allsel,all

kbc,1

rechnen=60*60*ANZAHL_STUNDEN
time,rechnen

nsubsteps,100,100,10

autots,on

outresults,all,all

solve
  
FINISH 

/POST1  
 
/EFACET,1   
PLNSOL, TEMP,, 0

finish