finish
/clear

! Variablen
! Session Name
name = 'testlauf5'
! Dauer der Simulation 
t1 = 1000
! Konvektion
convval = 200
! Maximum Zeitintervall
maxtime = 10
! Anzahl der Materialien
m1 = 3

/config,nres,100000

/filname,name

/show,name,grp
/show,term

/prep7

! Geometrie
sphere,0,11.5,0,360
sphere,0,11.5,90,180
sphere,0,11.5,180,270
sphere,0,11.5,270,360
numcmp,all

block,0,17,-17,17,0,-17
block,0,17,-17,17,0,17
block,0,-17,-17,17,0,-17
block,0,-17,-17,17,0,17
vovlap,all
numcmp,all

cylind,2,3,0,40,0,360
vsbv,all,13,,,

cylind,2,3,0,40,0,360
vglue,all
numcmp,all

! Element Attribute

! Zucker
vsel,s,volu,,1,7
vsel,a,volu,,9
vsel,a,volu,,11,12
vsel,a,volu,,17,18
vatt,1,1,1

! Alu
vsel,s,volu,,8
vsel,a,volu,,10
vsel,a,volu,,13,16
vatt,2,1,2

! Stiel
vsel,s,volu,,19
vatt,3,1,3

! Vernetzen
*do,i,1,m1,1
et,i,200,10
*enddo
allsel
lesize,all,3	
mshape,1,3D
mshkey, 0
vmesh, all

! Materialdaten Zucker
! Wärmeleitfähigkeit siehe Graph Materialparameter Datei
mptemp,1,20,40,60,80,100,120
mpdata,kxx,1,1,0.28E3,0.29E3,0.30E3,0.31E3,0.314E3,0.377E3
mptemp,7,140
mpdata,kxx,1,7,0.368E3
! Spezifische Wäme aus Physical properties of sugars and their solutions S.73 Table 1/6
! Polynom ersten Grades
mp,c,1,1165.1E6,3.494E6
! Dichte von Saccharose in Abhängigkeit der Temperatur
! Polynom ersten Grades
mp,dens,1,1.527E-6,-0.5643E-9
mptemp,1,20,30,39,43,50,125
!mpdata,ex,1,1,900E3,850E3,800E3,50E3,40E3,30E3
mpdata,ex,1,1,800E3,800E3,800E3,50E3,40E3,30E3
mp,nuxy,1,0.3
mptemp
mptemp,1,20,40,43,125
mpdata,alpx,1,1,2.5E-4,2.5E-4,5E-4,5E-4
mp,reft,1,25

!tb,biso,1,6
!tbtemp,20,1
!tbdata,1,450E3,50E3
!tbtemp,30,2
!tbdata,1,425E3,50E3
!tbtemp,39,3
!tbdata,1,400E3,50E3
!tbtemp,43,4
!tbtemp,50,5
!tbdata,1,4E3,1E3
!tbtemp,125,6
!tbdata,1,3E3,1E3

tb,biso,1,6
tbtemp,20,1
tbdata,1,50E3,50E3
tbtemp,30,2
tbdata,1,50E3,50E3
tbtemp,39,3
tbdata,1,50E3,50E3
tbtemp,43,4
tbdata,1,2E3,1E3
tbtemp,50,5
tbdata,1,1E3,1E3
tbtemp,125,6
tbdata,1,1E3,1E3

! Materialdaten Alu
mp,reft,2,20
mp,kxx,2,170E3
mp,c,2,896E6
mp,dens,2,2.7E-6
mp,ex,2,70000E3
mp,nuxy,2,0.35
mp,alpx,2,23.4E-6

! Materialdaten Stiel
mp,reft,3,20
mp,kxx,3,0.24E3 
mp,c,3,1000E6
mp,dens,3,1.4E-6
mp,ex,3,3000E3 
mp,nuxy,3,0.3
mp,alpx,3,7E-5

! Kontrolle der Materialdaten
/show,file
*do,i,1,m1,1
mpplot,kxx,i
mpplot,c,i
mpplot,dens,i
mpplot,ex,i
mpplot,alpx,i
mpplot,nuxy,i
*enddo
tbplot
/show,term

! Thermische Analyse
*do,i,1,m1,1
et,i,70,0,0,0,0,0,0,0
*enddo

finish

! Solution
/solu
allsel
antype,trans
toffst,273
tunif,25

! Randbedingungen
allsel
! Starttemperatur Alu + Stiel
vsel,s,volu,,8
vsel,a,volu,,10
vsel,a,volu,,13,16
vsel,a,volu,,19
nslv,s,1
ic,all,temp,50
! Starttemperatur Zuckermasse
asel,s,area,,1,5
asel,a,area,,7
asel,a,area,,24,26
asel,a,area,,28,30
asel,a,area,,33,35
asel,a,area,,50
asel,a,area,,54
asel,a,area,,61
asel,a,area,,64
asel,a,area,,78,81
nsla,s,1
sf,all,conv,convval,25
allsel
vsel,s,mat,,1
nslv,s,1
ic,all,temp,125

allsel
timint,on
kbc,1			
trnopt,full
autots,on
deltim,1,1,maxtime,on
outres,all,all
time,t1
solve

finish

! Grafiken erstellen
/show,file

/post26
/title,Temp.-Zeit.-Diagramm
nsol,2,node(0,0,0),temp,,Kern-Mitte
nsol,3,node(0,11,0),temp,,Kern-Aussen
nsol,4,node(0,17,0),temp,,Aussen-konv.-Uebergang
/grid,1
plvar,2,3,4
finish

/post1

esel,s,mat,,1,3
! Ansicht: Q-slized Z buffer 
/type,1,8
! Ansicht normal zu WP
/cplane,1
/view,1,-1,0,0
/ang,1,90
! WP um 90° um y-Achse drehen
wprota,,,90

set,near,,,,10
/contour,1,9,25,13
/title,Temperatur 10s
plnsol,temp,,0,1

set,near,,,,20
/contour,1,9,25,13
/title,Temperatur 20s
plnsol,temp,,0,1

set,near,,,,50
/contour,1,9,25,13
/title,Temperatur 50s
plnsol,temp,,0,1

set,near,,,,100
/contour,1,9,25,12
/title,Temperatur 100s
plnsol,temp,,0,1

set,near,,,,200
/contour,1,9,25,10
/title,Temperatur 200s
plnsol,temp,,0,1

set,near,,,,300
/contour,1,9,25,9
/title,Temperatur 300s
plnsol,temp,,0,1

set,near,,,,400	
/contour,1,9,25,7
/title,Temperatur 400s
plnsol,temp,,0,1

set,near,,,,500
/contour,1,9,25,7
/title,Temperatur 500s
plnsol,temp,,0,1

set,near,,,,600
/contour,1,9,25,6
/title,Temperatur 600s
plnsol,temp,,0,1

set,near,,,,700
/contour,1,9,25,5
/title,Temperatur 700s
plnsol,temp,,0,1

set,near,,,,800
/contour,1,9,25,4
/title,Temperatur 800s
plnsol,temp,,0,1

set,near,,,,900
/contour,1,9,25,3
/title,Temperatur 900s
plnsol,temp,,0,1

set,near,,,,1000
/contour,1,9,25,2
/title,Temperatur 1000s
plnsol,temp,,0,1

/show,term

finish

save,name,db

/EOF

! Mechanische Analyse
/prep7

lsclear,all

!*do,i,1,m1,1
!et,i,65
!*enddo
!et,1,65
!et,2,65
!et,3,65
etchg,tts

allsel

finish

! Solution
/solu

!vsel,s,volu,,8
!vsel,a,volu,,10
!vsel,a,volu,,13,16
!nslv,s,1
!d,all,ux
!d,all,uy
!d,all,uz

allsel 
antype,static
toffst,273
ldread,temp,,,900,,name,rth
outres,all,all
solve
finish

! Grafiken erstellen
/post1
/contour,1,12,auto
/title,Displacement 900s
plnsol,u,sum,2,1.0
/title,Stress 900s
plnsol,s,eqv,2,1.0
/title,elastic Strain 900s
plnsol,epel,eqv,2,1.0
/title,plastic Strain 900s
plnsol,eppl,eqv,2,1.0
finish