Méthode des éléments finis
Thomas SAYET
Site Vinci-Lagrange 104
Modélisation unidimensionnelle
1
17/09/2018
Plan
Détermination des fonctions de base
Propriétés des fonctions de base
Détermination directe
Calcul de 𝐾𝑒 et 𝐹𝑒
Quadratique (3 nœuds)
Linéaire (2 nœuds) Cubique (4 nœuds)
2
17/09/2018
Détermination des fonction de base
Propriétés des fonctions de base
Sur chaque élément, on approche le champ de déplacement par :
𝒏
𝒖𝒆 𝒙 = 𝐍 𝐞 𝒖𝒆𝒏 = 𝑵𝒆𝒊 𝒙 {𝒖𝒊 }
𝒊=𝟏
avec
𝐍 𝐞 matrice des fonctions de base
𝒖𝒆𝒏 vecteur des valeurs (inconnues) du champ de déplacement aux
nœuds de l’élément.
𝑵𝒆𝒊 𝒙 fonctions (polynomiales) de base (ou fonctions d’interpolation) de
l’élément
𝒏 nombre de nœuds de l’élément
{𝒖𝒊 } valeur du champ de déplacement au nœud 𝑖
3
17/09/2018
Détermination des fonction de base
M x
Propriétés des fonctions de base
Elément barre à 3 nœuds 1𝐷 i j k
𝑢 𝑀 = 𝑁𝑖 𝑥 𝑢𝑖 + 𝑁𝑗 𝑥 𝑢𝑗 + 𝑁𝑘 𝑥 𝑢𝑘
𝒖𝒊
= 𝑵𝒊 𝑵𝒋 𝑵𝒌 𝒖𝒋 = 𝑵𝒆 {𝒖𝒆𝒏 }
𝒖𝒌
Les fonctions de base 𝑁𝑖𝑒 𝑥 sont telles que : 𝑁𝑖𝑒 𝑥𝑗 = 𝛿𝑖𝑗
Détermination des fonctions de base (pour chaque élément) : trouver un
polynôme tel qu’il prenne des valeurs connues en un certain nombre de
points (les nœuds).
Détermination possible si le nombre de termes du polynôme est égal au
nombre de nœuds.
4
17/09/2018
Détermination des fonction de base
Détermination directe 1 u1 2 u2 x
Cas 1D: barre à 2 nœuds
x1 x2
Approximation du déplacement sur l’élément : 𝑢 = 𝑢 𝑥 𝑥Ԧ
2 nœuds, 2 termes pour le polynôme (polynôme de degré 1) :
𝑢 𝑥 = 𝑎𝑥 + 𝑏
• Nœud 1:
𝑢 𝑥1 = 𝑢1 = 𝑎𝑥1 + 𝑏
• Nœud 2:
𝑢 𝑥2 = 𝑢2 = 𝑎𝑥2 + 𝑏
𝑢1 − 𝑢2
𝑎=
𝑥1 − 𝑥2
֞ 𝑢2 𝑥1 − 𝑢1 𝑥2
𝑏=
𝑥1 − 𝑥2
5
17/09/2018
Détermination des fonction de base
Détermination directe
On obtient:
𝑢1 − 𝑢2 𝑢2 𝑥1 − 𝑢1 𝑥2
𝑢 𝑥 = 𝑥+
𝑥1 − 𝑥2 𝑥1 − 𝑥2
𝑥2 − 𝑥 𝑥 − 𝑥1
𝑢 𝑥 = 𝑢1 + 𝑢
𝑥2 − 𝑥1 𝑥2 − 𝑥1 2
= 𝑁1 𝑥 𝑢1 + 𝑁2 𝑥 𝑢2
𝑁1 𝑥 et 𝑁2 𝑥 sont les fonctions de bases
Cas général, pour un élément unidimensionnel à 𝑛 nœuds, les fonctions de
base sont définies par :
𝒏
𝒙 − 𝒙𝒋
𝑵𝒊 𝒙 = ෑ
𝒙𝒊 − 𝒙𝒋
𝒋=𝟏
𝒋≠𝒊
6
17/09/2018
Détermination des fonction de base
Calcul de 𝑲𝒆 et 𝑭𝒆 1 u1 2 u2 x
Cas 1D: barre à 𝟐 nœuds
x1 x2
𝑢 𝑥 = 𝑁1 𝑥 𝑢1 + 𝑁2 𝑥 𝑢2
𝑥2 − 𝑥 𝑥2 − 𝑥 𝑥 − 𝑥1 𝑥 − 𝑥1
𝑁1 𝑥 = = 𝑒𝑡 𝑁2 𝑥 = =
𝑥2 − 𝑥1 𝐿 𝑥2 − 𝑥1 𝐿
Interpolation des déformations
𝑑𝑢 1 1 1 1 𝑢1
𝜀 𝑥 = = − 𝑢1 + 𝑢2 = −
𝑑𝑥 𝐿 𝐿 𝐿 𝐿 𝑢2
𝐵𝑒
Remarque : la déformation est constante sur cet élément
7
17/09/2018
Détermination des fonction de base
Calcul de 𝑲𝒆 et 𝑭𝒆 1 u1 2 u2 x
x1 x2
𝑲𝒆 = න 𝑩𝒆 𝑻
𝑯 𝑩𝒆 𝒅𝑽
𝛀𝐞
avec
𝒅𝑽 = 𝑺𝒅𝒙 , 𝑺 étant la section de la barre
𝝈 = 𝑬𝜺 loi de comportement 1D
֜ 𝑯 =𝑬
D’où
𝒙𝟐 𝒙𝟐
𝑬𝑺 −𝟏 𝑬𝑺 𝟏 −𝟏
𝑲𝒆 = 𝟐 න −𝟏 𝟏 𝒅𝒙 = 𝟐 න 𝒅𝒙
𝑳 𝒙𝟏 𝟏 𝑳 𝒙𝟏 −𝟏 𝟏
𝑬𝑺 𝟏 −𝟏 𝑬𝑺 𝟏 −𝟏
= 𝟐 𝒙𝟐 − 𝒙𝟏 =
𝑳 −𝟏 𝟏 𝑳 −𝟏 𝟏
Remarque : matrice indépendante des coordonnées des nœuds 8
17/09/2018
Détermination des fonction de base
Calcul de 𝑲𝒆 et 𝑭𝒆 1 u1 2 u2 x
Cas 1D: barre à 𝟐 nœuds
x1 x2
• 𝒇: charge volumique dans la direction 𝒙
• Déplacement:
𝑥2 − 𝑥 𝑥 − 𝑥1 𝑢1
𝑢 𝑥 = 𝑁1 𝑥 𝑢1 + 𝑁2 𝑥 𝑢2 = 𝑢2
𝐿 𝐿
𝑁𝑒
Second membre élémentaire:
𝒙𝟐
𝟏 𝟐
− 𝒙𝟐 − 𝒙
𝒇𝑺 𝒙𝟐
𝒙𝟐 − 𝒙 𝒇𝑺 𝟐 𝒙𝟏
𝑭𝒆 = න 𝑵𝒆 𝑻 𝒇𝒆 𝒅𝑽 = න 𝒙 − 𝒙 𝒅𝒙 = 𝒙𝟐
𝛀𝐞 𝑳 𝒙𝟏 𝟏 𝑳 𝟏 𝟐
𝒙 − 𝒙𝟏
𝟐 𝒙𝟏
𝒇𝑺 𝟐 𝒇𝑺𝑳 𝟏
𝒙𝟐 − 𝒙𝟏
= 𝟐
= 9
17/09/2018 𝟐𝑳 𝒙 −𝒙 𝟐 𝟏
