Code_Aster ®
Version
5.0
Titre :
Thermique non linéaire
Date :
22/06/00
Auteur(s) :
C. DURAND
Clé :
R5.02.02-A
Page :
1/16
Organisme(s) : EDF/MTI/MMN
Manuel de Référence
Fascicule R5.02 : Thermique
Document : R5.02.02

Thermique non linéaire
Résumé
L'opérateur THER_NON_LINE [U4.54.02] permet de résoudre les problèmes de thermique transitoire dans les
solides en présence de non-linéarités des propriétés des matériaux (capacité calorifique et conductivité), ou des
conditions aux limites (échange thermique de type rayonnement). On présente ici la formulation et l'algorithme
employé, ce dernier étant proche de celui lié à l'opérateur STAT_NON_LINE [R5.03.01]. Les différentes options
de calcul nécessaires ont été présentées dans les éléments de structure plans, axisymétriques et
tridimensionnels [U3.22.01], [U3.23.01] et [U3.24.01].
Manuel de Référence
Fascicule R5.02 : Thermique
HI-75/99/013 - Ind A

Code_Aster ®
Version
5.0
Titre :
Thermique non linéaire
Date :
22/06/00
Auteur(s) :
C. DURAND
Clé :
R5.02.02-A
Page :
2/16
Table des matières
1 Expression de l'équation de la chaleur en thermique non linéaire.........................................................3
1.1 Equation de la chaleur pour un solide immobile..............................................................................3
1.2 Loi de Fourier ..................................................................................................................................4
1.3 Equation de la chaleur dans le cas du modèle de thermique transitoire non-linéaire .....................4
1.4 Avantage numérique de la formulation en enthalpie pour les problemes avec changement de
phase...............................................................................................................................................5
2 Conditions aux limites, chargement et condition initiale ........................................................................5
2.1 Flux normal non-linéaire..................................................................................................................5
2.2 Flux normal non-linéaire - condition de type rayonnement à l'infini.................................................6
3 Formulation variationnelle du problème.................................................................................................6
4 Discrétisation en temps de l'équation différentielle................................................................................8
4.1 Introduction de la -méthode...........................................................................................................8
4.2 Application à l'équation de la chaleur ..............................................................................................8
5 Discrétisation spatiale et adaptation de l'algorithme de Newton au problème.....................................10
5.1 Discrétisation spatiale....................................................................................................................10
5.2 Calcul stationnaire .........................................................................................................................11
5.3 Calcul transitoire............................................................................................................................12
5.4 Convergence .................................................................................................................................13
6 Principales options de thermique non-linéaire calculées dans le Code_Aster ....................................14
6.1 Conditions aux limites et chargements..........................................................................................14
6.2 Calcul des matrices élémentaires et terme transitoire ..................................................................14
6.3 Calcul du résidu.............................................................................................................................14
7 Bibliographie ........................................................................................................................................15
Manuel de Référence
Fascicule R5.02 : Thermique
HI-75/99/013 - Ind A

Code_Aster ®
Version
5.0
Titre :
Thermique non linéaire
Date :
22/06/00
Auteur(s) :
C. DURAND
Clé :
R5.02.02-A
Page :
3/16
1
Expression de l'équation de la chaleur en thermique non
linéaire

1.1
Equation de la chaleur pour un solide immobile
Dans tout ce document, on ne traite que la thermique des corps solides, même si le changement de
phase liquide/solide est pris en compte. Il n'y a donc pas d'échange de chaleur par convection mais
seulement par conduction.
Le premier principe de la thermodynamique relie la variation temporelle d'énergie totale dEtotale d'un
système compris dans un volume de contrôle au travail des efforts extérieurs W et à la chaleur
Q reçus par ce même système :
dE
= d(E
+ E
) = W
+ Q
totale
interne
cinétique

éq 1.1-1
En injectant le théorème de l'énergie cinétique dans cette équation, on fait apparaître ainsi la
puissance des efforts intérieurs, fonction du champ de vitesse [bib1] :
!Einterne = !Q - P (u)
i
éq 1.1-2
Pour la résolution du problème de thermique, le système est supposé sans mouvement. La puissance
des efforts intérieurs Pi (u) est donc nulle. En effet, dans la majorité des applications visées, les
phénomènes thermiques et mécaniques sont découplés ; la puissance volumique dissipée par les
déformations plastiques, Pi = c .!plastique , est négligée devant la chaleur échangée en surface ou
les autres sources de chaleur volumiques.
L'équation [éq 1.1-2] qui exprime la variation de la chaleur dans le volume s'écrit alors :
d
S

e d = Q! = (rvol - div q d
éq 1.1-3
dt

)
S
S
où on a noté :
e
l'énergie interne,

la masse volumique,
rvol
le taux volumique d'apport extérieur de chaleur,
q
le vecteur flux de chaleur.
En outre, puisque le solide est immobile, pour tout volume de contrôle (
t) = , on obtient alors
l'équation locale de conservation de la chaleur :
de

= r - div q
éq 1.1-4
dt
vol
Dans le cas où tout le système est animé d'un mouvement de corps rigide, un terme supplémentaire
apparaît dans le membre de gauche, faisant intervenir la vitesse du solide et le gradient de l'énergie.
Cette situation est traitée par la commande THER_NON_LINE_MO [R5.02.04].
Manuel de Référence
Fascicule R5.02 : Thermique
HI-75/99/013 - Ind A

Code_Aster ®
Version
5.0
Titre :
Thermique non linéaire
Date :
22/06/00
Auteur(s) :
C. DURAND
Clé :
R5.02.02-A
Page :
4/16
Dans le cas d'une transformation réversible, l'équation [éq 1.1-4] devient, avec l'aide du second
principe de la thermodynamique qui permet d'écrire dans notre cas dE
= TdS :
interne
T s! = r - divq
éq 1.1-5
vol
et enfin l'équation de la chaleur sous sa forme classique :
C T! = r - div q
éq 1.1-6
P
vol
s
avec la capacité calorifique à pression constante définie par : C = T
P
T
P
Comme il est expliqué au chapitre 1.4, il peut être avantageux d'écrire le terme de gauche de l'équation
[éq 1.1-6] avec l'enthalpie qui ne dépend alors que de la température :
!
= r - div q
vol
éq 1.1-7
T
où (T) =
C dT
P
T0
1.2
Loi de Fourier
En conduction thermique, la loi de FOURIER fournit une équation reliant le flux de chaleur au gradient
de la température (vecteur normal à la surface isotherme). Cette loi fait apparaître, dans sa forme la
plus générale, un tenseur de conductivité. Dans le cas d'un matériau isotrope, ce tenseur se réduit à un
coefficient (pouvant dépendre de la température), la conductivité thermique :
q(x,t) = - (T) T
(x,t)
éq 1.2-1
1.3
Equation de la chaleur dans le cas du modèle de thermique
transitoire non-linéaire

En combinant les équations [éq 1.1-5] et [éq 1.2-1], on obtient :
d
r
div ( (T) T
vol -
-
) =
éq 1.3-1
dt
ou, si la capacité calorifique ne dépend pas de la température :
dT
r
- div(-(T) T
) = C
vol
p
éq 1.3-2
dt
Manuel de Référence
Fascicule R5.02 : Thermique
HI-75/99/013 - Ind A

Code_Aster ®
Version
5.0
Titre :
Thermique non linéaire
Date :
22/06/00
Auteur(s) :
C. DURAND
Clé :
R5.02.02-A
Page :
5/16
1.4
Avantage numérique de la formulation en enthalpie pour les
problemes avec changement de phase.

La relation entre enthalpie et capacité calorifique est :

T
(T) =
C (u) du

p
T0
Quand cette fonction enthalpie présente des variations brusques, il est plus précis de manipuler (T)
que sa dérivée. En effet, les allures caractéristiques de ces fonctions au voisinage du point de fusion
sont les suivantes :
Enthalpie
Liquide
C
H(T )
C
h
2
p
a
l
eu
r
late
nte
C (T )
p
1
C (T )
p
2
Solide
Solide
H(T )
Liquide
1
Température
Température
Au cours d'une itération, soit parce que le transitoire thermique est violent, soit parce que la plage de
changement de phase est très petite (corps pur), les deux itérés successifs de la température peuvent
se situer de part et d'autre de la discontinuité. L'évaluation de la pente de la fonction enthalpie au
voisinage du point de fusion sera très fausse si on considère C (T
p
1), C (T
p
2 ) ou une moyenne
pondérée des deux. En revanche, la pente de la droite en pointillés est toujours une approximation
correcte de d dT au point de fusion.
2
Conditions aux limites, chargement et condition initiale
On se reportera à [R5.02.01] pour les conditions aux limites thermiques et les chargements conduisant
à des équations linéaires en température ainsi que pour la condition initiale.
2.1
Flux normal non-linéaire
Ce sont des conditions de type Neumann, définissant le flux entrant dans le domaine.
- q (x,t).n = g (x,T) sur la frontière
éq 2.1-1
g (x,T) est une fonction de la température et éventuellement de la variable d'espace x et/ou de
temps t et n désigne la normale extérieure à la frontière , q est le vecteur flux de chaleur (dirigé
suivant les températures décroissantes).
Cette expression permet d'introduire par exemple des conditions du type échange avec un coefficient
d'échange convectif dépendant de la température :
- q (x,t).n = g (x,T) = h(x,T) (T (x,t)
ext
- T)
éq 2.1-2
Manuel de Référence
Fascicule R5.02 : Thermique
HI-75/99/013 - Ind A

Code_Aster ®
Version
5.0
Titre :
Thermique non linéaire
Date :
22/06/00
Auteur(s) :
C. DURAND
Clé :
R5.02.02-A
Page :
6/16
2.2
Flux normal non-linéaire - condition de type rayonnement à l'infini
Un cas particulier des conditions aux limites précédentes est le rayonnement à l'infini de corps gris qui
se traduit par un cas particulier de fonction g (x,T) :
- q (x,t).n = [(T(x) + .
273 )
15 4 - (
4

T +
.
273 )
15 ]
éq 2.2-1
Les caractéristiques à définir lors de la définition de ce chargement sont l'émissivité , la constante de
Stefan-Boltzmann =5,73.10­8 usi et la température à l'infini.
T(r) et T sont alors exprimées en degrés Celsius. ­273.15°C est la température du zéro absolu.
3
Formulation variationnelle du problème
Nous nous bornerons ici à présenter le problème avec uniquement les conditions aux limites de
température imposée [R5.02.01 §2.1], de flux normal imposé [R5.02.01 §2.3], d'échange
[R5.02.01 §2.4], de flux non linéaire [§2.1] et de rayonnement [§2.2].
Soit un ouvert de R3, de frontière = 1 2 3 4 5.
On doit résoudre l'équation [éq 1.1-4] en T sur × ]0,t[ avec les conditions aux limites :
T
= Td (r,t)
sur

1
T
(T)
= f (r,t)
sur

n
2
T
(T)
= h(r,t) (T (r,t) - T

)
sur
n
ext
3
éq 3-1

T
(T)
= g(r,T)
sur
n
4

T
(T)
=
4
4

[(T + .
273 )
15 - (T +
.
273 )
15 ] sur

n
5
et avec, éventuellement, des conditions initiales T(t = )
0 . Si ces dernières ne sont pas précisées, on
résout au préalable le problème stationnaire, c'est à dire l'équation [éq 1.3-1] sans le terme d'évolution
temporelle.
Soit v une fonction suffisamment régulière s'annulant sur 1, en remarquant :
d

(T).v d
!(T).v d

=
dt


T

éq 3-2
(T) T
.v d = - div((T) T
).v d + (T)
.v d





n



Manuel de Référence
Fascicule R5.02 : Thermique
HI-75/99/013 - Ind A

Code_Aster ®
Version
5.0
Titre :
Thermique non linéaire
Date :
22/06/00
Auteur(s) :
C. DURAND
Clé :
R5.02.02-A
Page :
7/16
la formulation faible de l'équation de la chaleur peut alors s'écrire :
d

T

(T).v d
(T) T
. v
d - (T)
.v d = r
.v d

+


éq 3-3
dt
n
vol




On en déduit la formulation variationnelle du problème :
d(T).v d + (T) T .v d+ hT.vd



3 =
dt


3
r .v d + f .v d
éq 3-4
2 +
h T .v d



3 +
vol
ext

2
3
g.v d
4
4
27315
27315
4 +


[.(T + . ) -(T + . ) ].vd5
4
5
Manuel de Référence
Fascicule R5.02 : Thermique
HI-75/99/013 - Ind A

Code_Aster ®
Version
5.0
Titre :
Thermique non linéaire
Date :
22/06/00
Auteur(s) :
C. DURAND
Clé :
R5.02.02-A
Page :
8/16
4
Discrétisation en temps de l'équation différentielle
4.1
Introduction de la -méthode
Une façon classique de discrétiser une équation différentielle du premier ordre est la -méthode.
Considérons l'équation différentielle suivante :
!y(t) = (t, y(t))

éq 4.1-1
y( )

0 = y0
La -méthode consiste à discrétiser l'équation [éq 4.1-1] par un schéma aux différences finies
1 (y - y ) = (t ,y ) + 1
n 1
+
n
n 1
+
n 1
( - ) (t , y )

éq 4.1-2
t
+
n
n
yn+1 est une approximation de y(tn )
+1 , yn étant supposée connue et est le paramètre de la
méthode, [0 ]
1
, .
Remarque :
si = 0 le schéma est dit explicite,
si
0 le schéma est dit implicite.
4.2
Application à l'équation de la chaleur
Utilisons la -méthode dans la formulation variationnelle de l'équation de la chaleur, où l'on a posé :
T + = T(r,t + t
)
T - = T(r,t)
h+ = h(r,t + t
)
h- = h(r,t)
f + = f (r,t + t
)
f - = f (r,t)
T + = T (r,t + t
) T - = T (r,t)

ext
ext
ext
ext
r + = r (r,t + t
) r - = r (r,t) T d+ = T d (r,t + t
) T d- = T d (r,t
vol
vol
vol
vol
)
g + = g(r,T + )
g - = g(r,T - )
T d (r, t) représente la température imposée sur la frontière du domaine, en fonction du temps et
de l'espace.
Introduisons les espaces de fonctions suivants :
V
1
d
=

()
/ =
+
t+
{v H
v
T
1
}
V
1
d
=

()
/ =
-
t-
{v H
v
T
1
}
V
1
0 = {v H () v/ = 0
1
}
Manuel de Référence
Fascicule R5.02 : Thermique
HI-75/99/013 - Ind A

Code_Aster ®
Version
5.0
Titre :
Thermique non linéaire
Date :
22/06/00
Auteur(s) :
C. DURAND
Clé :
R5.02.02-A
Page :
9/16
Le champ T -
V

+
v V
t- étant supposé connu, on cherche T
Vt+ tel que 0 :
(T + ) - (T - ) v d
+
( (T + ) T +

. v
+ (1- ) (T - ) T -

. v
) d



t



-
( f + + (1 - ) f - ) v d2
-
( h+ T + + (
-
-
+
+
-
-


ext
1- ) h Text - h T - (1- ) h T ) v d3
2
3
-
( g + + (1- ) g - ) v d
-
4
=
( r + + (

vol 1-) rvol) v d
4

éq 4.2-1
Pour ne pas alourdir excessivement l'écriture et dans la mesure où le procédé est identique aux autres
termes, on n'a pas fait figurer le terme de rayonnement dans ces équations (intégrale sur 5 ).
En posant :
(hT ) = h+T + + (1 - ) h-T -
ext
ext
ext
f = f + + (1 - ) f -
r = r + + (1 - ) r +
vol
vol
,
on obtient finalement :
(T + ) v d + (T+) T+

. v
d + h+ T + v d




t

3


3
- g(T + ) v d4 = L (v,T - )

1
éq 4.2-2
4
v V
0
où on a posé :
(T - )
L (v,T - ) =
v d - 1
( - ) (T - ) T -

. v d + f v d
1




t
2



2

éq 4.2-3
+ ((hT ) - 1
( - ) h-T - ) v d + r v d + 1
( - ) g(T - ) v d
ext
3


4

3


4

A un instant de calcul donné, ce terme est connu. En effet, seule la température à l'instant précédent,
T - , ainsi que les valeurs à l'instant courant de fonction connues du temps, interviennent.
Dans le cas où la répartition de température dans le système à l'instant initial n'est pas fournie, on
résout le problème stationnaire. Les termes d'évolution disparaissent, = 1 ; le champ de température
à l'instant initial est donné par :
(T t=0 ) T t=

0 . v
d + ht=0 T t=0 v d
0
3
-
g(T t= ) v d



4


3

4

=
f t=0 v d
0
0
0
2
+ ht= T t= v d 3
+ rt= v d


ext


éq 4.2-4
2

3


v V
0
Manuel de Référence
Fascicule R5.02 : Thermique
HI-75/99/013 - Ind A

Code_Aster ®
Version
5.0
Titre :
Thermique non linéaire
Date :
22/06/00
Auteur(s) :
C. DURAND
Clé :
R5.02.02-A
Page :
10/16
Le problème s'écrit finalement sous la forme condensée :
Soi
t - - connu tr
ouver
,

+
T
V
T
V
+ tel que



t
t

éq 4.2-5
v V

a(v,T + ) = L (v,T


-
0
1
)
5
Discrétisation spatiale et adaptation de l'algorithme de
Newton au problème

Le principe de la méthode de Newton est très détaillé dans [R5.03.01], on n'exposera ici que les
adaptations propres à l'algorithme de thermique non linéaire.
5.1 Discrétisation
spatiale
Soit Ph un découpage de l'espace , désignons par N le nombre de noeuds du maillage, pi la
fonction de forme associé au noeud i . On désigne par J l'ensemble des noeuds appartenant à la
frontière 1.
Soient :
V h
d
=
=
( )
;
=
+
( , )

t+
{ v
v p x
v
T
x t
j
J
i
i
j
j
}
i= ,
1 N
V h
d
=
=
( )
;
=
-
( , )

éq 5.1-1
t-
{ v
v p x
v
T
x t
j
J
i
i
j
j
}
i= ,
1 N
V h
0 = { v = v p (x)
;
v = 0
j J
i
i
j
}
i= ,
1 N
Le problème [éq 4.2-5] peut être remplacé par le problème discrétisé à nombre d'inconnues fini
suivant :
Soi
t -
h
- connu tr
ouver
,

+
T
V
T
V h
+ tel que



t
t

éq 5.1-2
v
V h

a(v ,T +
0
) = L (v ,T


-
h
h
1
h
)
qu'on peut aussi écrire, avec le même formalisme que STAT_NON_LINE [R5.03.01], sous forme
vectorielle :

vT R(T+,t+ ) = vT L(T-,t+ ) v
tel que B
v = 0

éq 5.1-3
B T+


= Td(t + )
où l'opérateur B exprime la condition aux limites de température imposée T +
V h
t+ . Il est défini par :
0 si j J
(Bv) j = v

éq 5.1-4
j
si
j J

Le cas où l'application R est linéaire est traité par la commande THER_LINEAIRE [R5.02.01].
Manuel de Référence
Fascicule R5.02 : Thermique
HI-75/99/013 - Ind A

Code_Aster ®
Version
5.0
Titre :
Thermique non linéaire
Date :
22/06/00
Auteur(s) :
C. DURAND
Clé :
R5.02.02-A
Page :
11/16
La dualisation des conditions aux limites, détaillée dans [R3.03.01], conduit au problème non linéaire
en T+ :
R

(T+,t+ )+ BT +
= L(T-,t+)

éq 5.1-5
B T+
Td

+
=


(t )
Les inconnues sont le couple (T+ , +
) , où + représente les « multiplicateurs de Lagrange » des
conditions aux limites de Dirichlet.
Résoudre le système [éq 5.1-5] revient à annuler en (T+ , +
)
+
+
i
i
le vecteur F(T , ) , appelé résidu,
défini par :
L(T- , +
t ) - R(T+ , +
t ) - BT +

F(T+ , +
) =


Td
éq 5.1-6
( +
t ) - B T+

La méthode de Newton consiste à construire une suite de vecteurs {xn} convergeant vers la solution
n
de F(x) = 0 à l'aide de l'application linéaire tangente de F .
5.2 Calcul
stationnaire
Le problème variationnel est celui de l'équation [éq 4.2-4]. A noter : en calcul stationnaire, l'enthalpie
n'intervient pas dans l'application R .
On introduit la matrice de l'application linéaire tangente de la fonction R(Tn ) :
R
K n = T Tn
Celle de la fonction F(Tn , n
) est alors :
K n

BT


B
0



Dans le cas du calcul stationnaire, on doit itérer à partir d'une valeur uniforme d'initialisation du champ
de température ; en l'occurrence T0 = 0 en tout noeud. La première itération du calcul, dite itération
de prédiction, consiste à résoudre le système suivant :
K
(T ) BT T
- T L - R(T ) - BT
0
1
0
0
0


=

éq 5.2-1
B
0
d
1 - 0


T

- BT



0


Comme on peut le voir dans l'équation du problème stationnaire [éq 4.2-4], la température n'apparaît
pas au second membre : on écrit L et non L(T )
0 .
Si le problème est linéaire, R(T ) = K(T ) T = K .T
0
0
0
0 . Tous les termes en T0 disparaissent par
simplification. La solution est obtenue en une itération par inversion d'un système identique à celui
décrit dans [R5.02.01 §6].
Manuel de Référence
Fascicule R5.02 : Thermique
HI-75/99/013 - Ind A

Code_Aster ®
Version
5.0
Titre :
Thermique non linéaire
Date :
22/06/00
Auteur(s) :
C. DURAND
Clé :
R5.02.02-A
Page :
12/16
Les itérations suivantes sont des itérations de Newton, avec réactualisation ou non de la matrice
tangente K .
K
(T ) BT T
+ - T
T
i 1
i
L - R(T ) - B
(i)
i
i


=

éq 5.2-2
B
0 i
1
+ - i

0





Pour l'itération de prédiction, l'écriture du sous-système inférieur de l'équation [éq 5.2-1], après
simplification, nous assure que B T
T
1 =
d . Le premier itéré et tous les suivants vérifient donc les
conditions de Dirichlet.
Les parenthèses autour de l'indice d'itération dans l'expression K(T )
(i) signifient qu'on peut
réactualiser ou non la matrice tangente au fil des itérations.
Remarque :
La température d'initialisation T0 n'a d'influence que pour un calcul stationnaire non linéaire. En
étant de l'ordre de grandeur des températures attendues, elle permettrait de « partir » de moins
loin de la solution qu'un champ nul partout ; et ainsi diminuerait le nombre d'itérations. Aujourd'hui,
il n'est pas possible d'entrer une valeur de
T0 . Le vecteur température est initialisé, en dur, à zéro.
5.3 Calcul
transitoire
Pour la première itération du pas de temps, dite itération de prédiction, on « fait comme si » le
problème décrit dans [éq 5.1-5] était linéaire. Cette formulation doit permettre d'obtenir directement la
solution pour un problème de thermique linéaire. Mais ici, la situation est un peu différente du calcul
stationnaire à cause de la formulation en enthalpie. La linéarisation de [éq 5.1-5] donne :
R

(T- ,t+) + K(T- ,t+)(T+ - T- ) + BT +
= L(T-,t+)

éq 5.3-1
B T+
Td +


=
(t )
Ce qui revient à résoudre, pour le problème présenté sous forme matricielle :
K
(T- ) BT T+
1
L
(T- , +
t ) + K(T- )T- - R(T- )

=
d

B
+
0

éq 5.3-2
T ( +
t )

1

La fonction enthalpie est connue à une constante d'intégration près qui apparaît dans la relation liant
R(T- )
-
-
-
+
à K(T )T . Cette même constante se retrouve dans l'expression de L(T , t ) . On peut
alors l'éliminer en aboutissant au système d'équations suivant :
K
-
T
+
~
(T ) B T1
L
(T- , +
t )

= d

B
+
0

éq 5.3-3
T ( +
t )

1

~
L(T- , t + ) est le second membre calculé avec la capacité calorifique et non l'enthalpie (option
CHAR_THER_EVOLNI [§6.2]).
Enfin, comme pour le cas stationnaire vu au chapitre précédent, les itérations suivantes sont des
itérations de Newton :
Manuel de Référence
Fascicule R5.02 : Thermique
HI-75/99/013 - Ind A

Code_Aster ®
Version
5.0
Titre :
Thermique non linéaire
Date :
22/06/00
Auteur(s) :
C. DURAND
Clé :
R5.02.02-A
Page :
13/16
K(T ,t + ) BT T+
1 - T+ L(T- ,t+ ) - R(T ,t +

) - BT
(i)
i+
i
i
i


=

éq 5.3-4

B
+
+

0 i+1 - i
0


Cette fois-ci, par contre, L(T- , t + ) est calculé avec l'enthalpie et non la capacité calorifique pour être
cohérent avec R(T+ )
i
.
5.4 Convergence
Puisque le temps intervient dans l'expression de la matrice tangente, et également le pas de temps, on
préfère actualiser systématiquement celle-ci au début de chaque pas pour ne pas trop dégrader la
vitesse de convergence. En revanche, liberté est laissée à l'utilisateur de contrôler sa fréquence de
calcul au cours d'un pas de temps.
A chaque itération, on peut effectuer la recherche d'un pas de progression optimum vers la solution par
quelques itérations (2 ou 3) de recherche linéaire. Cette méthode est décrite en détail dans [R5.03.01].
Le calcul est réputé convergé quand le vecteur résidu est nul [éq 5.1-6] :
L(T- ,t + ) - R(T+ +
T +

i , t
) - B
i
F(T+ + +



i ,
i , t
) =
éq 5.4-1
Td (t + )

- B T+

i

La partie inférieure du vecteur est toujours nulle (conditions de Dirichlet). On vérifie donc :
L(T- , +
t ) - R(T+
T
i , +
t ) - B
+
i 2
éq 5.4-2
L(T- , +
t ) - BT +
i 2
L'utilisateur a également la possibilité d'arrêter les itérations sur un critère absolu :
L(T- , +
t ) - R(T+
T
i , +
t ) - B
+i


éq 5.4-3

Manuel de Référence
Fascicule R5.02 : Thermique
HI-75/99/013 - Ind A

Code_Aster ®
Version
5.0
Titre :
Thermique non linéaire
Date :
22/06/00
Auteur(s) :
C. DURAND
Clé :
R5.02.02-A
Page :
14/16
6
Principales options de thermique non-linéaire calculées
dans le Code_Aster

6.1
Conditions aux limites et chargements
On se reportera à [R5.02.01] pour les conditions aux limites et les chargements linéaires.
(1 - ) ( - )
Flux non
g T v d
CHAR_THER_FLUNL
4
linéaire
4


[(T + .
273 )
15 4 - (
-
4

1 - )(T +
.
273 )
15 ]
Rayonnement
v d
CHAR_THER_RAYO_R
4
CHAR_THER_RAYO_F
4

6.2
Calcul des matrices élémentaires et terme transitoire
Inertie
Cp
thermique,
MTAN_RIGI_MASS
v.v d +
(T ) v . v
d


+
conductivité
t



4 +
Rayonnement
(T + 27315 3
. ) v.v d


MTAN_THER_RAYO_R
4
MTAN_THER_RAYO_F
4

hv.v d
Coefficient

MTAN_THER_COEF_R
4
d'échange
MTAN_THER_COEF_F
4
dg
Flux non
MTAN_THER_FLUXNL
-
+
(T )v.v d4
linéaire
dT
4
Terme
(T - )
transitoire
CHAR_THER_EVOLNI
v d - (1 - ) (T - ) T - . v d



t



CpT - v d - (1-) (T-) T-
. v
d



t



6.3
Calcul du résidu
1
i
i
i
RESI_RIGI_MASS
(T ) v d + (( T )) T . v
d


t


i
Rayonnement
(T +
. ) v d

27315 4

RESI_THER_RAYO_R
4
RESI_THER_RAYO_F
4

( h+ Ti ) v d
Coefficient


RESI_THER_COEF_R
3
d'échange
RESI_THER_COEF_F
3
Flux non
RESI_THER_FLUXNL
-
i
(
)

linéaire
g T v d 3
3
Manuel de Référence
Fascicule R5.02 : Thermique
HI-75/99/013 - Ind A

Code_Aster ®
Version
5.0
Titre :
Thermique non linéaire
Date :
22/06/00
Auteur(s) :
C. DURAND
Clé :
R5.02.02-A
Page :
15/16
7 Bibliographie
[1]
SALENCON. Mécanique des milieux continus. Ellipses. 1988.
[2]
RUUP, PENIGUEL. SYRTHES - Conduction et rayonnement, Manuel théorique de la version
3.1. HE-41/98/048/A
Manuel de Référence
Fascicule R5.02 : Thermique
HI-75/99/013 - Ind A

Code_Aster ®
Version
5.0
Titre :
Thermique non linéaire
Date :
22/06/00
Auteur(s) :
C. DURAND
Clé :
R5.02.02-A
Page :
16/16
Page laissée intentionnellement blanche.
Manuel de Référence
Fascicule R5.02 : Thermique
HI-75/99/013 - Ind A

Document Outline