finish				! Beenden von Prozessor bzw. laufender Prozesse
/clear				! Löscht Datenspeicher

/prep7				! Starten von Preprocessor


!+++++++++++++++++ Eingabe von Geometriedaten, Temperaturen und Materialeigenschaften +++++++++++++++++++++++++++++

	Pi=3.141692654
dgesamt=8e-3			! Gesamtdurchmesser von Spiegel (m)

LBolzen=30e-3			! Höhe der Versteifung (m)

elementgroesse=0.5e-3		! elementgroesseerenz zweier übereinanderliegender Elemente
				! Taucht auf bei der Unterteilung der Linien zur Steuerung der Vernetzung

RadiusA=5e-3			! Definiert Lage und Durchmesser des projezierten Rings
RadiusB=6e-3
	ringfläche=(Pi*((RadiusB)*(RadiusB)-(RadiusA)*(RadiusA)))
	ringbreite=RadiusB-RadiusA

	bezugstemp=20			! Bezugstemperatur für Wärmeausdehnung (Grad Celsius)
	kuehltemp=20			! Kühltemperatur des angekoppelten Kühlkörpers (Grad Celsius)

heizleistung=15			! opt. Heizleistung in Watt
	optintensitaet=(heizleistung/ringfläche)	! In Bolzen absorbierte Wärmestromdichte (W/m^2)
				

! Materialeigenschaften von Versteifung (Material 1)
	m1kxx=39.6			! Thermische Leitfähigkeit (W/(m*K))
	m1ex=400e9         		! Elastizitätsmodul (N/m^2)
	m1nuxy=0.26			! Poisson, Querkontraktionszahl
	m1alpx=11.5e-6			! Thermischer Ausdehnungskoeffizient (1/K)
	m1dens=7810			! Dichte (kg/m^3)
	m1kapa=461			! Wärmekapazität in (pJ/kgK)


! ++++++++++++++++++++++ Beginn des Hauptprogramms ++++++++++++++++++++++++++++++++++++++
! Ab hier müssen keine Parameter mehr eingestellt werden

! Festlegen der Elementtypen

et,1,plane55			! Elementtyp 1, 2-D Thermal-Solid-Element Plane 55
keyopt,1,3,1			! Keyoption 3 von Element 1: axisymmetric (d.h. radialsymmetrischer Aufbau)


! Materialdefinition des Bolzens (Material 1)
	mp,kxx,1,m1kxx			! Thermische Leitfähigkeit (W/(m*K))
	mp,reft,1,bezugstemp		! Referenztemperatur für thermische Ausdehnung
	mp,ex,1,m1ex         		! Elastizitätsmodul (N/m^2)
	mp,nuxy,1,m1nuxy		! Poisson, Querkontraktionszahl
	mp,alpx,1,m1alpx		! Thermischer Ausdehnungskoeffizient (1/K)
	mp,dens,1,m1dens
	mp,C,1,m1kapa			! Spezif. Wärmekapazität


! Setzen der Keypoints
	k,1,0,0				! Keypoint Nr 1 an Position x=0,y=0
	k,2,(dgesamt/2),0		! Keypoint Nr 2 an Position x=Dgesamt/2, y=0


	kgen,2,1,2,1,0,LBolzen,0,10	! Kopiere Keypoints 1 bis 2 mit dy=LBolzen und benenne sie mit Nr 11 bis 12


/pnum,kp,1			! Einschalten der Keypointnummerierung
kplot				! Zeige Keypoints im Display



l,1,2		!Horizontallinien
l,11,12

l,1,11		!Vertikallinien
l,2,12


a,1,2,12,11	! Flächen


/pnum,line,0
/pnum,area,1
aplot


! Im Folgenden werden die horizontalen Linien selektiert und unterteilt, um die Vernetzung zu steuern

teilera=(dgesamt/2)/(elementgroesse)	! elementgroesse ist angestrebte Höhe und Breite aller Elemente
lsel,s,line,,1,2			! Selektiere Linie 1 bis 4
lesize,all,,,nint(teilera)		! Unterteile selektierte Linien in teilera-Teile


! Im Folgenden werden die vertikalen Linien selektiert und unterteilt, um die Vernetzung zu steuern

teilerc=LBolzen/elementgroesse
lsel,s,line,,3,4
lesize,all,,,nint(teilerc)


type,1
mat,1
amesh,1


outpr,basic,1


finish


! +++++++++++++++++++++++ Berechnung des Temperaturprofils ++++++++++++++++++++++++

/solu				! Start des Solution Prozessors

allsel				! Alles selektieren

antype,trans			! Analyseart: Statische Analyse

deltim,0.03,0.01,1.0

time,1

/solu

/pnum,elem,1
eplot

kbc,1

esel,s,,,477,478

nsle,s,,all
nplot

	sf,all,hflux,optintensitaet	! HGEN benötigt eine Leistung pro Volumen

time,5

esla,s
/pbf,hgen,,1			! Zeige Bodyforce-Load Hgen
eplot
allsel				! Alles selektieren


/pbc,temp,,1			! Zeige Boundary Condition Temperatur
nplot				! Zeige Linien auf Display

allsel

solve				! Startet mit der Lösung des Problems

finish				! Beendet Solution-Processor

/post1				! Startet Post-Processor
/triad,lbot			! Globales Koodinatensystem in Ecke unten Links zeigen
plnsol,temp			! Zeige Lösung Temperatur

finish				! Beendet Post-Processor

/eof