République du Bénin

*****
Ministère de l’Enseignement Supérieur et de la recherche
scientifique
*****
Ecole Polytechnique d’Abomey-Calavi (EPAC)
*****
Filière : Génie Civil
*****
Discipline : Mécanique des Eléments Finis

TP1 : TP1 : ELEMENT DE POUTRE DE

BERNOULLI

Membres du groupe : Chargé du cours :

1. ASSOGBA Expédit Dr. Rostand MOUTOU PITTI
2. CAKPO-BESSE Nestor Boris
3. SEHOUE Prince

ANNEE ACADEMIQUE : 2023-2024

TP1 : ELEMENT DE POUTRE DE BERNOULLI EPAC/GC 2023-2024

 TP1 : ELEMENT DE POUTRE DE BERNOULLI
Table des matières
Objectif :................................................................................................................................................. 2
Enoncé : .................................................................................................................................................. 2
A- Paramètres généraux : Opérateur Option .............................................................................. 3
B- Données du problème ................................................................................................................ 3
C- Géométrie : MAILLAGE.......................................................................................................... 4
D- Modèle de comportement du matériau.................................................................................... 5
E- Matériaux et caractéristiques des éléments............................................................................. 5
F- Matrice de rigidité : ....................................................................................................................... 6
G- Conditions aux limites ............................................................................................................... 6
H- Chargement................................................................................................................................ 6
I- Résolution du problème ................................................................................................................ 7
J- Solution de la théorie des poutres ................................................................................................ 7
K- Post-traitements : Calcul des déplacements ............................................................................ 8
L- Post-traitements : Calcul des efforts ........................................................................................ 9
MAILLAGE DE LA STRUCTURE .................................................................................................. 10
Epure des moments fléchissants ......................................................................................................... 11
Moment fléchissant sur la poutre discrétisée en 5 éléments ........................................................ 11
Moment fléchissant sur la poutre discrétisée en 10 éléments ...................................................... 12
Moment fléchissant sur la poutre discrétisée en 500 éléments .................................................... 12
Déformée de la structure poutre + treillis pour une discrétisation en 5, 10 et 500 éléments ........ 13
Effort tranchant sur la poutre discrétisée en 5 éléments ............................................................. 13
Effort tranchant sur la poutre discrétisée en 10 éléments ........................................................... 14
Effort tranchant sur la poutre discrétisée en 500 éléments ......................................................... 14
Commentaire : ..................................................................................................................................... 14

1
 TP1 : ELEMENT DE POUTRE DE BERNOULLI
Objectif :
Modélisation, étude comportementale d’une poutre hyperstatique
encastrée – articulée P1P2 (d’inertie Iz et de section S1) soumise à une
charge uniformément répartie d’intensité q soutenue par le treillis
formé des barres P2P3, P3P4 et P3P5 (de section S2) et simulation à
travers le logiciel Cast3M.

Enoncé :
Après avoir recopié (depuis l’ENT, Cours en ligne, 4AGC - MEF,
Documents, TP) et complété le fichier tp1 [Link]
permettant de résoudre ce problème, comparer les résultats obtenus
par ´éléments finis aux valeurs de la RdM (programmées dans le
fichier) en fonction de la finesse du maillage. Si, au lieu de la charge
uniforme, le chargement est constitué d’une force ponctuelle au centre
de la poutre P1P2, quel est le nombre d’´éléments conduisant à la
meilleure précision. Commentaires... Remarque : le fichier demande
seulement le nombre d’´éléments sur la poutre P1P

2
 TP1 : ELEMENT DE POUTRE DE BERNOULLI
A- Paramètres généraux : Opérateur Option

- Dimension et mode de calcul : la dimension de l’espace utilisé

est de 2 qui correspond à une étude du système suivant un plan à
savoir les axes x et y. Nous travaillons en contraintes plane et
évalueront les déformations correspondantes.

- "epsi 'LINEAIRE'" Cette ligne définit la précision de la solution.

Ici, la précision est définie comme linéaire..

Exemple :
*OPTI DIME 2 ELEM SEG2 ;
*OPTI EPSI LINEAIRE ;
Dans cet exemple, le calcul se fait en 2D. Le maillage génère des
segments à 2 nœuds ; éléments de type linéaire.

- Echo (0/1) indique le niveau de détail des messages : Echo 1

toutes les données apparaissent en sortie. Autrement dire cette
option active l'affichage des messages de diagnostic pendant
l'exécution du programme.
- Répétition de bloc (repeter bgene) :
Cette section marque le début d'une boucle qui semble répéter
l'analyse pour différentes configurations ou paramètres
B- Données du problème
Les dimensions et les propriétés matérielles de la poutre sont
spécifiées, y compris sa longueur (l), sa section transversale (a x a),
les moments d'inertie de la section transversale (Iz), le module de
Young (E), et la charge surfacique appliquée (p)
Exemple :
l = 1.; a = 0.01; s1 = a * a; iz = (a ** 4)/12.; s2 = 0.01 * s1;
E = 2.1e11; p = -1000. ;

3
 TP1 : ELEMENT DE POUTRE DE BERNOULLI

C- Géométrie : MAILLAGE
La géométrie représente la forme discrétisée du domaine à étudier.
Elle est composée d’objets de type MAILLAGE. La génération du
maillage sera faite à partir des géométries créées.
Il est utilisé essentiellement dans la résolution de ce TP en géométrie
les opérateurs points et droites (type d’objets lignes).
- Définition des coordonnées
P1 = 0. 0. ; P2 = l 0. ; P3 = l l ; P4 = 0. l ; P5 = 0. (2 * l)
- Création des lignes (droite)
La poutre P1P2 et les barres P2P3 ; P3P4 ; P3P5 sont désignées par les
éléments ligne L12, L23, L34, L35. L’ensemble des barres P2P3 ;
P3P4 et P3P5 forme un trélli lequel soutient la poutre P1P2.
Exemple :
▪ L23 = droi 1 P2 P3 ; L34 = droi 1 P3 P4 ; L35 = droi 1 P3 P5 ;
Treillis = l23 et l34 et l35.
▪ Poutre = droi nelem P1P2
La numérotation droi 1 : 1 spécifie le nombre d’élément soit par
définition L23 = droi 10 P1P2 est divisé en 10 éléments d’égale
longueur.
- Géométrie totale
geotot = treillis et poutre .
- L’option trace permet de générer une fenêtre graphique du
travail en cours d’exécution.
Exemple : trac QUAL geotot

4
 TP1 : ELEMENT DE POUTRE DE BERNOULLI

D- Modèle de comportement du matériau

Le modèle associe à un objet de type MAILLAGE un domaine
physique, une loi de comportement du matériau et la formulation
éléments finis.
Les noms des formulations et des modèles de comportement doivent
être écrits en toute lettre. Chaque modèle défini est affecté à une
géométrie (GEOi).
Exemple :
mo1 = mode treillis mecanique elastique isotrope barr;
mo2 = mode poutre mecanique elastique isotrope pout;
• L’objet mo1 définit une loi de comportement mécanique
élastique et affecte au maillage trellis des éléments de type barre.
• L’objet mo2 définit une loi de comportement mécanique
élastique et affecte au maillage poutre des éléments de type
poutre.

E- Matériaux et caractéristiques des éléments

On utilise l’opérateur MATE(RIAU) pour définir les propriétés
matérielles d’un modèle donné. Pour les propriétés géométriques, on
peut utiliser aussi l’opérateur MATE ou bien l’opérateur
CARA(CTERISTIQUE).
Exemple :
ma1 = mate mo1 youn E NU 0.3 sect s2
ma2 = mate mo2 youn E INRZ iz NU 0.3 sect s1
Dans ma1 et ma2, le module d’Young est décrit par une évolution
donnant la composante YOUN ; le coefficient de poisson est décrit
par NU et la section des éléments par sect.
Dans ma2 l’inertie de la poutre est décrite par INRZ.

5
 TP1 : ELEMENT DE POUTRE DE BERNOULLI

F- Matrice de rigidité :
L’opérateur RIGI construit la matrice de rigidité à partir du modèle et
des caractéristiques matérielles et géométriques
Exemple : k = (rigi mo1 ma1) et (rigi mo2 ma2)
G-Conditions aux limites
L’opérateur BLOQUER bloque un ou des degrés de liberté d’un objet
géométrique
cl1 = bloq Ux Uy Rz p1;
cl2 = bloq Ux Uy (p4 et p5);
cl = cl1 et cl2
• CL1 indique que le point P1 est encastré : les déplacements Ux
Uy et Rz sont bloqués.
• CL2 indique que seuls les déplacements UX et UY des points P4
et P5 sont bloqués

H-Chargement
Il existe des opérateurs spécifiques permettant de définir un
chargement (opérateurs FORCE, MOMENT, PRESSION...). Dans le
cadre de notre TP il s’agit de l’opérateur force surfacique : fsur
Exemple : char = fsur pout mo2 (0. P), indique que l’élément poutre
est chargé d’une force surfacique de composante P suivant la verticale.
Il s'agit d'une force surfacique (fsur) appliquée sur la poutre (pout) en
utilisant le modèle mécanique mo2 . La force surfacique est définie
comme ayant une composante en x nulle et une composante en y égale
à p.

6
 TP1 : ELEMENT DE POUTRE DE BERNOULLI

I- Résolution du problème
dep = resou (k et cl) char;
Cette ligne résout le système d'équations pour obtenir les déformations
de la structure. La fonction resou prend comme arguments la matrice
de rigidité `k`, les conditions aux limites `cl` et le chargement `char`.
J- Solution de la théorie des poutres
La section suivante est intitulée "Solution de la théorie des poutres" et
calcule les valeurs théoriques des déformations et des efforts dans la
poutre en fonction de la théorie des poutres. Voici ce que représente
chaque ligne fait :
-rac2 = 2. ** 0.5;: Calcule la racine carrée de 2.
- f = -3 * p * l /(8. * (1 + (6 * (1. + rac2) * Iz / (S2*l*l))));: Calcule
l'effort résultant (`f`) dans la poutre en utilisant la formule de la
théorie des poutres.
- u_rdm = ((f / 3.) + (p * l / 8.)) * ((l ** 3)/(e * iz));: Calcule le
déplacement (`u_rdm`) dans la poutre en utilisant la formule de la
théorie des poutres.
- N1_rdm = f; N2_rdm = -1. * f; N3_rdm = rac2 * f;: Calcule les
efforts normaux (N1_rdm, N2_rdm, N3_rdm) dans les différentes
barres du treillis en fonction de la théorie des poutres.

Ces calculs permettent de comparer les résultats obtenus par le

programme avec les valeurs théoriques prédites par la théorie des
poutres, ce qui permet de vérifier la validité du modèle et des calculs
effectués.

7
 TP1 : ELEMENT DE POUTRE DE BERNOULLI
K-Post-traitements : Calcul des déplacements
La section calcul des déplacements calcule les déplacements verticaux
et horizontaux de la poutre.
*deformée
- def0 = defo geotot dep 0. VERT; : Cette ligne crée une
déformation de la structure (`def0`) en utilisant la fonction
`defo`. Elle affiche la déformation initiale (`dep`) de la structure
`geotot` en vert.
- def1 = defo geotot dep 5. rouge; : Cette ligne crée une
déformation de la structure (`def1`) en utilisant la fonction
`defo`. Elle affiche la déformation déformée (`dep`) de la
structure `geotot` à l'étape 5 en rouge.
- titr 'Deformee'; : Cette ligne affiche le titre "Deformee" pour le
graphique des déformations.
- trac (def0 et def1); : Cette ligne trace les déformations initiale et
déformée sur le même graphique.

- u_mef = extr dep Mz (0 1); : Cette ligne extrait le déplacement

vertical à l'extrémité de la suspente à partir des résultats de
déplacement (`dep`).

- mess 'Deplacement vertical a l extremite de la suspente :';mess; :

Ces lignes affichent un message indiquant le déplacement
vertical à l'extrémité de la suspente.

- ch1 = CHAINE ' calcule = ' 'FORMAT' '(1pg12.5)'

u_mef', RdM = ' u_rdm 'FORMAT' '(1pg10.3)'
' (ecart:'(abs (100*(u_mef - u_rdm)/u_rdm))' %)'; : Cette
ligne crée une chaîne de caractères (`ch1`) pour afficher les
déplacements calculés (`u_mef`) et les déplacements analytiques
attendus (`u_rdm`), ainsi que l'écart relatif entre les deux.

8
 TP1 : ELEMENT DE POUTRE DE BERNOULLI
L- Post-traitements : Calcul des efforts
- sig = sigma (mo1 et mo2) (ma1 et ma2) dep; : Cette ligne calcule
les efforts internes (contraintes) (sig) dans les éléments de la
structure en utilisant les résultats de déplacement (dep) et les
propriétés matérielles.
- MZ = EXCO SIG momz; VY = EXCO SIG effy;: Ces lignes
extraient les moments fléchissants (MZ) et les forces axiales
VY) à partir des résultats des efforts internes (sig).
- XX = COOR 1 MZ; LX = PROG; ...: Ces lignes initialisent des
variables pour stocker les coordonnées (XX) des résultats des
moments fléchissant, et initialisent des variables (LX, LMZ,
LVY) pour stocker les données des résultats à afficher.

- REPETER BOU1 (NP-1) ; ... FIN BOU1;: Cette boucle parcourt

chaque nœud (`NP`) de la structure pour extraire les données des
résultats.

- EVMZef = EVOL ROUGE MANU 'X' LX 'Mz' LMZ;: Cette

ligne crée une évolution des moments fléchissants (`EVMZef`)
pour afficher les résultats calculés et les résultats analytiques
attendus.

- TITR 'Mz SUR POUTRE; DESS (EVMZef et evm_rdm)

AXES;: Ces lignes affichent le titre "Mz SUR POUTRE" et
tracent les résultats des moments fléchissants sur le même
graphique.

9
 TP1 : ELEMENT DE POUTRE DE BERNOULLI
Ces sections du code traitent les résultats des déplacements et des
efforts internes de la structure après l'analyse. Elles extraient ces
données à partir des résultats de l'analyse et les affichent
graphiquement pour une meilleure compréhension.

MAILLAGE DE LA STRUCTURE

10
 TP1 : ELEMENT DE POUTRE DE BERNOULLI

Epure des moments fléchissants

Moment fléchissant sur la poutre discrétisée en 5 éléments

11
 TP1 : ELEMENT DE POUTRE DE BERNOULLI
Moment fléchissant sur la poutre discrétisée en 10 éléments

Moment fléchissant sur la poutre discrétisée en 500 éléments

12
 TP1 : ELEMENT DE POUTRE DE BERNOULLI
Déformée de la structure poutre + treillis pour une discrétisation
en 5, 10 et 500 éléments

Effort tranchant sur la poutre discrétisée en 5 éléments

13
 TP1 : ELEMENT DE POUTRE DE BERNOULLI

Effort tranchant sur la poutre discrétisée en 10 éléments

Effort tranchant sur la poutre discrétisée en 500 éléments

Commentaire :
On peut remarquer en vert les graphes obtenus par les calculs de la rdm et
en rouge les résultats issus de la modélisation par élément fini par le
logiciel Cast3M.
Par comparaison de l’allure des épures, nous remarquons que plus on
augmente le nombre d’élément plus les résultats de la modélisation par
élément fini se rapprochent de ceux de la rdm.
Une plus grande discrétisation de la poutre conduit à un raffinement de
l’écart entre les deux éléments et à l’obtention d’une bonne précision. Dans
le cadre du nombre d’éléments qui ont font objet de comparaison dans
notre résolution ; 500 éléments nous offrent une meilleure précision.

14