Equilibre des solides thermoélastiques dans le

contexte infinitésimal
 Plan

1 Loi de comportement thermoélastique linéarisée

Cas des déformations infinitésimales
Cas des transformations infinitésimales
Solides thermoélastiques linéaires isotropes
Module de Young et coefficient de Poisson
Anisotropie
Couplage thermoélastique

2 Le problème aux limites d’équilibre thermoélastique linéarisé

Linéarisation des équations d’équilibre
Formulation du problème aux limites
Schémas de résolution
Existence et unicité
Stabilité

3 Bilan
 Plan

1 Loi de comportement thermoélastique linéarisée

Cas des déformations infinitésimales
Cas des transformations infinitésimales
Solides thermoélastiques linéaires isotropes
Module de Young et coefficient de Poisson
Anisotropie
Couplage thermoélastique

2 Le problème aux limites d’équilibre thermoélastique linéarisé

Linéarisation des équations d’équilibre
Formulation du problème aux limites
Schémas de résolution
Existence et unicité
Stabilité

3 Bilan
 Plan

1 Loi de comportement thermoélastique linéarisée

Cas des déformations infinitésimales
Cas des transformations infinitésimales
Solides thermoélastiques linéaires isotropes
Module de Young et coefficient de Poisson
Anisotropie
Couplage thermoélastique

2 Le problème aux limites d’équilibre thermoélastique linéarisé

Linéarisation des équations d’équilibre
Formulation du problème aux limites
Schémas de résolution
Existence et unicité
Stabilité

3 Bilan
Thermoélasticité linéarisée en petites déformations
• contexte des déformations thermoélastiques infinitésimales
| T − T0 |
kE
∼
k  1, 1
T0
les grandes rotations sont autorisées!
• potentiel thermoélastique quadratique
1 ρ0 Cε
ρ0 ψ0 (E
∼
, T ) = −ρ0 (T − T0 )ŝ 0 + Π ∼
0
:E∼
− (T − T0 )2
2 T0
1
− (T − T0 )P ∼
:E ∼
+ E :Λ:E
2∼ ∼ ∼
∼

• lois d’état linéaires (affines)

∂ψ0 0
Π∼
= ρ0 =Π ∼
+Λ ∼
:E
∼
−P ∼
(T − T0 )
∂E∼
∼

∂ψ0 ρ0 Cε
ρ0 s0 = −ρ0 = ρ0 ŝ 0 + (T − T0 ) + P ∼
:E
∼
∂T T0

Loi de comportement thermoélastique linéarisée 5/64

Thermoélasticité linéarisée en petites déformations
• contexte des déformations thermoélastiques infinitésimales
| T − T0 |
kE
∼
k  1, 1
T0
• lois d’état linéaires (affines)
∂ψ0 0
Π∼
= ρ0 =Π ∼
+Λ ∼
:E
∼
−P ∼
(T − T0 )
∂E∼
∼

∂ψ0 ρ0 Cε
ρ0 s0 = −ρ0 = ρ0 ŝ 0 + (T − T0 ) + P
∼
:E
∼
∂T T0
Λ
∼
tenseur des modules d’élasticité
∼
Cε chaleur massique à déformation constante
−1
α
∼
=Λ∼
:P
∼
tenseur des dilatations thermiques
∼
0
Π
∼
état de contrainte initial au point X
• déformations d’origine thermique
th
E
∼
= (T − T0 )α
∼
0 th
Π
∼
=Π
∼
+Λ
∼
: (E
∼
−E
∼
)
∼
Loi de comportement thermoélastique linéarisée 6/64
 Tenseur des modules d’élasticité
∂Πij
• rigidités élastiques d’un corps matériel Λijkl := ∂Ekl
∂Π
∼
Λ (E, T ) :=
∼ ∼
(E
∼
, T ) = Λijkl (E
∼
,T)ei ⊗ ej ⊗ ek ⊗ el
∼ ∂E
∼

tenseur euclidien d’ordre 4

• symétries mineures du tenseur des modules résultant de la
symétrie des tenseurs Π
∼
et E
∼

Λijkl = Λijlk = Λjikl

• symétries majeures du tenseur des modules résultant de
l’existence du potentiel d’élasticité
∂ 2 ψ0
 
∂ ∂ψ0 ∂Πkl
Λijkl = ρ0 = = ρ0 = Λklij
∂Ekl ∂Eij ∂Eij ∂Ekl ∂Eij
• tenseur des souplesses
th 0 −1
E
∼
=E
∼
+S
∼
: (Π
∼
−Π
∼
), S
∼
=Λ
∼
∼ ∼ ∼

Loi de comportement thermoélastique linéarisée 7/64

Thermoélasticité linéarisée en petites déformations
cas des matériaux isotropes
2
Π
∼
= α0 (IE , T )1
∼
+ α1 (IE , T )E
∼
+ α2 (IE , T )E
∼

IE symbolise les invariants principaux de E

∼

Loi de comportement thermoélastique linéarisée 8/64

Thermoélasticité linéarisée en petites déformations
cas des matériaux isotropes
2
Π
∼
= α0 (IE , T )1
∼
+ α1 (IE , T )E
∼
+ α2 (IE , T )E
∼

IE symbolise les invariants principaux de E

∼

Π
∼
= Π0 1
∼
+ λ(trace E)1 + 2µE
∼ ∼ ∼
− P(T − T0 )1
∼

λ, µ coefficients de Lamé (MPa)

Loi de comportement thermoélastique linéarisée 9/64

 Plan

1 Loi de comportement thermoélastique linéarisée

Cas des déformations infinitésimales
Cas des transformations infinitésimales
Solides thermoélastiques linéaires isotropes
Module de Young et coefficient de Poisson
Anisotropie
Couplage thermoélastique

2 Le problème aux limites d’équilibre thermoélastique linéarisé

Linéarisation des équations d’équilibre
Formulation du problème aux limites
Schémas de résolution
Existence et unicité
Stabilité

3 Bilan
 Thermoélasticité linéarisée en petites
transformations
• le contexte des transformations infinitésimales
H
∼
= Grad u = F
∼
−1
∼
=ε
∼
+ω
∼

| T − T0 |
kH
∼
k  1, 1
T0
kε
∼
k  1, kω ∼
k1
les grandes rotations ne sont pas autorisées!
1 ∂ui ∂uj ∂uk ∂uk
Eij = ( + + )
2 ∂Xj ∂Xi ∂Xi ∂Xj
E
∼
'ε
∼

Loi de comportement thermoélastique linéarisée 11/64

 Thermoélasticité linéarisée en petites
transformations
• le contexte des transformations infinitésimales
H
∼
= Grad u = F
∼
−1
∼
=ε
∼
+ω
∼

| T − T0 |
kH
∼
k  1, 1
T0
kε
∼
k  1, kω ∼
k1
les grandes rotations ne sont pas autorisées!
1 ∂ui ∂uj ∂uk ∂uk
Eij = ( + + )
2 ∂Xj ∂Xi ∂Xi ∂Xj
E
∼
'ε
∼

• solides thermoélastiques linéaires

0
Π
∼
=Π
∼
+Λ
∼
:ε
∼
− (T − T0 )P
∼
∼
0 th th
Π
∼
=Π
∼
+Λ
∼
: (ε
∼
−ε
∼
), ε
∼
= (T − T0 )α
∼
∼

Loi de comportement thermoélastique linéarisée 12/64

 Thermoélasticité linéarisée en petites
transformations
• retour sur le tenseur des contraintes de Cauchy
ρ ρ
σ∼
= F .Π.FT = (1
∼ ∼ ∼
+ H).(Π0 + Λ : (ε −ε th
)).(1 +H T
)
ρ0 ρ0 ∼ ∼ ∼ ∼
∼
∼ ∼ ∼ ∼

• linéarisation en ne gardant que les termes d’ordre 1 en H ∼

Loi de comportement thermoélastique linéarisée 13/64

 Thermoélasticité linéarisée en petites
transformations
• relation entre le tenseur des contraintes de Cauchy et le tenseur des
déformations infinitésimales
ρ ρ
σ
∼
= F .Π.FT = (1
∼ ∼ ∼
+ H).(Π0 + Λ : (ε −ε th
)).(1 +H T
)
ρ0 ρ0 ∼ ∼ ∼ ∼
∼
∼ ∼ ∼ ∼

• linéarisation en ne gardant que les termes d’ordre 1 en H

∼

th 0
σ
∼
=Λ
∼
: (ε
∼
−ε
∼
)+Π
∼
− (trace ε)Π0 + H
∼ ∼
.Π0 + Π
∼ ∼ ∼
0
.H
∼
T
∼

Loi de comportement thermoélastique linéarisée 14/64

 Thermoélasticité linéarisée en petites
transformations
• relation entre le tenseur des contraintes de Cauchy et le tenseur des
déformations infinitésimales
ρ ρ
σ
∼
= F .Π.FT = (1
∼ ∼ ∼
+ H).(Π0 + Λ : (ε −ε th
)).(1 +H T
)
ρ0 ρ0 ∼ ∼ ∼ ∼
∼
∼ ∼ ∼ ∼

• linéarisation en ne gardant que les termes d’ordre 1 en H

∼

th 0
σ
∼
=Λ
∼
: (ε
∼
−ε
∼
)+Π
∼
− (trace ε)Π0 + H
∼ ∼
.Π0 + Π
∼ ∼ ∼
0
.H
∼
T
∼

• en pratique, on a kΠ 0
k  kΛ k, σ0 ∼ 10−3
∼ ∼
∼
µ
th 0
σ
∼
=Λ
∼
: (ε
∼
−ε
∼
)+Π
∼
+ω .Π0 − Π
∼ ∼ ∼
0
.ω
∼
∼

Loi de comportement thermoélastique linéarisée 15/64

 Thermoélasticité linéarisée en petites
transformations
• relation entre le tenseur des contraintes de Cauchy et le tenseur des
déformations infinitésimales
ρ ρ
σ
∼
= F .Π.FT = (1
∼ ∼ ∼
+ H).(Π0 + Λ : (ε −ε th
)).(1 +H T
)
ρ0 ρ0 ∼ ∼ ∼ ∼
∼
∼ ∼ ∼ ∼

• linéarisation en ne gardant que les termes d’ordre 1 en H

∼

th 0
σ
∼
=Λ
∼
: (ε
∼
−ε
∼
)+Π
∼
− (trace ε)Π0 + H
∼ ∼
.Π0 + Π
∼ ∼ ∼
0
.H
∼
T
∼

• en pratique, on a kΠ 0
k  kΛ k, σ0 ∼ 10−3
∼ ∼
∼
µ
th 0
σ
∼
=Λ
∼
: (ε
∼
−ε
∼
)+Π
∼
+ω .Π0 − Π
∼ ∼ ∼
0
.ω
∼
∼
0
• linéarisation autour d’un état naturel Π
∼
=0
th
σ
∼
=Λ
∼
: (ε
∼
−ε
∼
)
∼

• linéarisation autour d’un état quasi–naturel kω .Π0 k  kΛ

∼ ∼ ∼
εk
:∼
∼

0 th
σ
∼
=Π
∼
+ Λ : (ε
∼
∼
−ε
∼
)
Loi de comportement thermoélastique linéarisée
∼ 16/64
 Plan

1 Loi de comportement thermoélastique linéarisée

Cas des déformations infinitésimales
Cas des transformations infinitésimales
Solides thermoélastiques linéaires isotropes
Module de Young et coefficient de Poisson
Anisotropie
Couplage thermoélastique

2 Le problème aux limites d’équilibre thermoélastique linéarisé

Linéarisation des équations d’équilibre
Formulation du problème aux limites
Schémas de résolution
Existence et unicité
Stabilité

3 Bilan
Solide thermoélastique isotrope linéaire dans l’état
quasi–naturel

e th
• décomposition de la déformation totale ε
∼
=ε
∼
+ε
∼
avec
th
ε
∼
= α(T − T )
0 ∼1
• loi de comportement thermoélastique linéaire isotrope (état
quasi–naturel)

σ
∼
= σ01
∼
+ λ(trace ε
∼
e
)1
∼
+ 2µ ε
∼
e

σ
∼
= σ01
∼
+ λ(trace ε) 1 + 2µ ε
∼ ∼ ∼
− 3κα(T − T0 )1
∼

avec
3λ + 2µ
κ=
3
module de compressibilité

Loi de comportement thermoélastique linéarisée 18/64

 Module de cisaillement µ
2 2

1 z 1
cisaillement simple
 γ
 u1 = 2 X2
   
0 γ/2 0 0 σ12 0


γ ] =  γ/2 0 0  = [ε ] =  σ12 0 0 
u2 = X1 [H ∼ ∼
], [σ
∼

 2 0 0 0 0 0 0
u3 = 0


σ12 = 2µ ε12 = µ γ
acier : µ = 80000 MPa, aluminium : µ = 26000 MPa
PU : µ = 400 MPa
Loi de comportement thermoélastique linéarisée 19/64
 Module de compressibilité κ

2 2

1 1
dilatation pure

u = AX, H
∼
=ε
∼
= A1
∼
, σ
∼
= 3κε
∼

trace σ
∼
= 3κ trace ε
∼

acier : κ = 165000 MPa, aluminium : κ = 70000 MPa

PU : κ = 1700 MPa, eau : κ = 2000 MPa

Loi de comportement thermoélastique linéarisée 20/64

 Les valeurs propres du tenseur des modules
d’élasticité
• Λ
∼
est un endomorphisme sur les tenseurs d’ordre 2
∼
symétriques. Il est symétrique (auto–adjoint). On peut le voir
aussi comme une forme bilinéaire symétrique :
? ? ?
σ
∼
:ε
∼
=ε
∼
:Λ
∼
:ε
∼
=ε
∼
:Λ
∼
:ε
∼
∼ ∼

• les valeurs propres du tenseur des modules d’élasticité isotrope

linéarisée Λ iso sont
∼
∼
? 3κ d’ordre 1
iso
Λ
∼
εsph = 3κε
:∼ ∼
sph
∼
? 2µ d’ordre 5
iso
Λ
∼
εdev = 2µε
:∼ ∼
dev
∼

• stabilité locale du milieu élastique autour de l’état naturel :

le tenseur des modules est défini positif
κ > 0, µ>0
(i.e. l’état naturel est observable)
Loi de comportement thermoélastique linéarisée 21/64
 Plan

1 Loi de comportement thermoélastique linéarisée

Cas des déformations infinitésimales
Cas des transformations infinitésimales
Solides thermoélastiques linéaires isotropes
Module de Young et coefficient de Poisson
Anisotropie
Couplage thermoélastique

2 Le problème aux limites d’équilibre thermoélastique linéarisé

Linéarisation des équations d’équilibre
Formulation du problème aux limites
Schémas de résolution
Existence et unicité
Stabilité

3 Bilan
Inversion de la loi de thermoélasticité linéaire
isotrope
• trace et déviateur σ = λ(trace ε e )1 + 2µεe
∼ ∼ ∼ ∼
=⇒ trace ε e = 1 trace σ , (εe )dev 1 dev
= 2µ σ
∼ 3κ ∼ ∼ ∼

e trace εe
∼ e dev trace σ∼
σ
∼
dev
ε
∼
= 1
∼
+ (ε
∼
) = 1 +
3 9κ ∼ 2µ
• relations d’élasticité linéaire isotrope (état naturel)
e 1+ν ν
ε
∼
= σ
∼
− (trace σ )1
∼ ∼
E E
E module de Young (MPa), ν coefficient de Poisson

Loi de comportement thermoélastique linéarisée 23/64

Inversion de la loi de thermoélasticité linéaire
isotrope
• trace et déviateur σ = λ(trace ε e )1 + 2µεe
∼ ∼ ∼ ∼
=⇒ trace ε e = 1 trace σ , (εe )dev 1 dev
= 2µ σ
∼ 3κ ∼ ∼ ∼

e trace εe
∼ e dev trace σ∼
σ
∼
dev
ε
∼
= 1
∼
+ (ε
∼
) = 1 +
3 9κ ∼ 2µ
• relations d’élasticité linéaire isotrope (état naturel)
e 1+ν ν
ε
∼
= σ
∼
− (trace σ )1
∼ ∼
E E
E module de Young (MPa), ν coefficient de Poisson
E ν λ
µ= , =
2(1 + ν) E 2µ(3λ + 2µ)
λ 3λ + 2µ 3λ + 2µ E
ν= , E =µ , κ= =
2(λ + µ) λ+µ 3 3(1 − 2ν)
• stabilité élastique E > 0, −1 < ν < 0.5

Loi de comportement thermoélastique linéarisée 24/64

Traction simple en élasticité linéaire isotrope

 
0 0 0
[σ
∼
]= 0 0 0 
0 0 σ
y
y  σ 
−ν 0 0
 E σ 
[ε ]= 0 −ν 0 
 
∼
 E σ 
2
x 0 0
E
1
1

Loi de comportement thermoélastique linéarisée 25/64

 • acier : E = 200000 MPa, ν = 0.3
aluminium : E = 70000 MPa, ν = 0.33,
PU E = 2000 MPa, ν = 0.4
en général 0 < ν < 0.5
• matériaux cellulaires

liège : ν ' 0

Loi de comportement thermoélastique linéarisée 26/64

 Matériaux auxétiques

coefficient de Poisson négatif (1987)!

Loi de comportement thermoélastique linéarisée 27/64

 Plan

1 Loi de comportement thermoélastique linéarisée

Cas des déformations infinitésimales
Cas des transformations infinitésimales
Solides thermoélastiques linéaires isotropes
Module de Young et coefficient de Poisson
Anisotropie
Couplage thermoélastique

2 Le problème aux limites d’équilibre thermoélastique linéarisé

Linéarisation des équations d’équilibre
Formulation du problème aux limites
Schémas de résolution
Existence et unicité
Stabilité

3 Bilan
 Notation de Voigt

représentation matricielle de l’algèbre des tenseurs d’ordre 2

symétriques
• matrice des rigidités
    
σ11 C11 C12 C13 C14 C15 C16 ε11
 σ22   C12 C22 C23 C24 C25 C26  ε22 
    
 σ33   C13 C23 C33 C34 C35 C36  ε33 
 σ23  =  C14 C24
    
   C34 C44 C45 C46 
 γ23 = 2ε23 

 σ31   C15 C25 C35 C45 C55 C56  γ31 = 2ε31 
σ12 C16 C26 C36 C46 C56 C66 γ12 = 2ε12

• le nombre de composantes indépendantes dépend du groupe

des symétries matérielles du milieu considéré
cas triclinique : 21 coefficients indépendants

Loi de comportement thermoélastique linéarisée 29/64

 Notation de Voigt

représentation matricielle de l’algèbre des tenseurs d’ordre 2

symétriques
• matrice des rigidités
    
σ11 C11 C12 C13 0 0 C16 ε11
 σ22   C12 C22 C23 0 0 C26  ε22 
    
 σ33   C13 C23 C33 0 0 C36  ε33 
 σ23  =  0
    
   0 0 C44 C45 0 
 γ23 = 2ε23 

 σ31   0 0 0 C45 C55 0  γ31 = 2ε31 
σ12 C16 C26 C36 0 0 C66 γ12 = 2ε12

• le nombre de composantes indépendantes dépend du groupe

des symétries matérielles du milieu considéré
cas monoclinique : 13 coefficients indépendants

Loi de comportement thermoélastique linéarisée 30/64

 Notation de Voigt

représentation matricielle de l’algèbre des tenseurs d’ordre 2

symétriques
• matrice des rigidités
    
σ11 C11 C12 C13 0 0 0 ε11
 σ22   C12 C22 C23 0 0 0  ε22 
    
 σ33   C13 C23 C33 0 0 0  ε33 
 σ23  =  0
    
   0 0 C44 0 0 
 γ23 = 2ε23 

 σ31   0 0 0 0 C55 0  γ31 = 2ε31 
σ12 0 0 0 0 0 C66 γ12 = 2ε12

• le nombre de composantes indépendantes dépend du groupe

des symétries matérielles du milieu considéré
cas orthotrope : 9 coefficients indépendants

Loi de comportement thermoélastique linéarisée 31/64

 Notation de Voigt

représentation matricielle de l’algèbre des tenseurs d’ordre 2

symétriques
• matrice des rigidités
    
σ11 C11 C12 C12 0 0 0 ε11
 σ22   C12 C11 C12 0 0 0  ε22 
    
 σ33   C12 C12 C11 0 0 0  ε33 
 σ23  =  0
    
   0 0 C44 0 0 
 γ23 = 2ε23 

 σ31   0 0 0 0 C44 0  γ31 = 2ε31 
σ12 0 0 0 0 0 C44 γ12 = 2ε12

• le nombre de composantes indépendantes dépend du groupe

des symétries matérielles du milieu considéré
cas isotrope : 2 coefficients C11 = λ + 2µ, C12 = λ, C44 = µ

Loi de comportement thermoélastique linéarisée 32/64

 Isotropie transverse

matériaux composites nids d’abeille

(fibre 100 µm) (taille de cellule 1 cm)

Loi de comportement thermoélastique linéarisée 33/64

 Plan

1 Loi de comportement thermoélastique linéarisée

Cas des déformations infinitésimales
Cas des transformations infinitésimales
Solides thermoélastiques linéaires isotropes
Module de Young et coefficient de Poisson
Anisotropie
Couplage thermoélastique

2 Le problème aux limites d’équilibre thermoélastique linéarisé

Linéarisation des équations d’équilibre
Formulation du problème aux limites
Schémas de résolution
Existence et unicité
Stabilité

3 Bilan
 Equation de la chaleur
• équation de la chaleur en thermoélasticité linéaire
ρT ṡ + div q = ρr
ρ0 Cε
ρ0 s = ρ0 s 0 + (T − T0 ) + 3ακ(trace ε
∼
)
T0
• conduction thermique linéarisée q = −k ∇T

ρCε Ṫ + 3καT0 (trace ε̇

∼
) −k∆T = ρr
| {z }
couplage thermomécanique
• régime permanent ∆T = 0
• régime homogène (adiabatique) ρCε Ṫ + 3καT0 (trace ε̇ ∼
)=0
traction simple :
quand on tire, le matériau se refroidit! T − T0 = − 3Eρ0αT
Cε ε33
0

• échauffements importants dans le cas de comportements

non–linéaires (viscoplasticité)
Loi de comportement thermoélastique linéarisée 35/64
 Plan

1 Loi de comportement thermoélastique linéarisée

Cas des déformations infinitésimales
Cas des transformations infinitésimales
Solides thermoélastiques linéaires isotropes
Module de Young et coefficient de Poisson
Anisotropie
Couplage thermoélastique

2 Le problème aux limites d’équilibre thermoélastique linéarisé

Linéarisation des équations d’équilibre
Formulation du problème aux limites
Schémas de résolution
Existence et unicité
Stabilité

3 Bilan
Problème d’équilibre d’un solide thermoélastique
• sollicitations (champ T (x ) donné)
? efforts volumiques

ρf (x ), x ∈ Ωt

? efforts surfaciques

∂Ωu t =σ
∼
.n = t d (x ), x ∈ ∂Ωtd

ud ? déplacements imposés

Ωt ∂Ωtd u = u d, x ∈ ∂Ωu

• équation de champ
td
div σ
∼
+ ρf = 0, x ∈ Ωt

• loi de comportement élastique

σ
∼
(x ) = FΩ0 (F
∼
(X ), T (X ))
Le problème aux limites d’équilibre thermoélastique linéarisé 37/64
 Plan

1 Loi de comportement thermoélastique linéarisée

Cas des déformations infinitésimales
Cas des transformations infinitésimales
Solides thermoélastiques linéaires isotropes
Module de Young et coefficient de Poisson
Anisotropie
Couplage thermoélastique

2 Le problème aux limites d’équilibre thermoélastique linéarisé

Linéarisation des équations d’équilibre
Formulation du problème aux limites
Schémas de résolution
Existence et unicité
Stabilité

3 Bilan
 L’hypothèse des petites perturbations (HPP)
les hypothèses présidant à la linéarisation du problème d’équilibre
sont :
• contexte infinitésimal : petite déformation+petite rotation

kGrad u k ∼ h  1

• variations limitées de température

• état initial naturel
0
Π
∼
=0
dans le cas d’un état initial quasi–naturel, les résultats sont
valides en remplaçant σ par σ −σ 0
∼ ∼ ∼

• solide thermoélastique homogène (et isotrope dans la suite,

par commodité). Si le milieu est homogène par morceaux, il
faudra écrire des conditions d’interface
Elles constituent le cadre HPP.
Le problème aux limites d’équilibre thermoélastique linéarisé 39/64
Problème d’équilibre thermoélastique linéarisé
• sollicitations (champ T (X )
donné)
? efforts surfaciques

t =σ
∼
.N = t d (X ), X ∈ ∂Ωtd

∂Ωu ? déplacements imposés

ud u = u d, X ∈ ∂Ωu

Ω0 ∂Ωtd • équation de champ

Div σ
∼
(X )+ρ0 f (X ) = 0, X ∈ Ω0
td
• loi de comportement élastique
(état naturel)

σ
∼
(X ) = Λ
∼
: (ε
∼
− (T − T0 )α
∼
)
∼
in fine, tout se passe dans les équations comme si X et x coı̈ncidaient.
Le problème aux limites d’équilibre thermoélastique linéarisé 40/64
 Plan

1 Loi de comportement thermoélastique linéarisée

Cas des déformations infinitésimales
Cas des transformations infinitésimales
Solides thermoélastiques linéaires isotropes
Module de Young et coefficient de Poisson
Anisotropie
Couplage thermoélastique

2 Le problème aux limites d’équilibre thermoélastique linéarisé

Linéarisation des équations d’équilibre
Formulation du problème aux limites
Schémas de résolution
Existence et unicité
Stabilité

3 Bilan
Problème d’équilibre thermoélastique linéarisé

le problème aux limites posé consiste à chercher les champs (u , σ

∼
)
sur Ω0 tels que
1


 ε
∼
= (grad u + (grad u )T ), ∀x ∈ Ω0 compatibilité


 2
Div σ + ρ0 f = 0, ∀X ∈ Ω0 équilibre



 ∼

(P) σ
∼
=Λ∼
: (ε
∼
− (T − T0 )α
∼
), ∀X ∈ Ω0 loi de comportement

 ∼

ui = uid , ∀X ∈ ∂Ωiu



 conditions aux limites


 t = σ N = t d , ∀X ∈ ∂Ωi

i ij j i td

conditions aux limites : en chaque point du bord, il faut donner

3 informations (déplacement ou effort) selon trois directions
orthogonales

Le problème aux limites d’équilibre thermoélastique linéarisé 42/64

 Plan

1 Loi de comportement thermoélastique linéarisée

Cas des déformations infinitésimales
Cas des transformations infinitésimales
Solides thermoélastiques linéaires isotropes
Module de Young et coefficient de Poisson
Anisotropie
Couplage thermoélastique

2 Le problème aux limites d’équilibre thermoélastique linéarisé

Linéarisation des équations d’équilibre
Formulation du problème aux limites
Schémas de résolution
Existence et unicité
Stabilité

3 Bilan
 Schéma de résolution par la méthode des
déplacements
déplacement
u
2


T


∇u


∇u


ε

déformations
ε

Le problème aux limites d’équilibre thermoélastique linéarisé 44/64

 Schéma de résolution par la méthode des
déplacements
déplacement
u
2



T


∇u




∇u





ε

déformations contraintes
ε σ
σ
 λ trace ε 1 2µε
loi de comportement
Le problème aux limites d’équilibre thermoélastique linéarisé 45/64
 Schéma de résolution par la méthode des
déplacements
équilibre


déplacement div σ ρ f  0
u σ n  t
2


T


∇u

Navier


∇u


ε

déformations contraintes
ε σ
σ λ  trace ε 1 2µε
loi de comportement
Le problème aux limites d’équilibre thermoélastique linéarisé 46/64
Schéma de résolution par la méthode des contraintes

équilibre


div σ ρ f  0
σ n  t

contraintes
σ

Le problème aux limites d’équilibre thermoélastique linéarisé 47/64

Schéma de résolution par la méthode des contraintes

équilibre


div σ ρ f  0
σ n  t

ε  
1 ν
σ ν
trace σ 1
E E

déformations contraintes
ε σ
loi de comportement

Le problème aux limites d’équilibre thermoélastique linéarisé 48/64

Schéma de résolution par la méthode des contraintes

équilibre


déplacement div σ ρ f  0
u σ n  t

compatibilité

ε  
1 ν
σ ν
trace σ 1
E E

déformations contraintes
ε σ
loi de comportement

Le problème aux limites d’équilibre thermoélastique linéarisé 49/64

 Plan

1 Loi de comportement thermoélastique linéarisée

Cas des déformations infinitésimales
Cas des transformations infinitésimales
Solides thermoélastiques linéaires isotropes
Module de Young et coefficient de Poisson
Anisotropie
Couplage thermoélastique

2 Le problème aux limites d’équilibre thermoélastique linéarisé

Linéarisation des équations d’équilibre
Formulation du problème aux limites
Schémas de résolution
Existence et unicité
Stabilité

3 Bilan
 la question de l’unicité
• Théorème des travaux virtuels dans le cas statique (en
l’absence de discontinuité)
Le champ de contraintes σ ∼
sur un corps matériel Ω soumis
aux conditions f , t d et u d vérifient les équations d’équilibre
div σ
∼
+ ρf = 0 en tout point si et seulement si, on a
Z Z Z Z
? ? d ?
σ
∼
:ε
∼
dV = ρf .u dV + t .u dS+ (σ
∼
.n ).u d dS
Ω Ω ∂Ωtd ∂Ωu

pour tout champ u ? tel que u ? = u d sur ∂Ωu .

Le problème aux limites d’équilibre thermoélastique linéarisé 51/64

 la question de l’unicité
• Théorème des travaux virtuels dans le cas statique (en
l’absence de discontinuité)
Le champ de contraintes σ ∼
sur un corps matériel Ω soumis
aux conditions f , t d et u d vérifient les équations d’équilibre
div σ
∼
+ ρf = 0 en tout point si et seulement si, on a
Z Z Z Z
? ? d ?
σ
∼
:ε
∼
dV = ρf .u dV + t .u dS+ (σ
∼
.n ).u d dS
Ω Ω ∂Ωtd ∂Ωu

pour tout champ u ? tel que u ? = u d sur ∂Ωu .

• Soient maintenant (σ , ε ) et (σ
∼ A ∼A
, ε ) deux solutions sur Ω
∼ B ∼B
du problème aux limites
Z Z Z Z
d
σ
∼A
:ε
∼A
dV = ρf .u A dV + t .u A dS+ (σ∼A
.n ).u d dS
Ω Ω ∂Ωtd ∂Ωu

Le problème aux limites d’équilibre thermoélastique linéarisé 52/64

 la question de l’unicité
• Théorème des travaux virtuels dans le cas statique (en
l’absence de discontinuité)
Le champ de contraintes σ ∼
sur un corps matériel Ω soumis
aux conditions f , t d et u d vérifient les équations d’équilibre
div σ
∼
+ ρf = 0 en tout point si et seulement si, on a
Z Z Z Z
? ? d ?
σ
∼
:ε
∼
dV = ρf .u dV + t .u dS+ (σ
∼
.n ).u d dS
Ω Ω ∂Ωtd ∂Ωu

pour tout champ u ? tel que u ? = u d sur ∂Ωu .

• Soient maintenant (σ , ε ) et (σ
∼ A ∼A
, ε ) deux solutions sur Ω
∼ B ∼B
du problème aux limites
Z Z Z Z
d
σ
∼A
:ε
∼A
dV = ρf .u A dV + t .u A dS+ (σ∼A
.n ).u d dS
Ω Ω ∂Ωtd ∂Ωu
Z Z Z Z
d
σ
∼A
:ε
∼B
dV = ρf .u B dV + t .u B dS+ (σ
∼A
.n ).u d dS
Ω Ω ∂Ωtd ∂Ωu

Le problème aux limites d’équilibre thermoélastique linéarisé 53/64

 la question de l’unicité
Z Z Z Z
σ
∼A
:ε
∼A
dV = ρf .u A dV + t d .u A dS+ (σ
∼A
.n ).u d dS
Ω Ω ∂Ωtd ∂Ωu
Z Z Z Z
σ
∼A
:ε
∼B
dV = ρf .u B dV + t d .u B dS+ (σ
∼A
.n ).u d dS
Ω Ω ∂Ωtd ∂Ωu

Le problème aux limites d’équilibre thermoélastique linéarisé 54/64

 la question de l’unicité
Z Z Z Z
σ
∼A
:ε
∼A
dV = ρf .u A dV + t d .u A dS+ (σ
∼A
.n ).u d dS
Ω Ω ∂Ωtd ∂Ωu
Z Z Z Z
σ
∼A
:ε
∼B
dV = ρf .u B dV + t d .u B dS+ (σ
∼A
.n ).u d dS
Ω Ω ∂Ωtd ∂Ωu
Z Z Z
σ
∼A
: ∆ε
∼
dV = ρf .∆u dV + t d .∆u dS
Ω Ω ∂Ωtd
∆u = u A − u B , ∆ε
∼
=ε
∼A
−ε
∼B

Le problème aux limites d’équilibre thermoélastique linéarisé 55/64

 la question de l’unicité
Z Z Z
σ
∼A
: ∆ε
∼
dV = ρf .∆u dV + t d .∆u dS
Ω Ω ∂Ωtd
Z Z Z
σ
∼B
: ∆ε
∼
dV = ρf .∆u dV + t d .∆u dS
Ω Ω ∂Ωtd

Le problème aux limites d’équilibre thermoélastique linéarisé 56/64

 la question de l’unicité
Z Z Z
σ
∼A
: ∆ε
∼
dV = ρf .∆u dV + t d .∆u dS
Ω Ω ∂Ωtd
Z Z Z
σ
∼B
: ∆ε
∼
dV = ρf .∆u dV + t d .∆u dS
Ω Ω ∂Ωtd
Z
∆σ
∼
: ∆ε
∼
dV = 0
Ω
∆σ
∼
=σ
∼A
−σ
∼B

Le problème aux limites d’équilibre thermoélastique linéarisé 57/64

 la question de l’unicité

Z
∆σ
∼
: ∆ε
∼
dV = 0
Ω
loi de comportement

σ
∼A
=Λ
∼
:ε
∼A
, σ
∼B
=Λ
∼
:ε
∼B
∼ ∼

d’où Z
∆ε
∼
:Λ
∼
: ∆ε
∼
dV = 0
Ω ∼

une condition suffisante d’unicité : Λ

∼
défini positif
∼

Le problème aux limites d’équilibre thermoélastique linéarisé 58/64

 Existence et unicité

• Problème aux limites de Dirichlet

u (X ) = u d (X ), ∀X ∈ ∂Ω0

l’existence et l’unicité sont assurées

le problème est bien posé
• Problème aux limites de Neumann

t =σ
∼
.n = t d , ∀X ∈ Ω0

l’existence est assurée ainsi que l’unicité à un mouvement de

corps rigide près (mais rotation infinitésimale)
• Conditions aux limites mixtes ∂Ωiu , ∂Ωitd
• Conditions aux limites non linéaires (contact)

Le problème aux limites d’équilibre thermoélastique linéarisé 59/64

 Théorème de superposition
Le problème (P) ainsi posé est linéaire. On entend par là que si
• (u A , σ
∼
A ) est une solution du problème (P) pour les données
A dA dA
(f , ui , ti );
• et si (u B , σ
∼
B ) est une solution du problème (P) pour les
B
données (f , uidB , tidB );
alors, pour tous les réels λA et λB ,
• (λA u A + λB u B , λA σ∼
A + λ σ B ) est une solution du problème
B∼
(P) pour les données
(λA f A + λB f B , λA uidA + λB uidB , λA tidA + λB tidB ).
Cette propriété est parfois appelée principe de superposition. Elle
permet de calculer l’état de contrainte et de déformation d’un solide pour
plusieurs chargement successifs. C’est tout à fait licite à condition
toutefois de veiller à chaque étape à la validité du cadre HPP
(transformation infinitésimale et états successifs quasi–naturels).
Pour l’appliquer, il faut que les domaines ∂Ωiu et ∂Ωid soient identiques
pour les deux chargements A et B.
Le problème aux limites d’équilibre thermoélastique linéarisé 60/64
 Plan

1 Loi de comportement thermoélastique linéarisée

Cas des déformations infinitésimales
Cas des transformations infinitésimales
Solides thermoélastiques linéaires isotropes
Module de Young et coefficient de Poisson
Anisotropie
Couplage thermoélastique

2 Le problème aux limites d’équilibre thermoélastique linéarisé

Linéarisation des équations d’équilibre
Formulation du problème aux limites
Schémas de résolution
Existence et unicité
Stabilité

3 Bilan
 la question de la stabilité

flambage
localisation
Warning! Le cadre HPP a un domaine de validité limité. Pour
déterminer ce domaine, il faut en sortir!
Lorsque l’unicité de la solution est perdue, certaines solutions
peuvent être instables. Exemple : état de compression homogène
dans une barre élancée (flambage d’Euler)

Le problème aux limites d’équilibre thermoélastique linéarisé 62/64

 Plan

1 Loi de comportement thermoélastique linéarisée

Cas des déformations infinitésimales
Cas des transformations infinitésimales
Solides thermoélastiques linéaires isotropes
Module de Young et coefficient de Poisson
Anisotropie
Couplage thermoélastique

2 Le problème aux limites d’équilibre thermoélastique linéarisé

Linéarisation des équations d’équilibre
Formulation du problème aux limites
Schémas de résolution
Existence et unicité
Stabilité

3 Bilan
 Bilan
équilibre


déplacement div σ ρ f  0
u σ n  t
2


T

compatibilité

∇u

Navier


∇u

ε  
1 ν
σ ν
 trace σ  1
 E E
ε

déformations contraintes
ε σ
σ λ  trace ε 1 2µε
loi de comportement

Bilan 64/64