* *nom fichier * TnomFich = 'Pb2_Convection+Rayonnement.res' ; * *Calcul thermique transitoire 2D ; *Proprietes thermique constantes ou variable * * OPTI DIME 2 ELEM QUA4 MODE 'PLAN' ; GRAPH = FAUX ; * *-----------------------Maillage--------------------------------------- * P1 = 0. 0. ; P2 = 0.04 0. ; P3 = 0.04 0.4; P4 = 0.0 0.4; * NELM1 = 40 ; NELM2 = 120 ; * L12 = DROI NELM1 P1 P2; L23 = DROI NELM2 P2 P3; L34 = DROI NELM1 P3 P4; L41 = DROI NELM2 P4 P1; LTOT = L12 ET L23 ET L34 ET L41; * GEO1 = SURF LTOT PLAN ; * *---------------Preparation donnees Materiaux------------------------- * ***Conductivite thermique Lam = EVOL 'MANU' 'T' (PROG 0. 100. 105. 200. 400. 600. 800. 1000. 1200.) 'K' (PROG 0.2 0.2 0.11 0.9 0.11 0.13 0.15 0.3 0.5) ; ***Capacite calorifique Cp = EVOL 'MANU' 'T' (PROG 0. 100. 150. 200. 1000. ) 'C' (PROG 1000. 1100. 2500. 1500. 1700. ) ; ***Masse volumique RHO = EVOL 'MANU' 'T' (PROG 0. 100. 105. 150. 200. 1000. ) 'RHO' (PROG 1500. 1500. 1400. 1390. 1400. 1430. ) ; ***Temperature de changement de phase T_fus = 150. ; ***Charleur latent Qlat = 19.E4 ; ***Emissivite Emis = 0.9 ; ***Coefficient d'echange thermique H = 25. ; ***Temperature ambiante T_ext = 18. ; ***Temperature initiale T_ini = 18. ; ***Temps de calcul final Temps_F = 180 * 60 ; *---------------Preparation Cycles temperature externe----------------- PrTem0 = PROG 0. PAS 1. (temps_f / 60.) ; PrText = (PROG ((Temps_F / 60) + 1)*T_ext) + (345. * (PROG 'LOGA' 'A' 8. 'B' 1. PrTem0)); * ***Temps en seconde et temperature en celsius * PrTem0 = PrTem0 * 60 ; PrText = PrText ; * *------------------------Modele et Materiaux--------------------------- * MOD_THER = MODE GEO1 'THERMIQUE' 'PHASE'; MAT_THER = MATE MOD_THER 'K' LAM 'C' Cp 'RHO' RHO 'TPHA' T_FUS 'QLAT' QLAT; * MOD_CON = MODE (L41 ET L23) 'CONVECTION' ; MAT_CON = MATE MOD_CON 'H' H ; * MOD_RAY = MODE (L41 ET L23) 'RAYONNEMENT' ; MAT_RAY = MATE MOD_RAY 'EMIS' Emis ; * *-----------------------Chargement thermique--------------------------- * EV1 = EVOL 'JAUN' 'MANU' PrTem0 PrText ; EV2 = EVOL 'MANU' (PROG 0. Temps_F) (PROG T_EXT T_EXT) ; ***Convection CHTCON1 = MANU 'CHPO' L41 1 'T' 1.; CHA1 = CHAR 'TECO' CHTCON1 EV1 ; * CHTCON2 = MANU 'CHPO' L23 1 'T' 1. ; CHA2 = CHAR 'TECO' CHTCON2 EV2 ; ***Rayonnement CHTRAY = MANU 'CHPO' L23 1 'T' 1.; CHA3 = CHAR 'TERA' CHTRAY EV2 ; *---------------------Champ de temperature initiale-------------------- * TemIni = MANU 'CHPO' GEO1 1 'T' T_Ini ; * *--------------Temps de calcul et temps de calcul final---------------- * TempCal = PROG 0. PAS 60. Temps_F ; TempSav = PROG 0. PAS 60. Temps_F ; * *------------------Creation de la table pour PASAPAS------------------- * TAB1 = TABLE ; TAB1 . 'MODELE' = MOD_THER et MOD_CON; TAB1 . 'CARACTERISTIQUES' = MAT_THER ET MAT_CON ET MAT_RAY; *TAB1 . 'BLOCAGES_THERMIQUES' = CL_TH ; TAB1 . 'TEMPERATURES' = TABLE ; TAB1 . 'TEMPERATURES' . 0 = TemIni ; TAB1 . 'CHARGEMENT' = CHA1 ET CHA2 ET CHA3; TAB1 . 'PROCEDURE_THERMIQUE' = 'NONLINEAIRE' ; TAB1 . 'RELAXATION_THETA' = 1. ; TAB1 . 'RAYONNEMENT' = TABLE ; TAB1 . 'RAYONNEMENT'. 1 = TABLE ; TAB1 . 'RAYONNEMENT'. 1 .'TYPE' = 'INFINI' ; TAB1 . 'RAYONNEMENT'. 1 .'MODELE' = MOD_RAY ; TAB1 . 'TEMPS_CALCULES' = TempCal ; TAB1 . 'TEMPS_SAUVES' = TempSav ; TAB1 . 'CELSIUS' = VRAI ; * *---------------------Appel a la procedure PASAPAS--------------------- * PASAPAS TAB1 ; * *---------------------------Post traitement---------------------------- LTimSav = TAB1.'TEMPS_SAUVES'; * I = 1; REPETER BOUCTE ('DIME' TAB1.'TEMPERATURES') ; * SI (I EGA 1) ; PrTemP1 = PROG ((EXTR (TAB1.TEMPERATURES.(I-1)) 'T' P1)) ; PrTemP2 = PROG ((EXTR (TAB1.TEMPERATURES.(I-1)) 'T' P2)) ; SINON; PrTemP1 = PrTemP1 ET (PROG ((EXTR (TAB1.TEMPERATURES.(I-1)) 'T' P1))) ; PrTemP2 = PrTemP2 ET (PROG ((EXTR (TAB1.TEMPERATURES.(I-1)) 'T' P2))) ; FINSI; * I = I + 1; * FIN BOUCTE ; * EvTemP1 = EVOL 'ROUG' 'MANU' 'Temps' LTimSav 'Temp_P1' PrTemP1 ; EvTemP2 = EVOL 'BLAN' 'MANU' 'Temps' LTimSav 'Temp_P2' PrTemP2; dess (ev1 et evtemP1 et evtemp2) ; OPTI SAUV FORMAT TNomFich; SAUV TAB1; * * FIN;