TP : MÉTHODE DES ÉLÉMENTS FINIS
Adrien Socié
Dans ce TP dédié à la méthode des éléments finis, nous nous intéresserons à la simulation du com-
portement d’une poutre soumise à son poids propre avec et sans traction. L’objectif est de voir la mise
en oeuvre numérique des notions abordées en cours ainsi que de faire le lien avec les différents cours
utilisant les éléments finis que vous suivez lors de votre 4ème année. Pour ce faire, nous utiliserons un
script python afin d’effectuer les simulations.
Vous devrez rendre sur moodle pour le 11 juin au plus tard un compte rendu de tp au format pdf
ainsi que le ou les fichiers pythons correspondant à votre travail.
Système Modèle Poutre Système éléments finis
Figure 1 – Schématisation de l’étude. La poutre de longueur L = 1 m ; de diamètre 0.1 m, de densité
900 kg.m−3 et de module de Young 0.1 GPa.
1 Cas d’une poutre soumise à son poids propre
Nous nous intéressons dans un premier cas à une poutre soumise à son poids propre encastré dans
le sol (cf Figure 1). La poutre a un comportement élastique et isotrope. Nous nous intéressons à la
résolution statique du système
1) Écrire la formulation forte du problème
2) Préciser la différence entre la formulation forte, faible (ou variationnelle) et éléments finis
3) Qu’est qu’une condition de type Dirichlet ?
4) Qu’est qu’une condition de type Neumann et Neumann homogène ? Par défaut quelle condition
est utilisée dans les logiciels éléments finis ?
5) Donner la solution analytique du système (déplacement et contrainte).
6) Qu’est ce qu’une fonction de forme ?
7) Qu’est ce que l’espace isoparamétrique ?
8) Donner la formulation faible du système
9) Donner la formulation élément finis. Spécifier la matrice de rigidité élémentaire ainsi que le vecteur
de force pour un élément.
10) Dans les codes éléments finis qu’est ce qui est calculé au niveau des points de gauss ? Qu’est ce
qui est calculé aux nœuds ?
11) Étudier la réponse du système pour 5 éléments. De même pour 50 éléments. Qu’observez vous ?
1
12) Question bonus : Mettre en place un algorithme permettant d’étudier la convergence aux maillages
c’est-à-dire la réduction de l’erreur entre la solution analytique et numérique en fonction du nombre de
maille. Présenter vos résultats. Vous êtes libre d’utiliser la formule qui vous semble la plus pertinente
pour calculer l’erreur, c’est-à-dire que vous pouvez utiliser une erreur relative ainsi que choisir la norme.
Vous devez tout de même expliciter votre calcul.
2 Cas d’une poutre soumise à son poids propre et en traction
Maintenant la barre est soumise à un effort de traction de valeur Fimp = 50N à l’extrémité de la
barre.
13) Donner la nouvelle formulation forte
14) Donner la nouvelle formule analytique donnant le déplacement et la contrainte
15) Pourquoi ne pouvons nous pas résoudre en quasi-statique un système ayant uniquement des
conditions de Neumann ?
16) Créer un nouveau fichier python et modifier l’algorithme
17) Étudier les résultats
18) Est-ce que la dépendance aux maillages est similaire ?
3 Cas d’une poutre hétérogène
Nous gardons les mêmes conditions limites que la section précédente. Dans cette section, nous nous
intéressons à la réponse mécanique d’une poutre hétérogène composée du matériau 1 et 2 ayant respecti-
vement un module de Young de E1 = 0.1 GPa et E2 = 1. GPa. Nous considérons que les deux matériaux
ont la même densité massique.
Dans un premier temps, la poutre est composée du matériau 1 de sa base (0) à la longueur L1, puis
elle est constituée du matériau 2 de la longueur L1 à longueur L où est appliquée la condition limite de
traction.
19) Créer un nouveau fichier python et modifier l’algorithme. Nous considérons que L1 = 0.5L et nous
considérons un maillage composé de 10 éléments. Expliciter dans votre compte rendu vos modifications
et analysez les résultats.
20) La densité du matériau 2 vaut maintenant 1000 kg.m−3 . Modifier le python, expliquer les modifica-
tions et étudier les résultats. Dans le cas où votre discrétisation est de 5 éléments, comment pourriez-vous
résoudre le système ?
Figure 2 – Schématisation des modèles de Voigt et Reuss (issue du cours d’homogénéisation de Farid
Benboujema en 2016).
Maintenant, nous nous intéressons à une poutre composée d’un matériau composite constitué des
matériaux 1 et 2. Afin d’estimer son module de Young nous utiliserons les modèles d’homogénéisation
de Reuss et Voigt présenté dans la Figure 2. La quantité du matériau 2 est représenté par sa fraction
volumique (Volume du matériau 2/Volume du composite), notée ϕ.
21) Donner le module de Young homogénéisé obtenu par le modèle de Reuss et Voigt.
22) Modifier le code. Étudier les deux modèles pour une fraction volumique de 0., 0.3 et 0.6.