Université de Nantes

UFR des Sciences et

Techniques

Calcul des Structures par Élements Finis

TP1 MATLAB: ANALYSE STATIQUE ET MODALE

AVELLANEDA MALAGON Gustavo

M1 Physique Parcours Mécanique
Mai 2017
Contenido
1. INTRODUCTION ........................................................................................................ 2
2. ÉTUDES SUR LA POUTRE ......................................................................................... 2
2.1. Étude Statique 1 : Étude d’une poutre en flexion ................................................... 2
2.1.1. Calcul de la flèche analytique.............................................................................. 3
2.1.2. Calcul de la flèche en utilisant le code des élements finis poutre, 1 élément.... 3
2.2. Analyse Modale 1 : Poutre encastré à son extrémité .............................................. 5
3. ÉTUDES SUR LE PORTIQUE ..................................................................................... 7
3.1. Étude Statique 2 : Étude de la déforme d’une structure portique ......................... 7
3.2. Analyse Modale 2 : Structure Portique ................................................................... 9
4. CONCLUSIONS..........................................................................................................10

1
1. INTRODUCTION
La méthode des élements finis, est une manière numérique de résoudre plusieurs
problèmes d’ingénierie et qui permet de déterminer une solution approchée sur un
domaine spatial. En prenant ça en compte, il est impératif de maitriser cette méthode
de calcul pour avoir un outil puissante à l’heure de faire le calcul de nos structures, soit
dans le cadre de la recherche ou de l’industrie, avec des codes tapes par nous-mêmes ou
avec des logiciels crées seulement pour ce but.

Le présent rapport est mon première approximation à la méthode de façon approfondie,

c’est-á-dire de façon académique et pas dans un bureau de calcul avec des logiciels
comme ETABS, SAP2000 ou MIDAS GEN, très connues et utilisées pour les ingénieurs
civiles. Dans la première partie, nous allons voir le développement d’un code sur
MATLAB pour l’analyse statique d’une poutre encastrée dans son extrémité et discrétisé
avec un seul élément fini poutre. Après, nous allons traiter le problème statique d’une
structure portique. Finalement, nous ferons le calcul des modes et fréquences propres
pour les structures poutre et portique discrétisés avec plusieurs éléments, en analysant
les résultats des calculs pour arriver à quelques conclusions et remarques.

2. ÉTUDES SUR LA POUTRE

Pour la solution des problemes traités dans cette rapport nous avons contruit une code
de calcul sur MATLAB qui permet son application dans le cas STATIQUE et le cas de
l’ANALYSE MODALE.

Le code de calcul est nourri avec les fichiers qui contiennent les propriétés mécaniques
et géométriques du problème traité, nommées de la façon suivant :

Information qui
Nom du Fichier Extension
contient le fichier
Tableau des connectivités ElementsN .csv
Coordonnés principales de
NodesN .csv
la structure
Forces imposées F_imp .csv
Déplacements imposées U_imp .csv
Propriétés de la Section MatGeo .csv

N : fait référence au nombre d’elements pris pour discretiser la structure.

2.1. Étude Statique 1 : Étude d’une poutre en flexion

Il s’agit de l’etude d’un poutre encastrée a son extremité soumis a un charge ponctuelle
de 300N á son extremite libre. La structure, discretisé avec un seul élement fini poutre,
a 6DDL (3DDL par noeud) :

2
 Figure 1. Poutre encastrée à son extrémité et soumis à une charge ponctuelle
Source : https://campus.mines-douai.fr/

Les données d’entrée pour commencer la résolution de ce problème sont les suivantes :

DATA VALUE
Longueur L (m) 1
Module de Young (Pa) 2,07 x 1011
Force Imposée (N) -300
Section 0,02m x 0,02m

L’extrémité encastrée est immobile, c’est-á-dire que dans cette extrémité les
déplacements sont imposées.

2.1.1. Calcul de la flèche analytique

D’abord il faut calculer la flèche analytique à l’aide de la théorie de poutres pour pouvoir
comparer avec les résultats données par le code de calcul et les valider. On trouve les
valeurs suivantes pour la flèche, la rotation, le moment et l’effort tranchant :

Formule analytique Résultat

flèche u
1 𝐹𝑙 3
𝑢= 𝑢 = −0,03632
3 𝐸𝐼
rotation θ 1 𝐹𝑙 2
𝜃= 𝜃 = −0,05448
2 𝐸𝐼
moment fléchissant max 𝑀 = −𝐹𝐿 300N
effort tranchant 𝑇 = −𝐹 300N

2.1.2. Calcul de la flèche en utilisant le code des

élements finis poutre, 1 élément
Une fois obtenus les résultats théoriques, il est moment de mettre en fonctionnement le
code écrit pendant les séances du TP pour voir l’exactitude de la méthode FEM. Pour
démarrer le calcul, un fois sur MATLAB il faut lancer le fichier Program.m, et il ensuite
nous affichera les graphs et résultats suivants :

3
 Figure 2. Calcul de la flèche en utilisant le code
Source : Auteur

Les valeurs affiches par MATLAB sont :

>> Program Où :
U_inc =
0 U_inc : Déplacements inconnues pour le
-0.0363 nœud de l’extrémité libre (nœud 2).
-0.0545
R_XYM : Vecteur des valeurs de la
R_XYM = réaction en X, Y et le Moment à
0 l’encastrement.
300.0000
300.0000 NTM : Vecteur avec les valeurs des
sollicitations internes (l’effort Normal,
NTM = Tranchant et le Moment Fléchissant)
0 pour le nœud 2 de la structure
-300
0

Ce qui confirme la validité de notre code de calcul EF.

A continuation l’affichage des diagrammes de sollicitations internes (effort tranchant et

moment fléchissant) fait par notre code de calcul :

4
 Figure 3. Diagrammes des sollicitations internes
Source : Auteur

2.2. Analyse Modale 1 : Poutre encastré à son

extrémité
A continuation, on va montrer le travail réalisé pour l’étude des 5 premiers modes de
vibration de la poutre décrite dans l’étude statique 1, et on va montrer l’allure obtenu
avec le logiciel pour chaque mode de vibration associé aux 5 premières fréquences
propres du système. Il faut tenir en compte que la masse volumique de l’acier constituant
la poutre est ρ = 7800 kg/m3. L’étude a été réalisé en discrétisant la poutre á l’aide de 10
éléments, 100 éléments et 1000 éléments. Les résultats de cette analyse sont résumés
dans le graph suivant :

5
 Figure 4. Poutre discrétisé á l’aide de 10 éléments

Figure 5. Poutre discrétisé á l’aide de 100 éléments

Figure 5. Poutre discrétisé á l’aide de 1000 éléments

Comme on peut voir dans les graphiques, la déviation des valeurs augmente avec les
fréquences plus grandes, mais le plus grand est le nombre d’éléments utilisé pour la
discrétisation, le plus la valeur de la fréquence tiens à converger. Pour les fréquences
plus petits, la convergence s’atteint avec moins nombre d’éléments.

6
3. ÉTUDES SUR LE PORTIQUE

3.1. Étude Statique 2 : Étude de la déforme d’une structure

portique
Il s’agit de l’etude d’un portique encastrée a son extremité soumis a un charge ponctuelle
de 19.62N sur le nœud 41 en direction x positive. La structure, discretisé avec de
élements finis poutre, a un appui simple á son extremité gauche et un appui type
« roller » a son extremité droite. Le affichage de la structure est le suivante :

Figure 6. Structure Portique discrétisé a l’aide de 60 éléments

Les données d’entrée pour commencer la résolution de ce problème sont les suivantes :

DATA VALUE
Longueur L (m) 0.5
Hauteur (m) 0.4
Module de Young (Pa) 74 x 109
Force Imposée (N) 19,62 suivant x positives
Section 40 mm x 5 mm

Avec les chargements imposés et les conditions limites dans les fichiers reçus par email,
on trouve la déformé suivante :

7
 Figure 7. Structure Portique déformé

La déformé est cohérente avec la déformée qualitative qui peut être produit par le
chargement imposé. Elle est aussi cohérente avec les calculs analytiques de la RDM.

Les valeurs des réactions trouvées sont :

Pour le nœud 1 : Rx=-19,62N Ry=-15.696N

Pour le nœud 61 : Ry=15,696N

Pour les diagrammes des sollicitations internes, on trouve le suivant :

0.07848*10 Nm

0.07848*10 Nm

Figure 8. Moment Fléchissant

8
 0.1962*10 N

0.1569*10 N

Figure 9. Effort Tranchant

3.2. Analyse Modale 2 : Structure Portique

A continuation, on va montrer le travail réalisé pour l’étude des 5 premiers modes de
vibration de la Structure Portique décrite dans l’étude statique 2, et on va montrer
l’allure obtenu avec le logiciel pour chaque mode de vibration associé aux 5 premières
fréquences propres du système. Il faut tenir en compte que la masse volumique de
l’aluminium constituant la structure est ρ = 2700 kg/m3.

Les résultats de cette analyse sont résumés dans la Figure 10, avec les fréquences trouvés
pour chaque mode analysé.

Figure 10. 5 premières modes propres pour l’analyse modale de la Structure Portique

9
4. CONCLUSIONS
 Pour l’analyse modale de la Structure Portique, on a trouvé que le risque de mise
en résonance est très haut, et ça peut être corrélé avec le type de matériau et la
géométrie de la section. Avec les fréquences plus hautes, les déplacements
trouvés sont excessives. Les fréquences trouvées pour la configuration de la
Structure Portique sont plus petites que pour la poutre de l’analyse modale 1,
c’est-à-dire que le risque de mise en résonance est plus haut pour le portique
étudié parce que sa fréquence propre pour le première mode est très petite, et il
est facile de la mettre en résonance avec une force aussi petite.
 La méthode des éléments finis est un outil mathématique très puissant pour le
calcul des structures de n’importe quel type. On a bien constaté qu’elle nous
permet de connaitre bien les déplacements dus à différentes types des
chargements sur nos structures, aussi comme les forces et efforts internes.
 Comme méthode d’analyse académique les TP Matlab sont très graphiques et
nous permettent de maitriser les principes et la théorie nécessaires pour pouvoir
utiliser après des logiciels commerciaux comme ADINA, ABACUS,
HYPERMESH, CATIA, SAP2000, entre outres. Le calcul des structures complexes
sur Matlab peut prendre beaucoup de temps et il est difficile de contrôler les
erreurs dans le code.
 Pour les analyses sur la poutre (statique et modale) on a trouvé les mêmes valeurs
analytiques et numériques, ce qui nous a permis de valider notre code.
 Pour les analyses sur la structure portique, la partie d’analyse statique est
d’accord à la théorie et on a trouvé les mêmes valeurs analytiques et numériques.
Pour l’analyse modale, ils existent des valeurs un peu différents, ce qui nous
donne l’idée de continuer à déboguer le code pour pouvoir l’utiliser comme il
faut. C’est l’affichage des diagrammes des modes sur Matlab qui m’a posé le plus
des problèmes à coder.

10