Code_Aster ®
Version
7.4
Titre :
Opérateur THER_NON_LINE_MO
Date :
20/01/05
Auteur(s) :
C. DURAND Clé
:
U4.54.03-D Page
: 1/6
Organisme(s) : EDF-R&D/AMA
Manuel d'Utilisation
Fascicule U4.5- : Méthodes de résolution
Document : U4.54.03
Opérateur THER_NON_LINE_MO
1 But
Calculer la réponse thermique stationnaire non linéaire avec chargement mobile.
Cet opérateur permet de résoudre l'équation de la chaleur stationnaire dans un référentiel mobile lié à
un chargement et se déplaçant dans une direction et à une vitesse données [R5.02.04].
Les non linéarités du problème proviennent aussi bien des caractéristiques du matériau qui dépendent
de la température, que des conditions aux limites de type rayonnement.
Les problèmes de ce type peuvent être traités avec des modèles utilisant des éléments finis de
structure plans, axisymétriques et tridimensionnels.
Le concept produit par cet opérateur est de type evol_ther comme pour une analyse par
THER_LINEAIRE [U4.54.01] ou THER_NON_LINE [U4.54.02], mais il ne contient qu'un seul numéro
d'ordre puisqu'on réalise un calcul stationnaire.
Manuel d'Utilisation
Fascicule U4.5- : Méthodes de résolution
HT-66/05/004/A
Code_Aster ®
Version
7.4
Titre :
Opérateur THER_NON_LINE_MO
Date :
20/01/05
Auteur(s) :
C. DURAND Clé
:
U4.54.03-D Page
: 2/6
2 Syntaxe
temper [evol_ther] = THER_NON_LINE_MO
(
MODELE
=
mo,
[modele]
CHAM_MATER =
chmat,
[cham_mater]
EXCIT =_F (
CHARGE
=
char,
[char_ther]
),
TEMP_INIT =_F (
EVOL_THER
=
temp,
[evol_ther]
NUME_INIT
= / nuini_evol, [I]
/
0,
[DEFAUT]
),
CONVERGENCE =_F(
CRIT_TEMP_RELA = / 1.E-3,
[DEFAUT]
/
testr,
[R]
CRIT_ENTH_RELA = / 1.E-2,
[DEFAUT]
/
testl,
[R]
ITER_GLOB_MAXI = / 10,
[DEFAUT]
/
iterl,
[R]
),
SOLVEUR =_F( voir document [U4.50.01])
),
TITRE
=
titre,
[l_Kn]
);
Manuel d'Utilisation
Fascicule U4.5- : Méthodes de résolution
HT-66/05/004/A
Code_Aster ®
Version
7.4
Titre :
Opérateur THER_NON_LINE_MO
Date :
20/01/05
Auteur(s) :
C. DURAND Clé
:
U4.54.03-D Page
: 3/6
3 Opérandes
3.1 Opérande
MODELE
MODELE = mo
Nom du modèle dont les éléments font l'objet du calcul thermique.
3.2 Opérande
CHAM_MATER
CHAM_MATER = chmat
Nom du champ de matériau affecté sur le modèle.
3.3 Mot
clé
EXCIT
EXCIT
=
Mot clé facteur permettant de définir plusieurs chargements. Pour chaque occurrence du mot clé
facteur, on définit une charge constante en temps.
3.3.1 Opérande
CHARGE
CHARGE = char
Concept de type charge produit par AFFE_CHAR_THER ou par AFFE_CHAR_THER_F
[U4.44.02].
Remarques importantes :
Pour chaque occurence du mot clé facteur EXCIT les différents concepts char utilisés
doivent être construits sur le même modèle mo.
Par ailleurs, l'opérateur AFFE_CHAR_THER ou AFFE_CHAR_THER_F doit obligatoirement
renseigner sous le mot clé facteur CONVECTION un concept de type cham_no_DEPL_R
définissant la vitesse de déplacement du domaine et produit par AFFE_CHAM_NO.
Il faut veiller à définir des conditions aux limites cohérentes avec le défilement de la pièce
(ou du chargement, selon le référentiel choisi) et l'hypothèse de stationnarité :
Si le solide est animé d'un mouvement de translation uniforme à vitesse constante, il est
équivalent de considérer le champ de température calculé :
· comme une photographie à un instant donné,
· ou comme une succession temporelle de la carte de température dans une section
d'abscisse fixe, orthogonale à la direction de défilement.
Chargement mobile
v = vitesse de défilement
C.L. température
C.L. adiabatique
imposée
temps
t = 0
t = infini
Ainsi, en aval dans l'espace (tôt dans le temps), la condition aux limites sur le bord de la
pièce doit imposer la température de la pièce avant que n'arrive le chargement mobile.
En amont dans l'espace (ou tard dans le temps), un flux thermique nul (adiabaticité)
garantit l'hypothèse de stationnarité.
Manuel d'Utilisation
Fascicule U4.5- : Méthodes de résolution
HT-66/05/004/A
Code_Aster ®
Version
7.4
Titre :
Opérateur THER_NON_LINE_MO
Date :
20/01/05
Auteur(s) :
C. DURAND Clé
:
U4.54.03-D Page
: 4/6
3.4 Mot
clé
TEMP_INIT
TEMP_INIT = litps
Permet de définir le champ initial à partir duquel le calcul en repère mobile sous condition de
stationnarité est effectué.
3.4.1 Opérande
EVOL_THER
EVOL_THER = temp
Nom de la structure de donnée de type evol_ther dont est extrait le champ de
température à partir duquel est calculée la solution au problème.
3.4.2 Opérande
NUME_INIT
NUME_INIT = / nuini_evol
/
0
[DEFAUT]
Numéro d'ordre du champ à extraire de cette structure de donnée.
Remarque :
Pour ce type de calcul, il n'y a pas de condition initiale à définir (elle se définit avec une
condition aux limites sur un bord du domaine défini par le modèle). Mais, on peut
néanmoins réaliser le calcul stationnaire à partir d'un champ de température déjà calculé
comme solution du problème pour une partie de la charge seulement. Cela peut être utile
dans le cas où l'application de la totalité du chargement rend la convergence difficile. On
peut alors enchaîner plusieurs résolutions avec un chargement qui évolue
progressivement et prendre comme premier itéré d'un calcul courant la solution du calcul
précédent. Concrètement, on utilise alors le mot clé TEMP_INIT pour définir le champ de
température avec lequel on veut initialiser le calcul courant.
3.5 Mot
clé
CONVERGENCE
CONVERGENCE =
Permet de définir les valeurs associées aux critères de convergence :
3.5.1 Opérande
CRIT_TEMP_RELA
CRIT_TEMP_RELA = / 1.e-3
/
testr
L'algorithme continue les itérations tant que la variation relative de l'itéré est supérieure à
testr.
où
désigne la valeur de la température au noeud i pour la nième itération.
Manuel d'Utilisation
Fascicule U4.5- : Méthodes de résolution
HT-66/05/004/A
Code_Aster ®
Version
7.4
Titre :
Opérateur THER_NON_LINE_MO
Date :
20/01/05
Auteur(s) :
C. DURAND Clé
:
U4.54.03-D Page
: 5/6
3.5.2 Opérande
CRIT_ENTH_RELA
CRIT_ENTH_RELA = / 1.e-2
/
testl
L'algorithme continue les itérations tant que la variation relative de l'itéré est supérieure à
testl.
où
désigne la valeur de l'enthalpie à la nième itération au iième point de Gauss.
3.5.3 Opérande
ITER_GLOB_MAXI
ITER_GLOB_MAXI = / 10
/
iterl
L'algorithme continue les itérations externes tant que le nombre d'itérations externes est
inférieur à iterl.
3.6 Mot
clé
SOLVEUR
SOLVEUR =
Ce mot clé facteur est facultatif : il permet de choisir un autre solveur de résolution de système.
Cet opérande est commun à l'ensemble des commandes globales [U4.50.01].
3.7 Opérande
TITRE
TITRE = titre
Titre que l'on veut donner au résultat stocké dans temp structure de données de type evol_ther
[U4.03.01].
4 Modélisation
Les problèmes de thermique non linéaire en repère mobile et sous condition de stationnarité peuvent
être traités avec des modèles utilisant les éléments finis 3D ou 2D plan décrits dans les documents
[U3.23.01] et [U3.24.01].
Manuel d'Utilisation
Fascicule U4.5- : Méthodes de résolution
HT-66/05/004/A
Code_Aster ®
Version
7.4
Titre :
Opérateur THER_NON_LINE_MO
Date :
20/01/05
Auteur(s) :
C. DURAND Clé
:
U4.54.03-D Page
: 6/6
5 Exemples
On présente ci-dessous les principales commandes utilisées pour effectuer un calcul de thermique
non-linéaire en repère mobile et sous condition de stationnarité.
conduc =
DEFI_FONCTION (
NOM_PARA
=
'TEMP',
VALE =
(
20.0,
95.0,
500.0, 325.0,
1000.0, 700.0
),
PROL_DROIT =
'CONSTANT',
PROL_GAUCHE=
'CONSTANT',
)
enthal =
DEFI_FONCTION
(
NOM_PARA
=
'TEMP',
VALE =
(
0.0, 0.0,
500.0, 1.40E9,
510.0, 3.50E9,
1200.0, 4.70E9
),
PROL_DROIT =
'LINEAIRE',
PROL_GAUCHE=
'LINEAIRE',
)
mat =
DEFI_MATERIAU
(THER_NL=
_F( LAMBDA =
conduc,
BETA =
enthal
),)
...
chavi = AFFE_CHAM_NO
( MAILLAGE = mail,
GRANDEUR = 'DEPL_R',
AFFE =_F( TOUT = 'OUI',
NOM_CMP = ('DX','DY'),
VALE_R = (1.0 , 0.0) )
char = AFFE_CHAR_THER
(
MODELE =
model,
FLUX_REP =_F( GROUP_MA = sou,
FLUN = 5.E6),
CONVECTION =_F( VITESSE =
chavi)
)
tempe
=
THER_NON_LINE_MO
(
MODELE
=
model,
CHAM_MATER =
chmat,
EXCIT
=_F(CHARGE =
char
),
)
...
FIN()
Manuel d'Utilisation
Fascicule U4.5- : Méthodes de résolution
HT-66/05/004/A
Document Outline