modélisation thermique

Bonjour.
En cette dernière année ou des sites liés je trouvai un script pour Salome-Meca avec le Code-Aster ayant licence GPL pour la solution d'un problème thermique non linéaire, à savoir le calcul des températures dans le béton soumis à un incendie courbe ISO783 normalisée l'auteur est M. F. Renou.
Le problème a été très bien fait, mais en 2D proposition de modélisation n'a pas la mise en œuvre de deux conditions aux limites, que je décris ci-dessous.
Le nom de fichier du script original est "thernonlin.comm", ici proposer à nouveau mais sous une forme réduite à concentrer le problème.

DEBUT();

MESH=LIRE_MAILLAGE(UNITE=20,
FORMAT='MED',);

#14400 secondi [240 minuti] con 48 step di [5.0 minuti]
LISTC=DEFI_LIST_REEL(DEBUT=0.0,
INTERVALLE=_F(JUSQU_A=14400.0,
NOMBRE=48,),);

#1400 passi di integrazione
LISTT=DEFI_LIST_REEL(DEBUT=0.0,
INTERVALLE=_F(JUSQU_A=1400.0,
NOMBRE=1400,),);

#costante di Stefan-Boltzmann [W/m2K4]
sig=DEFI_CONSTANTE(VALE=5.6703E-8,);

#fattore di emissività risultante [EC2]
eps=DEFI_CONSTANTE(VALE=0.8,);

#convezione pareti esposte al fuoco [Wm2C]
conv=DEFI_CONSTANTE(VALE=25.0,);

#convezione pareti esposte all'aria a temperatura ambiente [Wm2C]
convest=DEFI_CONSTANTE(VALE=9.0,);

#convezione cavità interne [Wm2C]
convint=DEFI_CONSTANTE(VALE=19.0,);

#Nominal standard curve
tempga = FORMULE(VALE='20+345*log10(8/60.0*INST+1)',
NOM_PARA='INST',);
TEMPGAZ=CALC_FONC_INTERP(FONCTION=tempga,
LIST_PARA=LISTC,);

tabtempg=CREA_TABLE(FONCTION=_F(FONCTION=TEMPGAZ,),);

IMPR_TABLE(TABLE=tabtempg,);

#conduttività termica ksup [W/mk]
lcsup = FORMULE(VALE='2.0-0.245*(TEMP/100.)+0.0107*(TEMP/100.)**2',
NOM_PARA='TEMP',);
LCSUPT=CALC_FONC_INTERP(FONCTION=lcsup,
LIST_PARA=LISTT,);
tablci=CREA_TABLE(FONCTION=_F(FONCTION=LCSUPT,),);
IMPR_TABLE(TABLE=tablci,);

#capacità termica in funzione del tenore di acqua
#cv15 per 15% default [EC2]
cv15=DEFI_FONCTION(NOM_PARA='TEMP',VALE=(0.0,2070000.0,
99.0,2070000.0,
100.0,3381000.,
115.,3381000.,
200.0,2254000.0,
400.0,2403500.0,
1200.0,2270675.0,
),);

tabcv15=CREA_TABLE(FONCTION=_F(FONCTION=cv15,),);

IMPR_TABLE(TABLE=tabcv15,);

#conduttività superiore [W/mk]
MATER=DEFI_MATERIAU(THER_NL=_F(LAMBDA=LCSUPT,
RHO_CP=cv15,),);

MODEL=AFFE_MODELE(MAILLAGE=MESH,
AFFE=_F(TOUT='OUI',
PHENOMENE='THERMIQUE',
MODELISATION='PLAN',),);

MATFIELD=AFFE_MATERIAU(MAILLAGE=MESH,
AFFE=_F(TOUT='OUI',
MATER=MATER,),);

loadF=AFFE_CHAR_THER_F(MODELE=MODEL,
RAYONNEMENT=_F(GROUP_MA='Grfeu',
SIGMA=sig,
EPSILON=eps,
TEMP_EXT=tempga,),
ECHANGE=_F(GROUP_MA='Grfeu',
COEF_H=conv,
TEMP_EXT=tempga,),);

RESULT=THER_NON_LINE(MODELE=MODEL,
CHAM_MATER=MATFIELD,
EXCIT=_F(CHARGE=loadF,),
INCREMENT=_F(LIST_INST=LISTC,),
ETAT_INIT=_F(VALE=20.0,),
NEWTON=_F(ITER_LINE_MAXI=15,),);

#risultati in formato binario [*.med]
IMPR_RESU(FORMAT='MED',
RESU=_F(RESULTAT=RESULT,),);

FIN();

Dans la pratique dans la modélisation que nous considérons que la convection dans les murs directement exposées aux gaz chauds de l'incendie qui sont inclus dans le groupe "Grfeu" que la ligne suivante:
conv=DEFI_CONSTANTE(VALE=25.0,);

Je voudrais ajouter la convection des murs exposés à l'environnement externe, appelé le groupe "Grest"
convest=DEFI_CONSTANTE(VALE=9.0,);

et par convection des parois de toutes les cavités internes, définies dans le groupe "Grint":
convint=DEFI_CONSTANTE(VALE=19.0,);

Comment puis-je ajouter dans le calcul de ces deux conditions?

Merci

Joindre des fichiers à un type de modélisation.
TRL_Forata