Ecole Nationale des Sciences Appliquées

OUJDA Année universitaire 2019-2020

Fiche des travaux dirigés des

Méthodes numériques calcul de structures en éléments finis avec solution

Chapitre 3- Méthode Elément finis 2D domaine élastique

Exercice1
Une membrane en traction simple, de dimension 1x2 et d’épaisseur l’unité charge par une force de
contacte d’intensité l’unité. On néglige son poids et la membrane travaille en contraintes planes. Le
matériau est élastique de module de Young l’unité et le coefficient de Poisson 0.4. Le domaine d’étude
est présenté dans la figure ci-dessous. Le maillage du barrage en éléments finis triangulaire est donné
sur la même figure ci jointe. On suppose que la membrane est encastrée au nœud 1 et en appui
simple en nœud 2 et 3. Le chargement extérieurs est appliqué aux nœud 3 et 4 voir figure ci-dessous.

1) Etablir la matrice de rigidité pour chaque élément de la membrane ?

2) Etablir la matrice assemblée globale de rigidité ?
3) Etablir le vecteur de déplacement ?
4) Déterminer les contraintes et la déformée dans chaque élément ?
5) Calculer les efforts de liaisons ?

1
Tableau des coordonnées :
N° de nœud X Y
1 0 0
2 1 0
3 0 2
4 1 2

Tableau de connectivité :
N° d’élément Connectivité
1 1 2 3
2 4 3 2

Tableau des conditions aux limites cinématiques :

N° de nœud U V
1 0 0
2 ? 0
3 0 ?
4 ? ?

Tableaux des conditions aux limites statiques :

N° de nœud Fx Fy
1 ? ?
2 0 ?
3 ? 0.5
4 0 0.5

2
Solution
1) Etablir la matrice de rigidité pour chaque élément de la membrane ?
Élément 1
𝑥1 X1 = 0
𝑦1 Y1 = 0
𝑥2 X2 = 1
𝑥, 𝑦 = =
𝑦2 Y2 = 0
𝑥3 X3 = 0
y3 Y3 = 2

𝑎, 𝑏, 𝑐 = 𝑏1 = 𝑦2 − 𝑦3 = −2 𝑏2 = 𝑦3 − 𝑦1 = 2 𝑏3 = 𝑦1 − 𝑦2 = 0
𝑐1 = 𝑥3 − 𝑥2 = −1 𝑐2 = 𝑥1 − 𝑥3 = 0 𝑐3 = 𝑥2 − 𝑥1 = 1
La section de l’élément 1 (triangle) : S1 = 1 m2

Avec: ϑ = 0,4

1 ϑ 0
E ϑ 1 0 1 1 0,4 0
C =
1−ϑ 2 1 − ϑ = 0,8 0,4 1 0
0 0 0 0 0,3
2

1 b1 0 b2 0 b3 0 1 −2 0 2 0 0 0
B1 = 0 c1 0 c2 0 c3 = 0 −1 0 0 0 1
2 2
c1 b1 c2 b2 c3 b3 −1 −2 0 2 1 0

La matrice de rigidité dans le repère global s’ écrit sous la formule suivante :

1 2 3 4 5 6
1 4.3 1.4 −4 −0.6 −0.3 −0.8
1 2 1.4 2,2 −0,8 −1.2 −0.6 −1.0
𝐾1 = 3 −4 −0.8 4 0 0 0.8
3.36
4 −0.6 −1.2 0 1,2 0.6 0
5 −0.3 −0.6 0 0.6 0.3 0
6 −0.8 −1 0.8 0 0 1
Elément 2
De même que l’élément 1

1 2 3 4 5 6
1 4.3 1.4 −4 −0.6 −0.3 −0.8
1 2 1.4 2,2 −0,8 −1.2 −0.6 −1.0
𝐾2 = 3 −4 −0.8 4 0 0 0.8
3.36
4 −0.6 −1.2 0 1,2 0.6 0
5 −0.3 −0.6 0 0.6 0.3 0
6 −0.8 −1 0.8 0 0 1
3
 2) Etablir la matrice assemblée globale de rigidité ?

L’assemblage des matrices globales se fait en utilisant l’expression suivante :

Kass = k1 + k2

1 2 3 4 5 6
1 4.3 1.4 −4 −0.6 −0.3 −0.8
1
2 1.4 2,2 −0,8 −1.2 −0.6 −1.0
Soit 𝐾1 = 3.36 3 −4 −0.8 4 0 0 0.8
4 −0.6 −1.2 0 1,2 0.6 0
5 −0.3 −0.6 0 0.6 0.3 0
6 −0.8 −1 0.8 0 0 1
Et

7 8 5 6 3 4
7 4.3 1.4 −4 −0.6 −0.3 −0.8
1 8 1.4 2,2 −0,8 −1.2 −0.6 −1.0
𝐾2 = 5 −4 −0.8 4 0 0 0.8
3.36
6 −0.6 −1.2 0 1,2 0.6 0
3 −0.3 −0.6 0 0.6 0.3 0
4 −0.8 −1 0.8 0 0 1

1
𝐾𝑎𝑠𝑠𝑏 = 3.36 x

4.3 1.4 -4 -0.6 -0.3 -08 0 0

1.4 2.2 -0.8 -1.2 -0.6 -1 0 0
-4 --0.8 4.3 0 0 1.4 -0.3 -0.6
-0.6 -1.2 0 2.2 1.4 0 -0.8 -1
-0.3 -0,6 0 1.4 4.3 0 -0.4 -0.8
-0.8 -1 1.4 0 0 2.2 -0,6 -1.2
0 0 -0.3 -0.8 -0.4 -0.6 4.3 1.4
0 0 -0.6 -1 -0.8 -1.2 1.4 2.2

4
 3) Etablir le vecteur de déplacement ?

Pour déterminer le vecteur déplacement, il faut éliminer la ligne et la colonne correspondantes aux
déplacements nulles.

U1 = 0
V1 = 0
U2 = ?
V 2= 0
U3 = 0
V3 = ?
U4 = ?
V4 = ?

1
𝐾𝑎𝑠𝑠𝑏 = 3.36 𝑥

4.3 1.4 -4 -0.6 -0.3 -08 0 0

1.4 2.2 -0.8 -1.2 -0.6 -1 0 0
-4 --0.8 4.3 0 0 1.4 -0.3 -0.6
-0.6 -1.2 0 2.2 1.4 0 -0.8 -1
-0.3 -0,6 0 1.4 4.3 0 -0.4 -0.8
-0.8 -1 1.4 0 0 2.2 -0,6 -1.2
0 0 -0.3 -0.8 -0.4 -0.6 4.3 1.4
0 0 -0.6 -1 -0.8 -1.2 1.4 2.2

Finalement il reste :
u2 0
4.3 1.4 −0.3 −0.6
1.4 2,2 −0,6 −1.2 v3
1 0.5
−0.3 −0.6 4.3 1.4 =
3.36
−0.6 −1.2 1.4 2.2 u4 0
v4 0.5

5
On trouve :
u2 = −0.4

v3 = 2

u4 = −0.4
v4 = 2

4) Déterminer les contraintes et la déformée de chaque élément ?

εxx
ε = εyy = B × U = L × N × U
2εxy

Elément 1

0
εxx −2 0 2 0 0 0 0 −0.4
1 −0.4
ε = εyy = × 0 −1 0 0 0 1 × = 1
2εxy 2 0
−1 −2 0 2 1 0 0 0
2

σxx 1 0.4 0 −0.4 0

σ
𝜎 = yy = 0.8 × 0.4 1 0 × 1 = 1
σxy 0 0 0.3 0 0
Elément 2

0
εxx −2 0 2 0 0 0 0 −0.4
1 −0.4
ε = εyy = × 0 −1 0 0 0 1 × = 1
2εxy 2 0
−1 −2 0 2 1 0 0 0
2

6
 σxx 1 0.4 0 −0.4 0
σ
𝜎 = yy = 0.8 × 0.4 1 0 × 1 = 1
σxy 0 0 0.3 0 0

La solution est homogène en déformation et contrainte, il s’agit de la solution exacte du problème de

traction simple d’une plaque rectangulaire.

5) Calculer les efforts de liaisons ?

Pour déterminer les forces de liaisons on reprend la matrice de rigidité assemblée et nous remplaçons
les déplacements calculés, après calcul on trouve :
Fx1= 0.0 Fy1= -0.5 Fy2= -0.5 Fx3=0.0
On peut vérifier l’exactitude de cette solution en appliquant le PFS.