!------------------------
!RANDBEDINGUNGEN & RECHNUNG
!------------------------

!Erstellen der Temp-Zeitfunktion als Parameterarray
anz=12
*dim,z,array,anz
*dim,t,array,anz
z(1)=0	!wird nicht benutzt Leerschritt
t(1)=20
!1. Ausgangspunkt der Transiente
z(2)=1
t(2)=349.2
!2. Anstieg von 0 bis 2 min
z(3)=2
t(3)=444.5
!3. Anstieg von 2 bis 8 min
z(4)=8
t(4)=645.5
!4. Anstieg von 8 bis 15 min
z(5)=15
t(5)=735
!4. Anstieg von 15 bis 24 min
z(6)=24
t(6)=808.5
!5. Anstieg von 24 bis 30 min
z(7)=30
t(7)=841.8
!5. Anstieg von 30 bis 45 min
z(8)=45
t(8)=902.3
!5. Anstieg von 45 bis 56 min
z(9)=56
t(9)=935
!6. Anstieg von 56 bis 77 min
z(10)=77
t(10)=982
!6. Anstieg von 77 bis 80 min
z(11)=80
t(11)=988.4
!6. Anstieg von 80 bis 90 min
z(12)=90
t(12)=1006

!Kontrolle
*stat,z(1)
*stat,t(1)

!------------------------
!RECHNUNG
!------------------------

/prep7
tunif,20		!Umgebungstemperatur t=0
ic,all,temp,20	!initial Conditions

asel,s,,,311	!Selektion der Belastungsfläche
asel,a,,,307
asel,a,,,303
asel,a,,,299
asel,a,,,295

/solu
antype,4
deltim,30,10,30,on
outres,all,all
kbc,0			
h_f2=13435

time,1*60
sfa,all,1,hflux,h_f2
solve	
			
*get,n_3,node,node(0,0,0),temp
h1=25*(t(3)-n_3)				!Temp. conv-Anteil
h2=(t(3)+273)*(t(3)+273)*(t(3)+273)*(t(3)+273)
h3=(n_3+273)*(n_3+273)*(n_3+273)*(n_3+273)
h4=h2-h3
h5=3.969e-8*h4				!Temp. rad-Anteil
h_3=h1+h5	
h_f3=h_3-h_f2				!Nettowärmestrom-Differenz H_f
*status

deltim,30,10,30,on
time,z(3)*60	!T in Minuten

sfa,all,1,hflux,h_f3	
solve
			
*get,n_4,node,node(0,0,0),temp
h1=25*(t(4)-n_4)				!Temp. conv
h2=(t(4)+273)*(t(4)+273)*(t(4)+273)*(t(4)+273)
h3=(n_4+273)*(n_4+273)*(n_4+273)*(n_4+273)
h4=h2-h3
h5=3.969e-8*h4				!Temp. rad
h_4=h1+h5	
h_f4=h_4-h_3				!Nettowärmestrom H_f
*status

deltim,30,10,30,on
time,z(4)*60	!T in Minuten

sfa,all,1,hflux,h_f4	
solve
		
*get,n_5,node,node(0,0,0),temp
h1=25*(t(5)-n_5)				!Temp. conv
h2=(t(5)+273)*(t(5)+273)*(t(5)+273)*(t(5)+273)
h3=(n_5+273)*(n_5+273)*(n_5+273)*(n_5+273)
h4=h2-h3
h5=3.969e-8*h4				!Temp. rad
h_5=h1+h5
h_f5=h_5-h_4				!Nettowärmestrom H_f
*status

deltim,30,10,30,on
time,z(5)*60	!T in Minuten

sfa,all,1,hflux,h_f5	
solve

			
*get,n_6,node,node(0,0,0),temp
h1=25*(t(6)-n_6)				!Temp. conv
h2=(t(6)+273)*(t(6)+273)*(t(6)+273)*(t(6)+273)
h3=(n_6+273)*(n_6+273)*(n_6+273)*(n_6+273)
h4=h2-h3
h5=3.969e-8*h4				!Temp. rad
h_6=h1+h5
h_f6=h_6-h_5				!Nettowärmestrom H_f
*status

deltim,30,10,30,on
time,z(6)*60	!T in Minuten

sfa,all,1,hflux,h_f6	
solve
			
*get,n_7,node,node(0,0,0),temp
h1=25*(t(7)-n_7)				!Temp. conv
h2=(t(7)+273)*(t(7)+273)*(t(7)+273)*(t(7)+273)
h3=(n_7+273)*(n_7+273)*(n_7+273)*(n_7+273)
h4=h2-h3
h5=3.969e-8*h4				!Temp. rad
h_7=h1+h5
h_f7=h_7-h_6				!Nettowärmestrom H_f
*status

deltim,10,10,10,on
time,z(7)*60	!T in Minuten

sfa,all,1,hflux,h_f7	
solve
			
*get,n_8,node,node(0,0,0),temp
h1=25*(t(8)-n_8)				!Temp. conv
h2=(t(8)+273)*(t(8)+273)*(t(8)+273)*(t(8)+273)
h3=(n_8+273)*(n_8+273)*(n_8+273)*(n_8+273)
h4=h2-h3
h5=3.969e-8*h4				!Temp. rad
h_8=h1+h5
h_f8=h_8-h_7				!Nettowärmestrom H_f
*status


deltim,10,10,10,on
time,z(8)*60	!T in Minuten

sfa,all,1,hflux,h_f8	
solve

			
*get,n_9,node,node(0,0,0),temp
h1=25*(t(9)-n_9)				!Temp. conv
h2=(t(9)+273)*(t(9)+273)*(t(9)+273)*(t(9)+273)
h3=(n_9+273)*(n_9+273)*(n_9+273)*(n_9+273)
h4=h2-h3
h5=3.969e-8*h4				!Temp. rad
h_9=h1+h5
h_f9=h_9-h_8				!Nettowärmestrom H_f
*status

deltim,10,10,10,on
time,z(9)*60	!T in Minuten

sfa,all,1,hflux,h_f9	
solve

			
*get,n_10,node,node(0,0,0),temp
h1=25*(t(10)-n_10)				!Temp. conv
h2=(t(10)+273)*(t(10)+273)*(t(10)+273)*(t(10)+273)
h3=(n_10+273)*(n_10+273)*(n_10+273)*(n_10+273)
h4=h2-h3
h5=3.969e-8*h4				!Temp. rad
h_10=h1+h5
h_f10=h_10-h_9				!Nettowärmestrom H_f
*status

deltim,30,10,30,on
time,z(10)*60	!T in Minuten

sfa,all,1,hflux,h_f10	
solve

		
*get,n_11,node,node(0,0,0),temp
h1=25*(t(11)-n_11)				!Temp. conv
h2=(t(11)+273)*(t(11)+273)*(t(11)+273)*(t(11)+273)
h3=(n_11+273)*(n_11+273)*(n_11+273)*(n_11+273)
h4=h2-h3
h5=3.969e-8*h4				!Temp. rad
h_11=h1+h5
h_f11=h_11-h_10				!Nettowärmestrom H_f
*status

deltim,30,10,30,on
time,z(11)*60	!T in Minuten

sfa,all,1,hflux,h_f11	
solve

			
*get,n_12,node,node(0,0,0),temp
h1=25*(t(12)-n_11)				!Temp. conv
h2=(t(12)+273)*(t(12)+273)*(t(12)+273)*(t(12)+273)
h3=(n_12+273)*(n_12+273)*(n_12+273)*(n_12+273)
h4=h2-h3
h5=3.969e-8*h4				!Temp. rad
h_12=h1+h5
h_f12=h_12-h_11				!Nettowärmestrom H_f
*status

deltim,30,10,30,on
time,z(12)*60	!T in Minuten

sfa,all,1,hflux,h_f12	
solve