!*** TEMPERATURVERTEILUNG IN EINEM BETONQUADER MIT HEIZKANAL *** !*************** UNTER KLIMATISCHEN EINWIRKUNGEN *************** FINISH !***Das Vormodell wird beendet /CLEAR !***säubert Bildschirm /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 !***Quaderkanten, Schriften von weiß auf schwarz /REPLOT ! !********************************************************** ZEIT = 2 !*** 1 stationäre Berechnung !*** 2 transiente Berechnung !----- [ZEITLICHE ABFOLGE ----- ZEITSCHRITTE = 10 ZEITENDE = 3600 !********************************************************** BreiteAnzahl = 10 !*** Anzahl an Feldern LängeAnzahl = 10 !*** Anzahl an Feldern !----- [ABMESSUNGEN QUADER]----- QLänge = 0.10 !*** Angabe in [m] QBreite = 0.10 !*** Angabe in [m] QDicke = 0.10 !*** Angabe in [m] !----- [Höhe des HEIZSTRANGES] ----- HHÖHE = 0.08 !*** Angabe in [m] !----- [TEMPERATUREINSTELLUNGEN] ----- TEMPoben = -10 !*** Umgebungstemperatur oben TEMPunten = -10 !*** Umgebungstemperaur unten HEIZtemp = 55 !*** Heiztemperatur !----- [WÄRMEÜBERGANGSKOEFFIZIENTEN] ----- ALPHAoben = 100 ALPHAunten = 100 !********************************************************** !------------------------------ BEGINN MIT DEM PREPROCESSOR ------------------------------ /PREP7 ET,1,SOLID70 !*** für Betonquader ET,2,FLUID66 !*** für Heizstrang !--------------- MATERIALEIGENSCHAFTEN --------------- !---------- Beton [Körper 1] ---------- MP,KXX,1,2.1 !*** Wärmeleitfähigkeit MP,KYY,1,2.1 !*** Wärmeleitfähigkeit MP,KZZ,1,2.1 !*** Wärmeleitfähigkeit MP,DENS,1,2200 !*** Dichte MP,C,1,900 !*** spez. Wärmekapazität !---------- thermische Materialdaten (Heizfluid) [Körper 2] MP,KXX,2,.60 !*** Wärmeleitfähigkeit MP,KYY,2,.60 !*** Wärmeleitfähigkeit MP,KZZ,2,.60 !*** Wärmeleitfähigkeit 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.02 !*** Durchmesser AHYD = 3.14*DHYD**2/4 !*** Querschnittsfläche DL = 0.01 !*** Länge der FLUID66-Parameter R,2,DHYD,AHYD,,9.81,0,1.6*DHYD*DL RMORE,1.6*DHYD*DL,,0.001 !------------------------------------------------------------------------- !----- Die Darstellung erfolgt aus dem folgenden Blickfeld ----- !------------ und ist bis zu einer Neuangabe gültig ------------ /PNUM,NODE,1 /VIEW,,-.5,-1,.8 /VUP,,Z NPLOT /CONTOUR,ALL,13,-10,5 !***Temperaturskala-Einteilung !-------------- Berechnung der Vernetzung des Quaders... [KÖRPERS 1] --------------- N,1 N,(BreiteAnzahl+1),,QBREITE !*** Y-Achse ist die Breite des Quaders FILL NGEN,(LängeAnzahl+1),100,1,(BreiteAnzahl+1),1,(QLÄNGE/ LängeAnzahl) *IF,QDICKE,EQ,.01,THEN NGEN,2,6000,1,(LängeAnzahl*100+(BreiteAnzahl+1)),1,0,0,(QDICKE/1) *ELSEIF,QDICKE,EQ,.02,THEN NGEN,3,6000,1,(LängeAnzahl*100+(BreiteAnzahl+1)),1,0,0,(QDICKE/2) *ELSEIF,QDICKE,EQ,.03,THEN NGEN,4,6000,1,(LängeAnzahl*100+(BreiteAnzahl+1)),1,0,0,(QDICKE/3) *ELSEIF,QDICKE,EQ,.04,THEN NGEN,5,6000,1,(LängeAnzahl*100+(BreiteAnzahl+1)),1,0,0,(QDICKE/4) *ELSEIF,QDICKE,EQ,.05,THEN NGEN,6,6000,1,(LängeAnzahl*100+(BreiteAnzahl+1)),1,0,0,(QDICKE/5) *ELSEIF,QDICKE,EQ,.06,THEN NGEN,7,6000,1,(LängeAnzahl*100+(BreiteAnzahl+1)),1,0,0,(QDICKE/6) *ELSEIF,QDICKE,EQ,.07,THEN NGEN,8,6000,1,(LängeAnzahl*100+(BreiteAnzahl+1)),1,0,0,(QDICKE/7) *ELSEIF,QDICKE,EQ,.08,THEN NGEN,9,6000,1,(LängeAnzahl*100+(BreiteAnzahl+1)),1,0,0,(QDICKE/8) *ELSEIF,QDICKE,EQ,.09,THEN NGEN,10,6000,1,(LängeAnzahl*100+(BreiteAnzahl+1)),1,0,0,(QDICKE/9) *ELSEIF,QDICKE,EQ,.10,THEN NGEN,11,6000,1,(LängeAnzahl*100+(BreiteAnzahl+1)),1,0,0,(QDICKE/10) *ELSEIF,QDICKE,EQ,.11,THEN NGEN,12,6000,1,(LängeAnzahl*100+(BreiteAnzahl+1)),1,0,0,(QDICKE/11) *ELSEIF,QDICKE,EQ,.12,THEN NGEN,13,6000,1,(LängeAnzahl*100+(BreiteAnzahl+1)),1,0,0,(QDICKE/12) *ELSEIF,QDICKE,EQ,.13,THEN NGEN,14,6000,1,(LängeAnzahl*100+(BreiteAnzahl+1)),1,0,0,(QDICKE/13) *ELSEIF,QDICKE,EQ,.14,THEN NGEN,15,6000,1,(LängeAnzahl*100+(BreiteAnzahl+1)),1,0,0,(QDICKE/14) *ELSEIF,QDICKE,EQ,.15,THEN NGEN,16,6000,1,(LängeAnzahl*100+(BreiteAnzahl+1)),1,0,0,(QDICKE/15) *ELSEIF,QDICKE,EQ,.16,THEN NGEN,17,6000,1,(LängeAnzahl*100+(BreiteAnzahl+1)),1,0,0,(QDICKE/16) *ELSEIF,QDICKE,EQ,.17,THEN NGEN,18,6000,1,(LängeAnzahl*100+(BreiteAnzahl+1)),1,0,0,(QDICKE/17) *ELSEIF,QDICKE,EQ,.18,THEN NGEN,19,6000,1,(LängeAnzahl*100+(BreiteAnzahl+1)),1,0,0,(QDICKE/18) *ELSEIF,QDICKE,EQ,.19,THEN NGEN,20,6000,1,(LängeAnzahl*100+(BreiteAnzahl+1)),1,0,0,(QDICKE/19) *ELSEIF,QDICKE,EQ,.20,THEN NGEN,21,6000,1,(LängeAnzahl*100+(BreiteAnzahl+1)),1,0,0,(QDICKE/20) *ENDIF !---------------------------------------------------------- E,1,2,6002,6001,101,102,6102,6101 EGEN,BreiteAnzahl,1,1 EGEN,LängeAnzahl,100,-BreiteAnzahl *IF,QDICKE,EQ,0.01,THEN EGEN,1,6000,-6000 *ELSEIF,QDICKE,EQ,0.02,THEN EGEN,2,6000,-6000 *ELSEIF,QDICKE,EQ,0.03,THEN EGEN,3,6000,-6000 *ELSEIF,QDICKE,EQ,0.04,THEN EGEN,4,6000,-6000 *ELSEIF,QDICKE,EQ,0.05,THEN EGEN,5,6000,-6000 *ELSEIF,QDICKE,EQ,0.06,THEN EGEN,6,6000,-6000 *ELSEIF,QDICKE,EQ,0.07,THEN EGEN,7,6000,-6000 *ELSEIF,QDICKE,EQ,0.08,THEN EGEN,8,6000,-6000 *ELSEIF,QDICKE,EQ,0.09,THEN EGEN,9,6000,-6000 *ELSEIF,QDICKE,EQ,0.10,THEN EGEN,10,6000,-6000 *ELSEIF,QDICKE,EQ,0.11,THEN EGEN,11,6000,-6000 *ELSEIF,QDICKE,EQ,0.12,THEN EGEN,12,6000,-6000 *ELSEIF,QDICKE,EQ,0.13,THEN EGEN,13,6000,-6000 *ELSEIF,QDICKE,EQ,0.14,THEN EGEN,14,6000,-6000 *ELSEIF,QDICKE,EQ,0.15,THEN EGEN,15,6000,-6000 *ELSEIF,QDICKE,EQ,0.16,THEN EGEN,16,6000,-6000 *ELSEIF,QDICKE,EQ,0.17,THEN EGEN,17,6000,-6000 *ELSEIF,QDICKE,EQ,0.18,THEN EGEN,18,6000,-6000 *ELSEIF,QDICKE,EQ,0.19,THEN EGEN,19,6000,-6000 *ELSEIF,QDICKE,EQ,0.20,THEN EGEN,20,6000,-6000 *ENDIF *IF,HHÖHE,EQ,0.08,THEN X=48000 NGEN,2,130000,1,(LängeAnzahl*100+BreiteAnzahl+1),1,0,0,HHÖHE TYPE,2 MAT,2 REAL,2 /TITLE,->RANDABSTAND:%(QBREITE/BreiteAnzahl*BreiteAnzahl/2)% [m] E,130001+(BreiteAnzahl/2),130001+(BreiteAnzahl/2)+100,X+1+(BreiteAnzahl/2),X+1+(BreiteAnzahl/2)+100 EGEN,LängeAnzahl,100,-1 EPLOT *ENDIF !-------------------------------------------------------------------- !********** STATIONÄRE BERECHNUNG ********** !******************************************* *IF,ZEIT,EQ,1,THEN /SOLU D,(130001+BreiteAnzahl/2),TEMP,HeizTEMP !*** KONVEKTION OBEN NSEL,S,LOC,Z,QDICKE !***Knoten an der Oberfläche SF,ALL,CONV,ALPHAoben,TEMPoben !***Aufbringen der Konvektion !*** KONVEKTION UNTEN NSEL,S,LOC,Z,0 !***Knoten an der Unterfläche SF,ALL,CONV,ALPHAunten,TEMPunten !***Aufbringen der Konvektion ALLS SOLVE !----------------------- /POST1 /PNUM,NODE PLNS,TEMP NSEL,S,,,130001,135051 !*** Ausgabe der Heizstrang-Temperaturen PRNS,TEMP NSORT,TEMP PRNS,TEMP ALLS NSEL,S,,,1,150000 !*** Darstellung des Quaders inkl. Heizstrang /PNUM,NODE PLNS,TEMP ALLS /AUTO,1 !***BILDGRÖßE ANS FENSTER /REP,FAST !***ANPASSEN /SHOW,JPEG JPEG,QUAL,75 JPEG,ORIENT,HORIZ JPEG,COLOR,2 JPEG,TMOD,1 !***BILD WIRD GESCHOSSEN /GFILE,1200 !* /REPLOT /SHOW,CLOSE /DEVICE,VECTOR,0 *ENDIF !-------------------------------------- !********** TRANSIENTE BERECHNUNG ********** !******************************************* *IF,ZEIT,EQ,2,THEN /SOLU ANTYPE,TRANS AUTOTS,ON TIMINT,ON TINTP,,,,1 KBC,1 *DO,i,1,ZEITSCHRITTE TIME,ZEITENDE DELTIM,(ZEITENDE/ZEITSCHRITTE) OUTRES,ALL,ALL D,(130001+BreiteAnzahl/2),TEMP,HeizTEMP !*** KONVEKTION OBEN NSEL,S,LOC,Z,QDICKE !***Knoten an der Oberfläche SF,ALL,CONV,ALPHAoben,TEMPoben !***Aufbringen der Konvektion !*** KONVEKTION UNTEN NSEL,S,LOC,Z,0 !***Knoten an der Unterfläche SF,ALL,CONV,ALPHAunten,TEMPunten !***Aufbringen der Konvektion ALLS *ENDDO SOLVE !------------------ /POST1 NSEL,S,,,1,150000 /PNUM,NODE PLNS,TEMP ALLS /AUTO,1 !***BILDGRÖßE ANS FENSTER /REP,FAST !***ANPASSEN /SHOW,JPEG JPEG,QUAL,75 JPEG,ORIENT,HORIZ JPEG,COLOR,2 JPEG,TMOD,1 !***BILD WIRD GESCHOSSEN /GFILE,1200 !* *DO,i,1,ZEITSCHRITTE ZEIT=ZEITENDE/ZEITSCHRITTE*i SET,,,,,ZEIT /REPLOT *ENDDO /SHOW,CLOSE /DEVICE,VECTOR,0 *ENDIF