fini
/clear
*abbr,hammel1,/input,hammel1,txt


!-------
/prep7
!------------------------
!MATERIALGENERIERUNG STAHL
!------------------------

et,1,solid70


mptemp,1,20,100,115,150,200,264
mptemp,7,314,367,414,434,570,600
mptemp,13,666,689,706,712,720,725
mptemp,19,729,735,743,760,779,800
mptemp,25,900,1200
mpdata,c,1,1,439.80,487.62,494.92,510.42,529.76,552.13
mpdata,c,1,7,569.82,590.80,612.89,623.64,727.66,759.92
mpdata,c,1,13,846.58,931.35,1072.31,1166.08,1388.33,1666.15
mpdata,c,1,19,2110.67,5000,2030,1159.48,916.25,803.26
mpdata,c,1,25,650.44,650



mptemp,1,20,100,115,150,200,264
mptemp,7,314,367,414,434,570,600
mptemp,13,666,689,706,712,720,725
mptemp,19,729,735,743,760,779,800
mptemp,25,900,1200
mpdata,kxx,1,1,53.33,50.67,50.17,49.01,47.34,45.21
mpdata,kxx,1,7,43.54,41.78,40.21,39.55,35.02,34.02
mpdata,kxx,1,13,31.82,31.06,30.49,30.29,30.02,29.86
mpdata,kxx,1,19,29.72,29.52,29.26,28.69,28.06,27.36
mpdata,kxx,1,25,27.3,27.3


mp,dens,1,7850

mp,emis,1,.70

mppl,c,1
mppl,kxx,1
mppl,dens,1
mppl,emis

!------------------------
!MATERIALGENERIERUNG BETON
!------------------------
et,2,solid70


mptemp,1,20,100,115,150,200,264
mptemp,7,314,367,414,434,570,600
mptemp,13,666,689,706,712,720,725
mptemp,19,729,735,743,760,779,800
mptemp,25,900,1200
mpdata,c,2,1,901.10,942.80,2020.00,966.65,988.80,1014.67
mpdata,c,2,7,1032.95,1050.46,1064.39,1069.87,1099.87,1104.80
mpdata,c,2,13,1113.48,1115.81,1117.30,1117.78,1118.38,1118.74
mpdata,c,2,19,1119.01,1119.39,1119.87,1120.74,1121.47,1122.00
mpdata,c,2,25,1120.40,1074.80



mptemp,1,20,100,115,150,200,264
mptemp,7,314,367,414,434,570,600
mptemp,13,666,689,706,712,720,725
mptemp,19,729,735,743,760,779,800
mptemp,25,900,1200
mpdata,kxx,2,1,1.33,1.23,1.21,1.17,1.11,1.04
mpdata,kxx,2,7,0.99,0.94,0.89,0.88,0.77,0.75
mpdata,kxx,2,13,0.71,0.69,0.68,0.68,0.68,0.67
mpdata,kxx,2,19,0.67,0.67,0.66,0.66,0.65,0.64
mpdata,kxx,2,25,0.60,0.55



mptemp,1,20,100,115,150,200,264
mptemp,7,314,367,414,434,570,600
mptemp,13,666,689,706,712,720,725
mptemp,19,729,735,743,760,779,800
mptemp,25,900,1200
mpdata,dens,2,1,2353.8,2353,2352.85,2352.5,2352,2351.36
mpdata,dens,2,7,2350.86,2350.33,2349.86,2349.66,2348.3,2348
mpdata,dens,2,13,2347.34,2347.11,2346.94,2346.88,2346.8,2346.75
mpdata,dens,2,19,2346.71,2346.65,2346.57,2346.4,2346.21,2346
mpdata,dens,2,25,2345,2342



mp,emis,2,.70

mppl,c,2
mppl,kxx,2		!Achtung! et=1 (type) aber mat=2
mppl,dens,2




!------Geometrie----
breit2=-95e-3  !Breite Beton
breit1=-100e-3	!Breite Stahl
hoch1=50e-3  !Höhe Balken
hoch2=100e-3
hoch3=155e-3
lang=50e-3    !Länge Balken  (betrachtete Länge)


block,0,lang,0,hoch1,0,breit2  !Symmetrie in x		!1
block,0,lang,0,hoch1,breit2,breit1  !Symmetrie in x 	!2
block,0,lang,hoch1,hoch2,0,breit2  !Symmetrie in x	!3
block,0,lang,hoch1,hoch2,breit2,breit1  !Symmetrie in x	!4
block,0,lang,hoch2,hoch3,0,breit2  !Symmetrie in x	!5
block,0,lang,hoch2,hoch3,breit2,breit1  !Symmetrie in x	!6

!------Vernetzung-----
nummrg,all,,,,low

ESIZE,25e-3,0,  			!STAHL!
TYPE,1   
MAT,1
REAL,1
ESYS,0   
SECNUM, 


vsweep,2
vsweep,4
vsweep,6

ESIZE,25e-3,0,  			!BETON!
TYPE,1   
MAT,2
REAL,2
ESYS,0   
SECNUM, 


vsweep,1
vsweep,3
vsweep,5


allsel

!------Randbedingungen-----
tunif,20
ic,all,temp,20

allsel
!--------Lösung----

/solu
endzeit=5400			!90 min
schritte=90
antype,trans
kbc,1

outres,all,all
*do,i,0,schritte
allsel
sfedele,all,all,all

time,endzeit/schritte*i+1e-12
nsubst,3,3,3

nsel,s,loc,x,lang	!Selektion Belastungsfläche
esln,s,0

!GEREIERUNG DER BEASTUNG
!SCHLEIFE ZUR BESTIMMUNG DER OBERFLÄCHENTEMPERATUR
!NOTWENDIG ZUR BESTIMMUNG DES KORREKTEN WÄRMESTROMS JE LASTSCHRITT

*do,j,7,8,1

n_1=nelem(j,2)		!ACHTUNG face 3 --> Knutenidentnummern: 5,6,7,8
n_2=nelem(j,3)
n_3=nelem(j,6)
n_4=nelem(j,7)
n_a1=temp(n_1)
n_a2=temp(n_2)
n_a3=temp(n_3)
n_a4=temp(n_4)

n_a=(n_a1+n_a2+n_a3+n_a4)/4		!Temperaturmittelwert

wsd=((25*((20+345*(log10(8*(i)+1)))-n_a))+(0.00000003969*(((20+345*(log10(8*(i)+1)))+273)**4-(n_a+273)**4)))
sfe,j,3,hflux,,wsd
*enddo


*do,j,31,38,1

n_1=nelem(j,1)		!ACHTUNG face 3 --> Knutenidentnummern: 1,4,5,8
n_2=nelem(j,4)
n_3=nelem(j,5)
n_4=nelem(j,8)
n_a1=temp(n_1)
n_a2=temp(n_2)
n_a3=temp(n_3)
n_a4=temp(n_4)

n_a=(n_a1+n_a2+n_a3+n_a4)/4		!Temperaturmittelwert

wsd=((25*((20+345*(log10(8*(i)+1)))-n_a))+(0.00000003969*(((20+345*(log10(8*(i)+1)))+273)**4-(n_a+273)**4)))
sfe,j,5,hflux,,wsd
*enddo

allsel
solve
sfedele,all,all,all

*enddo

/vsca,,,1			!Darstellung Belastung
/psf,hflux,,3,1,on
eplot


allsel
!------Auswertung------
/post1
plns,temp