Thermal analysis

Code_Aster, Salome-Meca course material

GNU FDL licence (http://www.gnu.org/copyleft/fdl.html)
Code_Aster : main features for thermics

Frequently used as a prerequisite to a mechanical calculation

Example : laser heating of a metal disk

Thermic calculation

Laser heating of a disk : isothermal values

Mechanical calculation

Final circumferential isostresses on the deformed geometry

2 - Code_Aster and Salome-Meca course material GNU FDL Licence
A simple example : pipe in thermal equilibrium
Infinite pipe

Load : a temperature field on each side

40°C

We want to know the temperature field

in the thickness of the pipe

15°C

3 - Code_Aster and Salome-Meca course material GNU FDL Licence

A simple example : pipe in thermal equilibrium
DEBUT();

mymesh=LIRE_MAILLAGE();

mymodl=AFFE_MODELE( MAILLAGE= mymesh,

AFFE=_F( TOUT='OUI',
PHENOMENE='THERMIQUE',
MODELISATION='AXIS',),);

thmat =DEFI_MATERIAU(THER=_F(LAMBDA=6,),);

mymat=AFFE_MATERIAU( MAILLAGE= mymesh,

AFFE=_F( TOUT='OUI',
MATER=thmat,),);

myload=AFFE_CHAR_THER( MODELE= mymodl,

TEMP_IMPO=(_F(GROUP_NO='SEXTERI',
TEMP=15,),
_F(GROUP_NO='SINTERI',
TEMP=40,),),);

resuth=THER_LINEAIRE( MODELE= mymodl,

CHAM_MATER= mymat,
EXCIT=_F(CHARGE= myload,),);

IMPR_RESU( FORMAT='MED',
RESU=_F(RESULTAT= resuth,),);

FIN();

4 - Code_Aster and Salome-Meca course material GNU FDL Licence

What is the problem to be solved?

AFFE_MODELE ( ...PHENOMENE = ‘THERMIQUE’... )

Heat equation : ρCT& − λ∆T − s = 0

Linear (THER_LINEAIRE) Stationnary

OR AND OR
Non-linear (THER_NON_LINE) Transient
If : Material parameters depend on T
AND/OR
Non linear boundary conditions (flux(T))

Temperature dependence? Time dependence?

5 - Code_Aster and Salome-Meca course material GNU FDL Licence

 Available finite elements

AFFE_MODELE ( … PHENOMENE = ‘THERMIQUE’

… MODELISATION = ‘XXXX’ ...) )
3D / AXIS / PLAN
3D_DIAG / AXIS_DIAG / PLAN_DIAG
regularization during heat shock
splitting of quadratic elements in 2D
splitting not avalaible in 3D
T 289°C
Example of a sudden cooling : 20°C

temps

COQUE / COQUE_PLAN / COQUE_AXIS (linear)

thermal shell "thin" structural elements
temperature field in the thickness : 3 degrees of freedom (Tsup, Tinf and Tmiddle)

6 - Code_Aster and Salome-Meca course material GNU FDL Licence

Material behavior
DEFI_MATERIAU
At least these two characteristics: Thermal conductivity λ
Heat capacity
ρCP
3 main materials :

THER Linear isotropic

THER_ORTH Linear orthotropic (definition of λ in 3 directions)

THER_NL Non-linear behaviour : one has to define ρCP (T ) or β (T ) and λ(T )

+ other specific materials

For concrete : drying, hydratation ... etc
For metals : metallurgical transformations, hardness ... etc

7 - Code_Aster and Salome-Meca course material GNU FDL Licence

Boundary conditions and loadings (1)

AFFE_CHAR_THER (...)

AFFE_CHAR_THER_F (...)

Boundary conditions (Dirichlet)

Imposed temperatures function of time and space TEMP_IMPO

LIAISON_DDL
Linear relationships between the nodal temperatures LIAISON_GROUP
LIAISON_MAIL

8 - Code_Aster and Salome-Meca course material GNU FDL Licence

Boundary conditions and loadings (2)
Loadings (Neumann)

dT
Natural convection (Fourier law) ECHANGE λ(T ) = h(t ) ⋅ (Text − T )
dn
dT1
Heat exchange between walls ECHANGE_PAROI λ1 = h(T2 − T1 )
dn1
Normal imposed flux : constant or dT
function of time and space FLUX_REP λ(T) = f(t, x)
dn

Non linear normal flux : function of FLUX_NL dT

the temperature RAYONNEMENT λ(T) = f(T)
dn
Non-linear only

Heat source SOURCE s(x, t)

9 - Code_Aster and Salome-Meca course material GNU FDL Licence

Is the problem well defined?
For a stationnary calculation : YES

For a transient calculation : NO

A condition is missing : the temperature at the initial time

Defined in THER_LINEAIRE or THER_NON_LINE with ETAT_INIT

4 possibilities :

Résults of a stationnary calculation STATIONNAIRE ='OUI'

A constant temperature VALE = T0

A known temperature field CHAM_NO = …

for example, created with CREA_CHAMP

The result of another thermal calculation A result concept (EVOL_THER = resu)

+ chosen time (NUME_ORDRE or INST )

10 - Code_Aster and Salome-Meca course material GNU FDL Licence

Solving methods

Linear thermics : THER_LINEAIRE

Non-linear thermics : THER_NON_LINE

Newton method for non-linéarities
Enthalpy formulation β (T ) Numerically easier to handle
phase changes

ρCP (T ) β (T )

T T
11 - Code_Aster and Salome-Meca course material GNU FDL Licence
Precautions for use in thermics
Time discretization with the θ-méthode.
ρC
(Tn+1 − Tn )=θ ∆Tn+1+(1−θ )∆Tn
λ.∆t
θ = 0.57, good compromise and the default value
Ability to change θ with PARM_THETA

Beware to inconsistent refinements in time and space.

Respect the time and distance characteristics of the material

ρC ( ∆x )
2
∆t ≤
(6θ )λ
Drive the non-linear convergence, if necessary.
Frequency of updating the tangent matrix
NEWTON
Linear search criteria
Convergence criteria (number of iterations, residues ...) CONVERGENCE

12 - Code_Aster and Salome-Meca course material GNU FDL Licence

Post-processing options

What's the result of THER_LINEAIRE / THER_NON_LINE ?

By default : only the nodal temperatures

How to get the thermal fluxes ?

With the operator THER_NON_LINE ( … OPTION … )

FLUX_ELNO (at nodes, by elements)
FLUX_ELGA (at the Gauss points)

With CALC_CHAMP
FLUX_NOEU (at nodes)

13 - Code_Aster and Salome-Meca course material GNU FDL Licence

Metallurgy

Possibilities :
Compute in moving frame (welding) : THER_NON_LINE_MO
With forced convection (welding) AFFE_CHAR_THER (...CONVECTION...)
Calculate the hardness, the metallurgical transformations CALC_META

14 - Code_Aster and Salome-Meca course material GNU FDL Licence

Thermo-hydration and drying
Common use : chained calculation :
Thermic and hydration calculation of concrete
Temperature field Hydration field

Calculation of concrete drying

Water concentration field

Mechanical calculation of
concrete shrinkage
Displacements, strains …
Stresses field

15 - Code_Aster and Salome-Meca course material GNU FDL Licence

And mechanics?

How to take into account the temperature field in the mechanical

calculation?

In AFFE_MATERIAU (...
AFFE_VARC=_F( NOM_VARC = 'TEMP',
EVOL = EVOTH,
VALE_REF = 20. )
...)
The mesh to be used in mechanics must be different from that used for
thermics. For two main reasons:
A linear mesh is better in thermics, while rather quadratic in mechanics
Areas of interest are different => the refined areas are different

One must project the temperature field on the mesh used for mechanics

PROJ_CHAMP ( … )

16 - Code_Aster and Salome-Meca course material GNU FDL Licence

And mechanics?
Example of the cylinder
DEBUT()
thmesh=LIRE_MAILLAGE()
thmodl=AFFE_MODELE(...)
thmater=DEFI_MATERIAU(...)
fthmater=AFFE_MATERIAU(...)
thload=AFFE_CHAR_THER(...)
thresu=THER_LINEAIRE(...)

Projection sur le
nouveau maillage
We want to make a mechanical calculation with a different mesh

memesh=LIRE_MAILLAGE()
thmodl2=AFFE_MODELE(MAILLAGE= memesh,
AFFE=_F( TOUT='OUI',
PHENOMENE='THERMIQUE',
MODELISATION='AXIS',),)
resproj=PROJ_CHAMP( RESULTAT= thresu,
MODELE_1= thmodl,
MODELE_2= thmodl2,)
...

17 - Code_Aster and Salome-Meca course material GNU FDL Licence

And mechanics?
Mechanical calculation of the expansion of the cylinder
mymodl=AFFE_MODELE( MAILLAGE=mymesh,
AFFE=_F( TOUT='OUI',
PHENOMENE='MECANIQUE',
MODELISATION='AXIS',),)
Temperature field
steel=DEFI_MATERIAU( ELAS=_F( E=2.1e11, + reference value
NU=0.2,
ALPHA=12e-6,),)

mymat=AFFE_MATERIAU( MAILLAGE=mymesh,
AFFE=_F( TOUT='OUI',
MATER=steel,),
AFFE_VARC=_F(TOUT='OUI',
NOM_VARC='TEMP',
EVOL=projres,
VALE_REF=20,),)

mecload=AFFE_CHAR_MECA(MODELE=MODME,
DDL_IMPO=_F(GROUP_NO='FACEB',
DX=0,
DY=0,),)

mecres=MECA_STATIQUE(MODELE=mymodl,
CHAM_MATER=mymat,
EXCIT=_F(CHARGE= mecload,),)

IMPR_RESU(...)

FIN()
18 - Code_Aster and Salome-Meca course material GNU FDL Licence
Analytical temperature field
Here we know the analytical formula of the temperature field in the thickness of the cylinder
(X):  X  T −T
Tinte − ln  exte inte
 Rinte  ln Rexte 
DEBUT(); R 
// Thermic calculation  inte 
thmesh=LIRE_MAILLAGE() memesh=LIRE_MAILLAGE()
thmodl=AFFE_MODELE(...)
T0 = FORMULE(VALE='40-log(X/20)*(40-15)/log(21/20)',
thmater=DEFI_MATERIAU(...) NOM_PARA=('X',),)
fthmater=AFFE_MATERIAU(...) fT0 =CREA_CHAMP( TYPE_CHAM='NOEU_TEMP_F',
thload=AFFE_CHAR_THER(...) OPERATION='AFFE',
thresu=THER_LINEAIRE(...) MAILLAGE= memesh,
// Projection AFFE=_F( TOUT='OUI',
memesh=LIRE_MAILLAGE() NOM_CMP='TEMP',
thmodl2=AFFE_MODELE(...) VALE_F=T0,),)
resproj=PROJ_CHAMP(...)
thresu=CREA_RESU(OPERATION='AFFE',
TYPE_RESU='EVOL_THER',
// Mechanic calculation NOM_CHAM='TEMP',
memodl=AFFE_MODELE(...) AFFE=_F(CHAM_GD=fT0,
steel=DEFI_MATERIAU(...) INST=0,),)
memat=AFFE_MATERIAU( MAILLAGE=memesh,
AFFE=_F( TOUT='OUI',
MATER=steel,),
AFFE_VARC=_F(TOUT='OUI',
NOM_VARC='TEMP',
EVOL= resproj,
VALE_REF=20,),)
meload=AFFE_CHAR_MECA(...)
meresu=MECA_STATIQUE(...)
IMPR_RESU(...)
FIN()

19 - Code_Aster and Salome-Meca course material GNU FDL Licence

Example of use

Laser treatment of a disc in linear thermics

Flux(r) = A + exp ( -r²/S²)

during 30 secondes and extinction

L_ECHAN

LBC (40 mm) T(t=0) = 20°C

L_ECHAN
Exchange by natural convection :
Text = 20°C
h = 5 W/m²°C

20 - Code_Aster and Salome-Meca course material GNU FDL Licence

Example of use (2)

Model and material definitions

DEBUT( )
mymesh = LIRE_MAILLAGE ( )
# Definition of the model
mymodl = AFFE_MODELE(MAILLAGE = mymesh ,
AFFE = _F( TOUT = 'OUI',
PHENOMENE = 'THERMIQUE',
MODELISATION = 'AXIS', ) )
# Definition of the thermal properties of the material
steel = DEFI_MATERIAU( THER = _F ( RHO_CP = 3.6e6,
LAMBDA = 14, ) )

# ssignment of the material on the finite element model

mymat= AFFE_MATERIAU(MAILLAGE= mymesh,
AFFE=_F( TOUT = 'OUI',
MATER = steel , ) )

21 - Code_Aster and Salome-Meca course material GNU FDL Licence

Example of use (3)
Definition of data needed for the loadings definition
# "mathematical" definition of the F(R) formula
A = 5e6
S = 0.0138
F = FORMULE(NOM_PARA='X',VALE="(A + EXP(-X**2/S**2))")

# Interpretation of the F(R) formula and transformation in a function Creation of a flux

# (tabulation X ; F(X)) function of X (=R) in
order to represent
the heat due to the
# Definition of the abscissa for tabulation laser
LISTR=DEFI_LIST_REEL(DEBUT = 0.,
INTERVALLE=_F( JUSQU_A = 0.04, NOMBRE = 100, ) ,)
# Calculation (and tabulation) of the formula for the 'LISTR' values of the radius
FF = CALC_FONC_INTERP( FONCTION = F,
LIST_PARA = LISTR,
PROL_DROITE ='LINEAIRE',
PROL_GAUCHE ='LINEAIRE' )

# Temperature of the outside

T_ext = DEFI_CONSTANTE ( VALE=20., )
# Coefficient for exchange by convection
H = DEFI_CONSTANTE(VALE=5 ,)
22 - Code_Aster and Salome-Meca course material GNU FDL Licence
Example of use (4)
Definition of two loadings
# Exchange condition on the edge (during heating)
thload1 = AFFE_CHAR_THER_F( MODELE=mymodl,
FLUX_REP =_F( GROUP_MA = 'LBC', FLUN = FF,),
ECHANGE =_F( GROUP_MA = ('BAS','HAUT','COTE',), HAUT
COEF_H = H, LBC
COTE
TEMP_EXT = T_EXT, ) )
BAS

# Exchange condition on the edge (during the free cooling) : convection everywhere
thload2 = AFFE_CHAR_THER_F ( MODELE = mymodl,
ECHANGE=_F( GROUP_MA =('LBC','HAUT','COTE','BAS',),
COEF_H = H ,
HAUT
TEMP_EXT = T_EXT , ),)
LBC
COTE

BAS

23 - Code_Aster and Salome-Meca course material GNU FDL Licence

Example of use (5)
Calculation
# Definition of time steps
list = DEFI_LIST_REEL( DEBUT=0., INTERVALLE=_F( JUSQU_A = 150., NOMBRE = 300, ) )

# Calculation for the heating phase

temper =THER_LINEAIRE( MODELE= mymodl,
CHAM_MATER = mymat,
EXCIT = _F( CHARGE = thload1, ),
INCREMENT = _F( LIST_INST = LIST,
INST_FIN = 30, ), Stop at t=30s
ETAT_INIT = _F( VALE = 20.0, ) ) at t=0s : T=20°C everywhere

# Calculation for the cooling phase

temper = THER_LINEAIRE( reuse = temper ,
MODELE = mymodl,
CHAM_MATER = mymat,
PARM_THETA = 0.60,
EXCIT = _F( CHARGE = thload2, ),
INCREMENT = _F( LIST_INST = list ,
INST_INIT = 30,), at t=30s : T= T° at the end
ETAT_INIT = _F( EVOL_THER = temper ,)) of the previous calculation
IMPR_RESU( ... )
FIN ( )

24 - Code_Aster and Salome-Meca course material GNU FDL Licence

Results

During the heating t=4s

at t=30s

During the cooling t=60s

25 - Code_Aster and Salome-Meca course material GNU FDL Licence

 End of presentation

Is something missing or unclear in this document?

Or feeling happy to have read such a clear tutorial?

Please, we welcome any feedbacks about Code_Aster training materials.

Do not hesitate to share with us your comments on the Code_Aster forum
dedicated thread.

26 - Code_Aster and Salome-Meca course material GNU FDL Licence