UNIVERSITÉ HASSAN II DE CASABLANCA

FACULTÉ DES SCIENCES BEN M’SIK

DÉPARTEMENT DE PHYSIQUE

Elasticité non linéaire 3D

Master Mécanique-Energétique

Réalisé par :
EL FINTI ISMAIL
EL GARRAB ASSIA Encadré: [Link]
 PLAN
01 02
Introduction Elasticité non
linéaire 3D

03 04
Comparaison Entre CONCLUSION
Calcul Linéaire Et Non
Linéaire
 01
INTRODUCTION
Qu’est-ce que l’´élasticité NON LINEAIRE ?
L'élasticité non linéaire décrit le comportement réversible des matériaux sous de
grandes déformations, lorsque la relation contrainte-déplacement n'est plus
proportionnelle. Cela englobe des effets géométriques et matériels complexes qui
ne peuvent être modélisés par l'élasticité linéaire.

Il existe trois types de non-linéarité que l'on peut observer dans le comportement des structures:
 Non-linéarités Géométriques

se manifestent lorsque les déplacements deviennent

importants, entraînant une déformation visible de la
structure et une variation de sa rigidité

Non-linéarités Matérielles
se produisent lorsque les déformations ne sont plus
proportionnelles aux contraintes appliquées. Ces déformations
peuvent devenir irréversibles, et les propriétés des matériaux
peuvent évoluer avec le temps.
non-linéarités des conditions limites

surviennent lorsque la structure est contrainte dans un seul

sens, ce qui soulève des problèmes de contact, parfois
accompagnés de frottement.
Quelles sont les applications de l'élasticité non linéaire?

Aéronautique Automobile Génie civil Biomécanique

 02
Elasticité non linéaire
3D
 DESCRIPTION DU MOUVEMENT ET DE LA DEFORMATION

A l'instant t=0 A l'instant t

Application d’un
chargement défini sur Ω𝟎

M N
n
m
3

Configuration de Lagrange 2 Configuration d’Euler

0
1

✓ La description lagrangienne est défini par : 𝐗 𝒙(𝑿,t)

✓ Le vecteur déplacement est défini par : 𝑼 (𝑿)=𝒙(𝑿,t) − 𝑿

 le tenseur gradient de la transformation :

On définit le tenseur gradient de la transformation par: ധ

𝐝𝐱 = 𝐅𝐝𝐗
𝝏𝒙𝒊
On a: 𝒙 = 𝒙 𝐗, 𝑡 Soit: 𝒅𝒙𝒊 = 𝒅𝑿𝑱
𝝏𝑿𝑱
Alors:
𝑭𝒊𝑱 =
𝝏𝒙𝒊
𝝏𝑿 𝑱
→C’est le tenseur gradient de la transformation
On relie le tenseur gradient de la transformation au gradient de déplacement

ന𝐅 = 𝜕𝒙 = 𝟏
ന + 𝛻𝐮
ധ ou 𝑭𝒊𝒋 =
𝝏𝒙𝒊
= 𝜹𝒊𝒋 + 𝑼𝒊,𝒋
𝜕𝐗 𝝏𝑿𝑱
 Le tenseur des déformations de Green-Lagrange

Au début, nous allons établir la définition de l'allongement:

𝒕∥𝑴𝑵∥(𝟏(𝑪−𝟏) )∥𝑴𝑵∥
𝟏 ∥𝒎𝒏∥𝟐 −∥𝑴𝑵∥𝟐 𝟐
𝜶= =
𝟐 ∥𝑴𝑵∥𝟐 𝒕∥𝑴𝑵∥.∥𝑴𝑵∥

ധ
On observe l'apparition d'un tenseur d'ordre deux, appelé tenseur des dilatations :𝐂 = 𝐭𝐅ധ𝐅ധ

Ce qui définit le tenseur des déformations de Green-Lagrange

𝟏
𝛄ധ = 𝟐(𝐂ധ − 𝟏
ന)
Expression en fonction du gradient du déplacement:

𝟏 𝟏
𝛄ധ = ( 𝒕𝑭𝑭 ന =
ന − 𝟏) ന + 𝑡𝛁
{(𝟏 ന+ 𝛁
ന 𝒖). (𝟏 ന = 𝟏 (𝛁
ന 𝒖)-𝟏} ന 𝒖) + 𝟏( 𝑡𝛁
ന 𝒖+ 𝑡 𝛁 ന 𝒖. 𝛁
ന 𝒖)
𝟐 𝟐 𝟐 𝟐

Terme linéaire Terme non-linéaire

 Energie de déformation

➢ L'énergie de déformation représente la quantité de travail nécessaire pour déformer un solide.

L’ énergie élastique stockée définit par : 𝝆𝟎 𝝎(𝜸)
La masse volumique La densité de l’énergie

L énergie de déformation : 𝐖 𝐮 = ‫𝛄 𝛚 𝟎𝛒 𝛀׬‬ധ 𝐮 . 𝐝𝛀𝟎

𝟎

ധ ∶ 𝐋ӖӖ ∶ 𝛄ധ ). 𝐝𝛀𝟎
𝟏
Pour un matériau élastique ‫𝛄 𝛚𝟎𝛒 𝛀׬‬ധ 𝐮 . 𝐝𝛀𝟎 = ‫𝛄 𝟐( 𝛀׬‬
𝟎 𝟎

D’où l énergie de déformation d un solide s écrit :

ന ∶ 𝑳ധധ ∶ 𝜸
𝟏
𝑾 𝒖 = ‫𝜸 ( 𝜴׬‬ ന). 𝒅𝜴𝟎
𝟎 𝟐

ധധ :Le tenseur des propriétés élastiques

Avec 𝑳
 Energie potentielle

𝛛𝐮𝛀𝟎
𝐅Ԧ 𝛛𝑭𝛀𝟎
Le travail développé par ces efforts pendant un déplacement u est:
Ԧ Ω0
𝑓,
Ԧ 𝐮 . 𝐝𝛀𝟎 + ‫׬‬
𝐏𝐝 𝐮 = ‫𝐟 𝛀׬‬. Ԧ 𝐮 . 𝐝𝐒𝟎
𝐅.
𝟎 𝛛 𝐅𝛀𝟎

L'Energie potentielle est donnée par :

𝐏 𝐮 = 𝐖 𝐮 − 𝐏𝐝 𝐮 Sur 𝛺0 en définit les forces de

volume 𝑓 par unités de volume
donc: non déformée et de surface 𝐹
par unité de surface non
déformée.
𝐏 𝐮 = ‫𝛄 ( 𝛀׬‬ധ ∶ 𝐋ӖӖ ∶ 𝛄ധ ). 𝐝𝛀𝟎- ‫𝐟 𝛀׬‬.
𝟏 Ԧ 𝐮 . 𝐝𝛀𝟎 − ‫׬‬ Ԧ 𝐮 . 𝐝𝐒𝟎
𝐅.
𝟎 𝟐 𝟎 𝛛 𝐅 𝛀𝟎
 Équations d’équilibre
Théorème:
Les positions d'équilibre correspondent aux points stationnaires de l'énergie potentielle sur
l'ensemble des champs de déplacements cinématiquement admissibles.

𝜹𝑷 𝒖 = 𝟎 ∀𝜹𝒖 C.A à 0

𝛿𝑷 𝒖 = ‫𝐋( 𝛀׬‬ӖӖ ∶ 𝛄ധ ) ∶ 𝛅ധ𝛄𝐝𝛀𝟎 - ‫𝐟 𝛀׬‬.

Ԧ 𝜹𝐮 . 𝐝𝛀𝟎 − ‫׬‬
𝛛
Ԧ 𝛅𝐮 . 𝐝𝐒𝟎
𝐅.
𝟎 𝟎 𝐅𝛀𝟎

On définit le second tenseur des contraintes de Piola-Kirchhoff S comme un nouveau tenseur.

𝐒ധ = 𝐋ӖӖ ∶ 𝛄ധ 𝐮 ou 𝑺𝒊𝒋 = 𝑳𝒊𝒋𝒌𝒍 𝜸𝒌𝒍

Ainsi, la variation de l'énergie potentielle peut s'exprimer sous la forme suivante:

Ԧ 𝛅𝐮 . 𝐝𝛀𝟎 + ‫׬‬
‫𝐒 𝛀׬‬ധ ∶ 𝛅ധ𝛄𝐝𝛀𝟎 = ‫𝐟 𝛀׬‬. Ԧ 𝛅𝐮 . 𝐝𝐒𝟎
𝐅. ∀𝜹𝒖. 𝑪. 𝑨 à 𝟎
𝟎 𝟎 𝛛 𝐅 𝛀𝟎
 Equations d’équilibre locales

𝟏 𝟏 𝐭
on a : 𝛅 𝛄ധ 𝐮 ന 𝛅𝐮+ 𝐭𝛁
= (𝛁 ന 𝛅𝐮)+ (𝛁ന 𝛅𝐮. ന 𝐮+ 𝐭𝛁
𝛁 ന 𝐮. 𝛁
ന 𝛅𝐮)
𝟐 𝟐

𝟏 𝟏
Soit : 𝑺𝒊𝒋 𝜹𝜸𝒊𝒋 = 𝑺𝒊𝒋 (𝛅𝒖𝒊,𝒋 + 𝜹𝒖𝒋,𝒊 )+ 𝑺𝒊𝒋 (𝜹𝒖𝒌,𝒊 𝒖𝒌,𝒋 + 𝒖𝒌,𝒊 𝜹𝒖𝒌,𝒋 )
𝟐 𝟐

𝟏
𝐒𝐢𝐣 𝛅𝛄𝐢𝐣 = 𝐒𝐢𝐣 𝛅𝐮𝐢,𝐣 + 𝛅𝐮𝐣,𝐢 +𝐒𝐤𝐣𝛅𝐮𝐢,𝐣 𝐮𝐢,𝐤 𝐒𝐢𝐣 𝛅𝛄𝐢𝐣 = (𝐒𝐢𝐣 + 𝐒𝐤𝐣𝐮𝐢,𝐤 )𝛅𝐮𝐢,𝐣
𝟐

En introduisant le 1é𝑚𝑒 tenseur des contraintes de Piola kirchhoff 𝑇ധ comme nouveau tenseur:

𝐓𝐢𝐣 = 𝐒𝐢𝐣 + 𝐒𝐤𝐣𝐮𝐢,𝐤

On fait intégration par partie , on trouve :

‫𝐒 𝛀׬‬ധ ∶ 𝛅ധ𝛄𝐝𝛀𝟎 = ‫ 𝒋𝒊𝑻 𝛀׬‬. 𝛅𝒖𝒊,𝒋 . 𝐝𝛀𝟎 = − ‫𝒋𝒊𝑻 𝛀׬‬,𝒋 . 𝛅𝒖𝒊 . 𝐝𝛀𝟎 + ‫𝛛׬‬ 𝑻𝒊𝒋 . 𝑵𝒋 . 𝛅𝒖𝒊 . 𝐝𝑺𝟎
𝟎 𝟎 𝟎 𝐅 𝛀𝟎

− න 𝐓𝐢𝐣,𝐣 . 𝛅𝐮𝐢 . 𝐝𝛀𝟎 + න 𝐓𝐢𝐣. 𝐍𝐣 . 𝛅𝐮𝐢 . 𝐝𝐒𝟎 = න 𝐟𝐢 . 𝛅𝐮𝐢 . 𝐝𝛀𝟎 + න 𝐅𝐢 . 𝛅𝐮𝐢 . 𝐝𝐒𝟎
𝛀𝟎 𝛛𝐅 𝛀𝟎 𝛀𝟎 𝛛𝐅 𝛀𝟎

‫𝐣𝐢𝐓( 𝛀׬‬,𝐣 +𝐟𝐢 ). 𝛅𝐮𝐢 . 𝐝𝛀𝟎 + ‫𝛛׬‬ (−𝐓𝐢𝐣. 𝐍𝐣 + 𝐅𝐢 ). 𝛅𝐮𝐢 . 𝐝𝐒𝟎 = 𝟎 ⇒ ∀𝛅𝐮𝐢
𝟎 𝐅 𝛀𝟎

D’après l’équilibre on trouve :

𝛛𝐓𝐢𝐣
+ 𝐟𝐢 = 𝟎 ⟹ 𝒅𝒂𝒏𝒔 𝛀𝟎 𝐝𝐢𝐯𝐓 ന + 𝐟Ԧ = 𝟎 ⟹ 𝐝𝐚𝐧𝐬𝛀𝟎
ቐ 𝛛𝐗𝐢 OU ቐ
𝐓𝐢𝐣𝐍𝐣 = 𝐅𝐢 ⟹ 𝐬𝐮𝐫 𝛛𝐅𝛀𝟎 𝐓ന . 𝐍 = 𝐅Ԧ ⟹ 𝐬𝐮𝐫 𝛛𝐅𝛀𝟎
On peut articuler le problème d'élasticité comme suit :
𝛛𝐮𝛀𝟎

Ԧ Ω0
𝑓,
𝑈 𝑒𝑠𝑡 𝐶. 𝐴 ⟺ 𝑈 = 𝑈𝑑 𝑠𝑢𝑟 𝜕𝑢 𝛺0
൞ 𝑆Ӗ 𝑒𝑠𝑡 𝑆. 𝐴
𝑆Ӗ = 𝐿ധധ : 𝛾Ӗ
𝛛𝑭𝛀𝟎
 03
Comparaison Entre
Calcul Linéaire Et
Non Linéaire
 DEFFIRENCE

Elasticité linéaire Elasticité non linéaire

Hypothèse des petites perturbations Toutes les perturbations sont acceptées

Méthodes numériques avancées nécessitant

Méthodes analytiques et itératives simples.
des itérations multiples.

Structures complexes (poutres plastiques,

Structures simples (poutres, dalles, etc.)
coques, etc.).

Relation de Hooke est proportionnel à la Comportement matériel plus complexe,

déformation incluant plasticité, viscoélasticité, etc.
 Tenseur de Green-Lagrange

𝛄ധ
linéaire

𝟏 𝒕
Non linéaire 𝛄ധ = ( 𝑭𝑭 ന)
ന− 𝟏
𝟐
𝟏 𝒕
𝛄ധ = ( 𝑭𝑭 ന)
ന− 𝟏 𝟏
𝟐 ന+
{(𝟏 𝑡
𝛁 ന+ 𝛁
ന𝒖). (𝟏 ന} = 𝟏 (𝛁
ന𝒖)-𝟏 ന𝒖+ 𝑡𝛁
ന𝒖) + 𝟏 𝑡
(𝛁 ന𝒖. ന𝒖)
𝛁
𝟐 𝟐 𝟐
𝟏
ന + 𝑡𝛁
= 𝟐 {(𝟏 ന+ 𝛁
ന𝒖). (𝟏 ന} = 𝟏 (𝛁
ന𝒖)-𝟏 ന𝒖) + 𝟏( 𝑡𝛁
ന𝒖+ 𝑡𝛁 ന𝒖. 𝛁
ന𝒖)
𝟐 𝟐 ➢ On va négligé le terme encadré en rouge qui
représente la non linéarité du tenseur
✓ Cette formule adaptée à n’importe quelle problème
𝟏
ന𝒖+ 𝑡𝛁
ധ = (𝛁
𝛄 ന𝒖)
𝟐 Tenseur de Green Lagrange
linéarisé
✓ Cette formule est valable juste pour l’HPP
 Energie de déformation

Non linéaire 𝑾 = න 𝝆𝟎𝒘𝒅𝜴𝟎 linéaire

𝜴𝟎

𝟏 ധധ ∶ 𝜸
ന∶𝑳
𝑾 𝒖 = ‫𝜸 𝟐( 𝜴׬‬ ന). 𝒅𝜴𝟎
𝟎
𝟏 ധധ ∶ 𝜺ധ ). 𝒅𝜴
𝑾 𝒖 = ‫𝜺 𝟐( 𝜴׬‬ധ ∶ 𝑳
• 𝐒ധ = 𝐋ӖӖ ∶ 𝛄ധ 𝐮
• ധധ ∶ 𝜺ധ
ന =𝑳
𝝈
• En notation indicielle 𝑺𝒊𝒋 = 𝑳𝒊𝒋𝒌𝒍𝜸𝒌𝒍
Sij : le 2nd tenseur des contraintes de Piola-Kirchhoff • En notation indicielle : 𝝈𝒊𝒋 = 𝑳𝒊𝒋𝒌𝒍𝜺𝒌𝒍
𝟏
ന∶𝜸
𝑾 𝒖 = 𝟐 ‫𝑺( 𝜴׬‬ ന). 𝒅𝜴𝟎 𝟏
𝟎 • ന ∶ 𝜺ധ 𝒅𝜴
𝑾 𝒖 = 𝟐 ‫𝝈𝜴׬‬
𝟏
𝑾𝒖 = ‫𝛀𝒅 𝜸 𝑺 ׬‬
𝟐 𝛀𝟎 𝒊𝒋 𝒊𝒋 𝟎 •
𝟏
𝑾 𝒖 = 𝟐 ‫𝜴𝒅 𝒋𝒊𝜺 𝒋𝒊𝝈𝜴׬‬
Avec 𝑳ധധ :Le tenseur des propriétés élastiques
 Energie potentielle
𝐏 𝐮 = 𝐖 𝐮 − 𝐏𝐝 𝐮
Non linéaire linéaire

• Le travail développé par les efforts de • Le travail développé par les efforts de
déplacement u : déplacement u :
Ԧ 𝐮 . 𝐝𝛀𝟎 + ‫׬‬
𝐏𝐝 𝐮 = ‫𝐟 𝛀׬‬. Ԧ 𝐮 . 𝐝𝐒𝟎
𝐅. Ԧ 𝐮 . 𝐝𝛀𝟎+ ‫׬‬
𝐏𝐝 𝐮 = ‫𝐟 𝛀׬‬. Ԧ 𝐮 . 𝐝𝐒𝟎
𝐅.
𝟎 𝛛 𝐅 𝛀𝟎 𝟎 𝛛 𝐅 𝛀𝟎

• L’énergie de déformation : • L’énergie de déformation :

𝟏 ധധ ∶ 𝜸
ന∶𝑳
𝑾 𝒖 = ‫𝜸 𝟐( 𝜴׬‬ ന)𝒅𝜴𝟎
𝟏
ധധ ∶ 𝜺ധ )𝒅𝜴
𝑾 𝒖 = ‫𝜺 ( ׬‬ധ ∶ 𝑳 𝟎
𝟎 𝜴𝟎 𝟐

donc donc

P u = ‫׬‬Ω (2 γധ ∶ Lധധ ∶ γധ )dΩ0- ‫׬‬Ω Ԧf. u . dΩ0 − ‫𝜕׬‬ ധധ ∶ 𝜺ധ )𝒅𝜴 - ‫ ׬‬Ԧf. u . dΩ − ‫׬‬
1 𝟏
F. u . dS0 P u = ‫𝜺 𝟐( 𝜴׬‬ധ ∶ 𝑳 𝟎 Ω 0 F. u . dS0
0 0 F Ω0 𝟎 0 𝜕 F Ω0
 Les équations d’équilibre

Non linéaire linéaire

𝛛𝐓𝐢𝐣 𝛛𝝈𝐢𝐣
+ 𝐟𝐢 = 𝟎 ⟹ 𝒅𝒂𝒏𝒔 𝛀𝟎 + 𝐟𝐢 = 𝟎 ⟹ 𝒅𝒂𝒏𝒔 𝛀
ቐ 𝛛𝐗𝐢 ቐ 𝛛𝐗𝐢
𝐓𝐢𝐣𝐍𝐣 = 𝐅𝐢 ⟹ 𝐬𝐮𝐫 𝛛𝐅𝛀𝟎 𝝈𝐢𝐣𝐍𝐣 = 𝐅𝐢 ⟹ 𝐬𝐮𝐫 𝛛𝐅𝛀
Ou sous forme matricielle Ou sous forme matricielle
𝐝𝐢𝐯𝐓ന + 𝐟Ԧ = 𝟎 ⟹ 𝐝𝐚𝐧𝐬 𝛀𝟎 𝐝𝐢𝐯ന𝛔 + 𝐟Ԧ = 𝟎 ⟹ 𝐝𝐚𝐧𝐬 𝛀
ቐ • ൝
ന . 𝐍 = 𝐅Ԧ ⟹ 𝐬𝐮𝐫 𝛛𝐅𝛀𝟎
𝐓 ന. 𝐍 = 𝐅Ԧ ⟹ 𝐬𝐮𝐫 𝛛𝐅𝛀
𝛔
 04
CONCLUSION
 Prédire la réponse
du structures et du
matériaux utilisées

Afin d’analyser le
problème par un
sens physique
 MERCI POUR
VOTRE ATTENTION