Bonjour,
J'ai réalisé une plaque carré, avec une convection à 15 degré sur un coté et une à 250 sur le coté d'en face. Mon problème est que dans la procédure pasapas, je me retrouve avec des températures supérieures à 250 degré (dans les 4000 degrés) par endroit (température initiale : 15 degré partout).
Je ne comprends pas d'ou cette erreur peut venir, Avez vous une idée ?
Merci d'avance
Sébastien
TITR ‘Plaque’ ;
OPTION DIME 2 ELEM TRI3 MODE PLAN ;
************************************************************************
************************* GEOMETRIE ************************************
************************************************************************
* POINTS
E= 15. -15. ;
F= 15. 15. ;
G= -15. 15. ;
H= -15. -15. ;
* LIGNES
DENS1= 0.5 ;
DENS2= 2.0 ;
LEF= DROI E F DINI DENS2 DFIN DENS2 ;
LFG= DROI F G DINI DENS2 DFIN DENS2 ;
LGH= DROI G H DINI DENS2 DFIN DENS2 ;
LHE= DROI H E DINI DENS2 DFIN DENS2 ;
* SURFACE
conte= LEF ET LFG ET LGH ET LHE ;
omega1= conte;
SURF1=SURF omega1 ;
* VISUALISATION
TITR ‘THERM1 MAILLAGE’ ;
TRAC QUAL SURF1 ;
************************************************************************
*************** CHOIX DU COMPORTEMENT ET DU MODELE EF ******************
************************************************************************
MOD1 = MODE SURF1 THERMIQUE ISOTROPE ;
MOD2= MODE LFG CONVECTION ;
MOD3= MODE LHE CONVECTION ;
* DONNEES MATERIAU ET ELEMENTS
MAT1=MATE MOD1 'K' 1.5 'C' 500. 'RHO' 1000.;
MAT2=MATE MOD2 'H' 10.0 ;
MAT3=MATE MOD3 'H' 10.0 ;
* CONDITIONS AUX LIMITES
* FLUX
FLU1=CONV MOD2 MAT2 T 15. ;
FLU2=CONV MOD3 MAT3 T 250. ;
**************************************************************************
************************** RESOLUTION ************************************
**************************************************************************
* MATRICES DE CONDUCTIVITE
COND1=COND MOD1 MAT1 ;
COND2=COND MOD2 MAT2 ;
COND3=COND MOD3 MAT3 ;
CONDTOT=COND1 ET COND2 ET COND3 ;
* MATRICE DE CAPACITE
CAP1 = CAPACITE MOD1 MAT1 ;
* FLUX
FLUTOT=FLU1 ET FLU2 ;
* RESOLUTION
CHPO1=RESO CONDTOT FLUTOT ;
***************************************************************************
*************************** PAS-A-PAS *************************************
***************************************************************************
MOD4 = MOD1 ET MOD2 ET MOD3 ;
MAT4 = MAT1 ET MAT2 ET MAT3 ;
TEMPS;
ETAB = TABLE;
CHTER0 = MANU CHPO SURF1 1 T 15. ;
ETAB.'TEMPERATURES' = TABLE;
ETAB.'TEMPERATURES' . 0 = CHTER0 ;
ETAB.’CONDUCTIVITE’ = CONDTOT ;
ETAB.’CAPACITE’ = CAP1 ;
ETAB.'MODELE' = MOD4;
ETAB.'CARACTERISTIQUES' = MAT4;
PP1 = PROG 0. PAS 100000. NPAS 50 ;
PP2 = PROG 51. * 1. ;
MEVOL1 = EVOL MANU ‘TEMPS’ PP1 ‘COEFF’ PP2 ;
DESS MEVOL1;
CHA1 = CHAR 'TECO' FLU1 MEVOL1;
CHA2 = CHAR 'TECO' FLU2 MEVOL1;
CHA3 = CHA1 ET CHA2;
ETAB.'CHARGEMENT' = CHA3;
ETAB.'TEMPS_CALCULES' = PROG 100000. 500000. 1000000. ;
ETAB.'PROCEDURE_THERMIQUE' = LINEAIRE;
PASAPAS ETAB ;
***************************************************************************
******************** POST-TRAITEMENT PASAPAS ******************************
***************************************************************************
*début de la boucle d'animation
REPETER bloc0 3 ;
i = &bloc0 ;
te1 = (ETAB.temperatures.i) ;
TRAC te1 'NCLK' SURF1 ;
hjgj
FIN bloc0 ; |