/Prep7
emunit,mks

!Daten allgemein

alpha=10		! Wärmeübergangskoeffizient
xnetz=2000		! Vernetzung in x-richtung
ynetz=5			! Vernetzung in y-richtung
Luftdicke=0.075		! Größe des Luftraumes um das blech
breite=1		! 1 m 
dicke= 0.0005		! 0.5 mm

!Daten Luft

teta_luft=20		! Temperatur
lambda_luft=0.025	! Wärmeleitfähigkeit
rho_luft=1.2		! Dichte	
c_luft=287		! Wärmekapazität			

!Daten Aluminium

teta_alu=80
lambda_alu=200
rho_alu=2.7
c_alu=900



!Geometrie

x0= -Luftdicke
x1= 0
x2= breite
x3= breite+Luftdicke

y0=-dicke/2-Luftdicke
y1=-dicke/2 
y2=0
y3=dicke/2
y4=dicke/2+Luftdicke 

rectng,x0,x3,y0,y4		! REchteck für die Luft
rectng,x1,x2,y1,y3		! REchteck für das Blech	
aovlap,1,2			! Überlappen der beiden Flächen
adele,1,,,1			! Löschen der Fläche unterm Blech
asel,all
aglue,all			! Verkleben der Flächen


!Materialdaten zuweisen

Et,1,Plane55
Et,2,Plane55

Mp,dens,1,rho_luft
Mp,kxx,1,lambda_luft
Mp,c,1,c_luft

Mp,dens,2,rho_alu
Mp,kxx,2,lambda_alu
Mp,c,2,c_alu

asel,s,area,,3
aatt,1,1,1

asel,s,area,,2
aatt,2,1,2

asel,all 

! Vernetzung
mshkey,1		!Mapmashing
mshape,1,2d		!Dreiecke

lsel,s,loc,x,x1
lsel,r,loc,y,y1,y2
lesize,all,,,ynetz

lsel,s,loc,x,x1
lsel,r,loc,y,y2,y3
lesize,all,,,ynetz

lsel,s,loc,x,x1,x2
lsel,r,loc,y,y1
lesize,all,,,xnetz
amesh,2

mshkey,0		! Freemeshing
amesh,3 

/SOLU 

ANTYPE,4                    ! Transiente Berechnung
SOLCONTROL,ON,0             ! Weiß nicht wofür steht bei den Übungen auch immer

toffst,273.15 		    ! Umrechnung ° -> K

! Blech hat ein Startwert von 80°
nsel,s,loc,x,x1,x2
nsel,r,loc,y,y1,y3 
ic,all,temp,teta_alu

! Luft hat eine Starttemperatur von 20°
NSEL,INVE                    
ic,all,temp,teta_luft

TIME,0.5                        ! Länge der Berechnung 1/2 Stunde
AUTOTS,-1                    	! ???
DELTIM,0.001,0.001,0.025,1	! Startpunkt und Zeitschritte
KBC,1                        !???

!Heatfluxes
nsel,s,loc,x,x1
nsel,r,loc,y,y1,y3
SF,all,hflux,c_alu	! Was muss als letzter Parameter stehen? Ist die spezifische Temperatur da richtig?

nsel,s,loc,x,x2
nsel,r,loc,y,y1,y3
SF,all,hflux,c_alu

nsel,s,loc,x,x1,x2
nsel,r,loc,y,y1
SF,all,hflux,c_alu

nsel,s,loc,x,x1,x2
nsel,r,loc,y,y3
SF,all,hflux,c_alu

!Konvektion um das Blech 
nsel,s,loc,x,x1
nsel,r,loc,y,y1,y3
sf,all,conv,alpha,teta_luft
 
nsel,s,loc,x,x2
nsel,r,loc,y,y1,y3
sf,all,conv,alpha,teta_luft

nsel,s,loc,x,x1,x2
nsel,r,loc,y,y1
sf,all,conv,alpha,teta_luft

nsel,s,loc,x,x1,x2
nsel,r,loc,y,y3
sf,all,conv,alpha,teta_luft

OUTRES,ALL,ALL               ! Ausgabe 
SAVE
/STAT,SOLU                   
/REPLOT                      
APLOT                        

alls 
solve 
/POST1 
!PLVAR,2 
nsel,all 
plns,temp 
save