Finish /clear /CONFIG,NRES,5000 ! ___________________________________________________________________________________ ! ! Definition von Parametern ! ___________________________________________________________________________________ *set,d1,5e-3 ! Durchmesser von Topfboden ; d1=5mm *set,d2,2e-3 ! Durchmesser von Topfseite *set,dc,3e-3 ! Durchmesser von Keramikplatte *set,h,50e-3 ! Höhe von Topf *set,b,120e-3 ! Breite von Topf *set,l,160e-3 ! Breite von Herdplatte TCC=2 !zugehörigen Kontakteigenschaften TCC (thermische Kontaktleitfähigkeit) zeitende=1000 ETWASSER=1 ! Elementtype für Wasser, ETC ETTOPF=2 ETSCHICHT=3 ETKERAMIK=4 ETCONTACT=5 ! Contact Element ETTARGE=6 ! Traget Element ! CONTA171/TARGET169 als Kontaktpaarung ! ____________________________________________________________________________________ ! ! Modellerstellung ! ____________________________________________________________________________________ /PREP7 /TITLE,Temperaturfeldberechnung ! ____________________________________________________________________________________ ! ! Geometrie ! ____________________________________________________________________________________ rectng,0,b/2,0,h ! Rechteckflaeche A1 wpoffs,0,d1 ! Arbeitskoordinatensystem verschieben rectng,0,b/2-d2,0,h-d1 ! Flaeche A2 asba,1,2 ! Differenzflaeche A1-A2=A3 rectng,0,b/2-d2,0,h-d1 ! Flaeche A1 wpoffs,0,-d1 ! Arbeitskoordinatensystem verschieben rectng,0,l/2,0,-dc ! Flaeche A2 /pnum,area,1 ! Schalter fuer die graphische Flaechennummerierung (EIN) aplot ! Flaechenplot aller selektierter Flaechen Klist ! Auflistung aller selektierter Punkte (Keypoints) /pnum,kp,1 ! Schalter fuer die graphische Punktnummerierung (EIN) kplot ! Punkteplot aller selektierter Punkte (Keypoints) k,15,0,d1/4 ! Definition eines Punktes ! 15 = Punktnummer ! 0 = x-Koordinate ! d1/4 = y-Koordinate k,16,0,3*d1/4 k,17,b/2-d2,d1/4 k,18,b/2-d2,3*d1/4 /pnum,line,1 ! Schalter fuer die graphische Liniennummerierung (EIN) lplot ! Linienplot aller selektierter Linien l,15,16 ! Liniendefinition (Liniennummer = 16) ! 15 = Anfangspunkt ! 16 = Endpunkt l,16,18 l,17,18 l,15,17 AL,15,16,17,18 ! Flaechendefinition ueber Linien ! 15 = 1.Linie ! 16 = 2.Linie ! 17 = 3.Linie ! 18 = 4.Linie /pnum,kp,0 /pnum,Line,0 /pnum,area,1 ! Schalter fuer die graphische Flaechennummerierung (EIN) Aplot ! Flaechenplot aller selektierter Flaechen ! ____________________________________________________________________________________ ! ! Überlappung & Verschmeltzung ! ____________________________________________________________________________________ aovlap,3,4 ! Flaeche A3 und A4 Überlappen aglue,5,2 ! Flaeche A5 und A2 verschmelzen Aplot /triad,off ! koordinatensystem ausblenden WPSTYL,,,,,,,,0 ! Display Working Plane deaktivieren /replot SAVE ! Sicherung der momentanen Datenbasis auf file.db ! ____________________________________________________________________________________ ! ! ELEMENTVEREINBARUNG ! ____________________________________________________________________________________ ET,ETWASSER,55,,,1 ! Elementvereinbarung ! ETWASSER = Elementtypnummer ! PLANE55 = Elementtyp ! 2-D Thermal Solid (Thermal Mass Solid, Quad 4Node 55) ! 1 = Axisymmetric ! keyopt(3) = 3, plane with thickness ET,ETTOPF,55,,,1 ET,ETSCHICHT,55,,,1 ET,ETKERAMIK,55,,,1 ET,ETCONTACT,171,2 ! Elementvereinbarung ! ETCONTACT = Elementtypnummer ! CONTA171 = Elementtyp ! 2-d 2-node surface-to-surface contact ! KEYOPT(1) = 2 ;Temp keyopt,ETCONTACT,12,5 ET,ETTARGE,TARGE169 ! Elementvereinbarung ! ETTARGE = Elementtypnummer ! TARGE169 = Elementtyp ! 2-d target segment R,14 !Real Konstante ! ____________________________________________________________________________________ ! ! MATERIALVEREINBARUNG ! ____________________________________________________________________________________ MPTEMP,1,0,5,10,15,20,25 MPTEMP,7,30,35,40,45,50,60 MPTEMP,13,70,80,90,99.63,100 MPDATA,ENTH,ETWASSER,1,0.06e3,21.1e3,42.1e3,63.0e3,83.9e3,104.8e3 MPDATA,ENTH,ETWASSER,7,125.8e3,146.7e3,167.6e3,188.5e3,209.4e3,251.2e3 MPDATA,ENTH,ETWASSER,13,293.1e3,335e3,377e3,417.5e3,2675e3 MPPlot,ENTH,ETWASSER ! Specific enthalpy (J/kg) MPDATA,KXX,ETWASSER,1,561e-3,570.6e-3,580e-3,589.4e-3,598.4e-3,607.2e-3 MPDATA,KXX,ETWASSER,7,615.5e-3,623.3e-3,630.6e-3,637.3e-3,643.6e-3,654.4e-3 MPDATA,KXX,ETWASSER,13,663.1e-3,670e-3,675.3e-3,679e-3,25.08e-3 MPPlot,KXX,ETWASSER MPDATA,VISC,ETWASSER,1,1792e-6,1518e-6,1306e-6,1137e-6,1001e-6,890.4e-6 MPDATA,VISC,ETWASSER,7,797.7e-6,719.6e-6,653.3e-6,596.3e-6,547.1e-6,466.6e-6 MPDATA,VISC,ETWASSER,13,404.1e-6,354.5e-6,314.6e-6,283e-6,12.27e-6 MPPlot,VISC,ETWASSER MPTEMP MPTEMP,1,20,30,40,50,60,70 ! 6 temperatures (temps 1-6) MPTEMP,7,80,90,100 ! 6 more temps (temps 7-12) MPDATA,HF,ETWASSER,1,580,780,980,1180,1380,1580 MPDATA,HF,ETWASSER,7,1780,1980,2180 MPPlot,HF,ETWASSER ! Film Koeffizient (W/m^2*k), Konvektion mit Alpha(T) !------------------------------------------------------------------------ MP,KXX,ETTOPF,25 ! Materialvereinbarung ! KXX = Thermal conductivity ! ETWASSER = Materialgruppennummer ! 25 = Wert MP,C,ETTOPF,460 ! Materialvereinbarung ! C = Specific Heat ! ETTOPF = Materialgruppennummer ! 460 = Wert MP,DENS,ETTOPF,7700 ! Materialvereinbarung ! DENS = Density ! ETTOPF = Materialgruppennummer ! 7700 = Wert !------------------------------------------------------------------------ MP,KXX,ETSCHICHT,401 MP,C,ETSCHICHT,385 MP,DENS,ETSCHICHT,8920 !------------------------------------------------------------------------ MP,KXX,ETKERAMIK,40 MP,C,ETKERAMIK,0.5e3 MP,DENS,ETKERAMIK,3200 ! ____________________________________________________________________________________ ! ! Vernetzung der Elemente ! ____________________________________________________________________________________ ESIZE,1e-3 ! Globale Elementgrösse MSHAPE,0,2d ! Festlegung der Vernetzungsart ! 0 = Viereck-Elemente ! 2d = Area Mesh MSHKEY,2 ! Vernetzungstype, 2 = Mapped Mesh, Falls nicht Free Mesh Type,ETWASSER ! Festlegung der Materialtyp MAT,ETWASSER ! Materialset 1 aktivieren AMESH,1 ! Es wird die Flaeche 1 vernetzt Type,ETTOPF ! Festlegung der Materialtyp MAT,ETTOPF ! Materialset 2 aktivieren AMESH,6 ! Flaeche mit Material 6 vernetzen Type,ETSCHICHT ! Festlegung der Materialtyp MAT,ETSCHICHT ! Materialset 3 aktivieren AMESH,3 ! Flaeche mit Material 3 vernetzen Type,ETKERAMIK ! Festlegung der Materialtyp MAT,ETKERAMIK ! Materialset aktivieren AMESH,4 ! Flaeche mit Material 4 vernetzen /PNUM,MAT,1 ! Schalter fuer die graphische Materialnummerierung (EIN) eplot ! Elementplot aller selektierter Elemente /number,1 ! Nur Color auswaehlen /replot ! Noch mal plot ! ____________________________________________________________________________________ ! ! Selektieren ! ____________________________________________________________________________________ LSEL,S,LINE,,3 ! Linienselektierung ! S = Bildung einer neuen Linienmenge ! aus der Menge aller Linien. ! LINE = Selektierungsattribut ! 3 = Liniennummer NSLL,S,1 ! Bildung einer neuen Teilmenge Knoten aus der Menge aller Knoten ! die auf den selektierten Linien liegen. ! 1 = alle Knoten auch die Randknoten der Linien werden selektiert. Nplot cm,Xwater,node ALLS Nplot LSEL,S,LINE,,4 NSLL,S,1 ! Bildung einer neuen Teilmenge Knoten aus der Menge aller Knoten ! die auf den selektierten Linien liegen. ! 1 = alle Knoten auch die Randknoten der Linien werden selektiert. cm,Ywater,node ALLS Nplot cmsel,s,Xwater cmsel,a,Ywater cm,Water,node cmsel,s,water ! Beruehrungsbereich des Wassers mit dem Topf nplot !----------------------------------------------------------------------------------- LSEL,S,LINE,,5 NSLL,S,1 Nplot cm,Xtopf,node ALLS Nplot LSEL,S,LINE,,6 NSLL,S,1 cm,Ytopf,node ALLS Nplot cmsel,s,Xtopf cmsel,a,Ytopf cm,topf,node cmsel,s,topf ! Innere Bereich der Topf nplot !------------------------------------------------------------------------------------- nsel,s,loc,y,0,h nsel,r,loc,x,b/2 cm,auy,node nsel,s,loc,x,b/2,l/2 nsel,r,loc,y,0 cm,aux,node cmsel,s,auy cmsel,a,aux nplot cm,aussen,node ! Aeussere Bereich der Topf !------------------------------------------------------------------------------------- nsel,s,loc,y,-dc nsel,r,loc,x,0,b/2 cm,unten,node ! Untere Bereich der Glaskeramik-Platte nsel,all nplot ! ____________________________________________________________________________________ ! ! Kontaktelemente ! ____________________________________________________________________________________ type,ETCONTACT real,14 cmsel,s,topf esurf type,ETTARGE cmsel,s,water esurf ALLS nplot ! ____________________________________________________________________________________ ! ! Thermische Randbedingungen (Konvektionen) ! ____________________________________________________________________________________ cmsel,s,aussen sf,all,conv,10,20 allsel ! alles auswaehlen nplot finish ! Schiesst Prepocessor ! ____________________________________________________________________________________ ! ! Aufruf des Loesungsmoduls ! ____________________________________________________________________________________ /solu Antype,Trans ! Transiente Analyse THOPT,QUASI EQSLV,ICCG !JCG toffst,293 ! sets offset from absolute 20°c tunif=20 ! Ausgangtemperatur 20°c !Festlegen der zeitlichen Randbedingungen !---------------------------------------- AUTOTS,ON ! Zeitschrittautomatik einschalten outres,all,all ! Alle Ergebnisse abspeichern kbc,0 ! Sprungfoermiger Lastverlauf ! 0 = Schluesselwert für die Randbedingung (Ramp Chang) !Festlegen der Integrationsparameter !----------------------------------- TINTP,,,,1.0 ! Parameter für transiente Temperaturfeldberechnungen ist Theta ! Standardwert von 0.5 auf 1.0 erhoeht !Lastschritt 1 : Stationaerer Anfangszustand !------------------------------------------- TIMINT,OFF ! transiente Effekte werden voruebergehend ausgeschaltet TIME,1.0e-6 ! keine Lastzustaende bei t = 0 erlaubt !Lastschritt 2 : Transiente Rechnung (0 < t < 1000 s) !---------------------------------------------------- TIMINT,ON !Randbedingungen Waermestromdichte !--------------------------------- *dim,Tbulk1,Table,11,1 Tbulk1(1,1)=20,20,30,40,50,60,70,80,90,100,110 Tbulk1(1,0)=0,100,200,300,400,500,600,700,800,900,1000 Tbulk1(0,1)=1 TM_START_1=1E-6 ! Starting time (must be > 0) TM_END_1=zeitende ! Ending time of the transient TM_INCR=10 ! Time increment *DO,TM1,TM_START_1,TM_END_1,TM_INCR ! Do for TM from TM_START to TM_END in ! steps of TM_INCR TIME,TM1 ! Time value cmsel,s,water sf,all,conv,-ETWASSER,%Tbulk1% NSEL,ALL ! Activiere alle nodes cmsel,s,unten sf,all,hflux,15000 ! Waermestromdichte nsel,all solve *ENDDO !Lastschritt 3 : Abkühlung (1000 < t < 2000 s) !--------------------------------------------- *dim,Tbulk2,Table,10,1 Tbulk2(1,1)=100,90,80,70,60,50,40,30,20,20 Tbulk2(1,0)=1000,1200,1300,1400,1500,1600,1700,1800,1900,2000 Tbulk2(0,1)=1 TM_START_2=zeitende TM_END_2=zeitende*2 TM_INCR=10 *DO,TM,TM_START_2,TM_END_2,TM_INCR TIME,TM cmsel,s,water sf,all,conv,-ETWASSER,%Tbulk2% NSEL,ALL sfdele,all,hflux solve *ENDDO finish ! schliessen von SOLVE-TEIL