RDforum.free.fr/Cast3M   Index du Forum
Site Web Cast3M
RDforum.free.fr/Cast3M
Forum d'entraide des utilisateurs de Cast3M
Répondre au sujet
Problème de convection
sebastien


Inscrit le: 27 Mai 2010
Messages: 2
Répondre en citant
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

Code:
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 ;
Voir le profil de l'utilisateurEnvoyer un message privé
Duong-ParisEst


Inscrit le: 04 Mar 2008
Messages: 9
Localisation: Paris-Est
Répondre en citant
Bonjour,
En fait la solution "CHPO1=RESO CONDTOT FLUTOT;" est déjà la réponse du problème stationnaire pourquoi il faut encore passer PASAPAS?
Voir le profil de l'utilisateurEnvoyer un message privéEnvoyer un e-mailYahoo MessengerMSN Messenger
Problème de convection
Vous ne pouvez pas poster de nouveaux sujets dans ce forum
Vous ne pouvez pas répondre aux sujets dans ce forum
Vous ne pouvez pas éditer vos messages dans ce forum
Vous ne pouvez pas supprimer vos messages dans ce forum
Vous ne pouvez pas voter dans les sondages de ce forum
Vous pouvez joindre des fichiers
Vous pouvez télécharger des fichiers
Toutes les heures sont au format GMT + 2 Heures  
Page 1 sur 1  

  
  
 Répondre au sujet