!!!----- BETONQUADER UNTER KLIMATISCHEN EINWIRKUNGEN -----!!!

FINISH					!***DAS VORMODELL WIRD BEENDET		
/CLEAR
!/OUTPUT,BETONQUADER,TXT
!/FILENAME,betonquader,1
/TITLE,- EIN BEHEIZTER BETONQUADER UNTER KLIMATISCHEN EINWIRKUNGEN -

!/RGB,INDEX,100,100,100,0		!
!/RGB,INDEX,80,80,80,13			!***Ändern des Hintergrundes
!/RGB,INDEX,60,60,60,14			!***von schwarz auf weiß
!/RGB,INDEX,0,0,0,15			!
!/REPLOT

!------------------------------ PARAMETER --------------------------------
QLÄNGE=0.5				!***Quaderlänge  in [m]
QBREITE=0.5				!***Quaderbreite in [m]
QHÖHE=0.2				!***Quaderhöhe   in [m]

HHÖHE=.18				!***Heizkanalhöhe in [m]
THEIZ=1000				!***Heizwassertemperatur

TUMGoben=10				!***Umgebungstemperatur OBEN
TUMGunten=5				!***Umgebungstemperatur UNTEN
ALPHAoben=100				!***Wärmeübergangskoeffizient OBEN
ALPHAunten=100				!***Wärmeübergangskoeffizient UNTEN


ZEITENDE=60*60*4			!***Solution-Dauer 100 Stunden

!-------------------------------------------------------------------------

/PREP7

ET,1,SOLID70
ET,2,FLUID116

!------------------------ MATERIALDATEN DER KÖRPER -----------------------

!---------- Beton [Körper 1] ----------

MP,DENS,1,2200				!***Dichte
MP,KXX,1,2.1				!***Wärmeleitfähigkeit
MP,C,1,920				!***spez. Wärmekapazität

!----------- thermische Materialdaten (Heizwasser)  [Körper 2] -----------
MP,KXX,2,80
MP,DENS,2,1000				!***Dichte
MP,VISC,2,1e-8				!***Viskosität
MP,HF,2,300				!***Wärmeübergangskoeffizient
MP,C,2,4000				!***spez. Wärmekapazität

!--------------------- Querschnittsdaten für FLUID66 ---------------------

DHYD=0.08				!***Durchmesser
AHYD= 3.14*DHYD**2/4			!***Querschnittsfläche
DL=0.1					!***Länge der FLUID66-Parameter

!-- Mit diesen Parametern werden die ersten real constant festgelegt --

R,2,DHYD,AHYD,,9.81,0,1.6*DHYD*DL
RMORE,1.6*DHYD*DL,,0.001


!------------------ EINGABE DER KNOTEN DES QUADERS -----------------------

N,1					!***Knoten 1 wird im Koordinatenursprung angeordnet

N,11,,QLÄNGE				!***auf einer Länge von QLÄNGE werden 11 Punkte angeordnet
FILL
NGEN,11,20,1,11,1,QBREITE/10		!***11 Knoten in x-y-Ebene,jeweils in QBREITE/10  
					!***in x-Richtung versetzt
NGEN,5,300,1,211,1,0,0,QHÖHE/4
NGEN,2,2000,1,211,1,0,0,HHÖHE		!***ANORDNUNG DER NUMMERN FÜR DEN HEIZSTRAHL. 
					!***ES MÜSSEN ANDERE SEIN,DA IM HEIZSTRANG IN HÖHE VARIABEL 
					!*** UND IM BETON UNTERSCHIEDLICHE TEMP. HERRSCHEN !!!


!----- Die geschaffene Anordnung der Punkte soll nun grafisch dargestellt werden. -----

/PNUM,NODE,1
/VIEW,,-.5,-1,.8
/VUP,,Z
NPLOT

!------------------------- Elemente des Betonquaders festlegen -------------------------

E,1,2,302,301,21,22,322,321
EGEN,10,1,1			
EGEN,10,20,-10
EGEN,4,300,-100

!---------- Der Heizwasserkanal soll durch folgende Elemente gebildet werden ----------

TYPE,2				!***Der Strömungskanal wird mit Elementen des zweiten 
MAT,2				!***Elementtyps, des zweiten Materials und der zweiten
REAL,2				!***Zusatzgruppe aufgebaut.

E,2009,2029,609,629
EGEN,8,20,-1
E,2169,2168,669,668
EGEN,2,-1,-1
E,2167,2147,667,647
EGEN,6,-20,-1
E,2047,2046,647,646
EGEN,2,-1,-1
E,2045,2065,645,646
EGEN,6,20,-1
E,2165,2164,665,664
EGEN,2,-1,-1
E,2163,2143,663,643
EGEN,8,-20,-1

EPLOT				!***Kanal grafisch darstellen

!--------------------------------- Der Lösungsabschnitt ---------------------------------

/SOLU

antype,trans
outres,all,all

SCHRITTE=20				!***zeitliche Abfolge ???????????????
deltime,ZEITENDE/SCHRITTE
TIME,ZEITENDE
SOLVE



D,2009,TEMP,THEIZ		!***Die Temperatur der Flüssigkeit beim Einströmen


!*** KONVEKTION OBEN

NSEL,S,LOC,Z,0.2			!***Knoten an der Oberfläche
SF,ALL,CONV,ALPHAoben,TUMGoben		!***Aufbringen der Konvektion


!*** KONVEKTION UNTEN

NSEL,S,LOC,Z,0				!***Knoten an der Unterfläche
SF,ALL,CONV,ALPHAunten,TUMGunten	!***Aufbringen der Konvektion

ALLS
SOLVE


!----------------------------------------------------------------------------------------

/POST1
NSEL,S,,,2001,2211
PRNS,TEMP
NSORT,TEMP
PRNS,TEMP
ALLS

NSEL,S,,,1,2211			!***Konturplot des Modells
/PNUM,NODE
PLNS,TEMP