Génie des systèmes industriels

Travaux Pratiques
ÉLÉMENTS FINIS

PROFESSEUR :
Jean-Mathieu Mencik

ÉTUDIANT :
Quoc Tan NGUYEN
4TD1 - 4TP2

24 juillet 2024
Table des matières
1 Introduction 3

2 Description de l’idée 4
2.1 Caractéristique et Dimension . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4

3 Modélisation 5
3.1 Matrice de rigidité . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
3.2 Matrice de masse . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
3.3 Matrice de rotation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
3.4 Assemblage . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7

4 Méthode d’analyse de la réponse harmonique 7

4.1 Analyse harmonique de la déformée . . . . . . . . . . . . . . . . . . . . . . . . . . 8
4.2 Mise en place du maillage . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10

5 Réduction de modèle 13

6 Conclusion 14

1
Table des figures
Fig 1.L’avion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3
Fig [Link]̀le simplifié . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3
Fig [Link] des sous-structures (peint à la main) . . . . . . . . . . . . . . . . . . 4
Fig [Link]́placement du S3 (au-dessus) et du S6 (en dessous) . . . . . . . . . . . . . . . 8
Fig [Link]́placement du S3 (au-dessus) et du S6 (en dessous) . . . . . . . . . . . . . . . 9
Fig [Link]́formée de la structure en f = 7.4 Hz . . . . . . . . . . . . . . . . . . . . . . . 9
Fig [Link]́formée de la structure en f = 11.9 Hz et f = 24.7 Hz . . . . . . . . . . . . . . 10
Fig [Link]́formée de la structure en f = 10.3 Hz . . . . . . . . . . . . . . . . . . . . . . . 10
Fig [Link]́diaire 2 (Noir), Intermédiaire 1 (Rouge), Grossier (Bleu) du S3 . . . . . . 11
Fig [Link]́ponse harmonique pour maillage fin (Noir) et maillage grossier (Bleu) du S3 . 12
Fig [Link]́ponse harmonique pour maillage fin (Noir) et intermédiaire 1 (Rouge) du S3 . 12
Fig [Link]́ponse harmonique pour maillage fin (Noir) et intermédiaire 2 (Vert) du S3 . . 13

2
1 Introduction

Dans ce rapport, on va faire un rechercher a le comportement dynamique d’un avion. Ce étude

s’appuie sur un modèle simplifié d’un avion qui est composé de 7 sous-structures. Ensuite, nous
allons construire ses propriétés et les simuler sur la base du logiciel Matlab, dans le but d’étudier
ses changements lorsqu’il est appliqué sous force opposée.

Figure 1 – L’avion

Ce rapport est construit selon les idées et le code Matlab suggérés par le professeur Jean-
Mathieu Mencik. Il sera présenté dans l’ordre suivant :
— Modéliser la structure
— Tracer les réponses en fréquence en différents points de mesure
— Analyser la sensibilité des solutions
— Mettre en évidence les erreurs associées aux maillages grossier et intermédiaires par rapport
au maillage fin
— Tester la méthode de réduction de modèles basée sur les modes de vibrations

Figure 2 – Modèle simplifié

3
2 Description de l’idée

Un modèle simplifié d’avion est composé de 7 sous-structures Sj (j = 1,2,..,7), lesquelles sont

elles-mêmes constituées de plusieurs éléments de poutre à 2 nœuds et 6 degrés de liberté par
nœud. Par ailleurs, le modèle étudié dans ce rapport comprend deux types de sous-structures,
celles qui forment un corps à section cylindrique (S1 et S2 ) et les autres à section rectangulaire
(de S3 à S7 ).

Figure 3 – Assemblage des sous-structures (peint à la main)

2.1 Caractéristique et Dimension

1. Caractéristiques géométriques
— La longeur L
— Le nombre d’élements N
— La longueur des éléments Le = L/N
2. Caractéristiques matériaux
— Module de Young : E = 70 (Gpa)
— Coefficient de poisson : ν = 0.3
E
— Module de cisaillement : G = 2(1+v)
— La densité : ρ = 2700 (kg/m3 )
3. Caractéristiques géométriques des sous-structures
Section cylindrique
— Rayon extérieur et intérieur : Rext et Rint

4
 — Épaisseur : e
R4 −R4
— Moment quadratique selon ⃗y , ⃗z : Iz = Ix = π. ext 4 int
— Moment quadratique polaire : Jpol = Iz + Iy

Structure Rayon extérieur(m) Rayon intérieur(m) Longueur(m)

S1 1.001 1 10
S2 1.001 1 20

Section rectangulaire
— Largeur, hauteur : b et h
— Section : S = b.h
3
— Moment quadratique selon ⃗y , ⃗z : Iz = Ix = π. b.h
12
— Moment quadratique polaire : Jpol = Iz + Iy

Structure Hauteur(m) Largeur(m) Longueur(m)

S3 0.1 0.6 20
S4 0.1 0.6 20
S5 0.6 0.1 5
S6 0.1 0.6 10
S7 0.1 0.6 10

Les valeurs retenues sont référencées à partir des consignes de l’enseignant. Concernant les
valeurs de hauteur et de largeur, on choisit des valeurs identiques mais différentes dans S5 car
la position de S5 est verticale. Dans le but de faire de LLElement
P outre
est un entier, pour que la force
soit appliquée sur un nœud de l’aile.

3 Modélisation

3.1 Matrice de rigidité

La matrice de masse de barre en traction/compression :

" #
ρSL e 2 1
[Me1oc ]b = (1)
6 1 2

La matrice de masse d’une poutre en torsion :

" #
ρJθ Le 2 1
[Meloc ]t = (2)
6 1 2

La matrice de masse d’une poutre en flexion dans le plan (xe , y e ) ), faisant intervenir les
déplacements/rotations (v1e , θz1
e ) et (v e , θ e ) :
2 z2
 
156 22Le 54 −13Le
ρSLe  22Le 4 (Le )2 13Le −3 (Le )2 
 
e
[M1◦c ]fxy = 
e e 
 (3)
420   54 13L 156 −22L 
−13Le −3 (Le )2 −22Le 4 (Le )2

5
La matrice de rigidité d’un élément de poutre en flexion dans le plan (z e , xe ), faisant intervenir
les déplacements/rotations w1e , θy1
e et w2e , θy2
e



:
 
12 6Le −12 6Le
6Le 4 (Le )2 −6Le 2 (Le )2 
 
e EIy 
[Kloc ]fzx =   (4)
(Le )3 
 −12 −6L
e 12 −6Le  
6Le 2 (Le )2 −6Le 4 (Le )2

3.2 Matrice de masse

La matrice de masse de barre en traction/compression :

" #
e ρSLe 2 1
[M1oc ]b = (5)
6 1 2
La matrice de masse d’une poutre en torsion :
" #
ρJθ Le 2 1
[Meloc ]t = (6)
6 1 2
La matrice de masse d’une poutre en flexion dans le plan (xe , y e ) ), faisant intervenir les
déplacements/rotations (v1e , θz1
e ) et (v e , θ e ) :
2 z2
 
156 22Le 54 −13Le
ρSLe  22Le 4 (Le )2 13Le −3 (Le )2 
 
e
[M1◦c ]fxy = 
e e 
 (7)
420   54 13L 156 −22L 
−13Le −3 (Le )2 −22Le 4 (Le )2
La matrice de masse d’un élément de poutre en flexion dans le plan (z e , xe ), faisant intervenir
e
les déplacements/rotations w1e , θy1 et w2e , θy2
e



:
 
156 22Le 54 −13Le
e e )2 e e )2 
ρSLe 
 
e 22L 4 (L 13L −3 (L
[M1◦c ]fzx =   (8)
420   54 13Le 156 −22Le  
−13Le −3 (Le )2 −22Le 4 (Le )2

3.3 Matrice de rotation

En raison des différences entre certaines structures du système de coordonnées local et du

système de coordonnées global, il est impératif de calculer la matrice de rotation.
Matrice de rotation autour l’axe ⃗z :
 
cos θ − sin θ 0
Rz (θ) =  sin θ cos θ 0  (9)
 

0 0 1
L’angle de rotation de chaque structure suivant :

S1 S2 S3 S4 S6 S7
π
θ 0 π 2 − π2 π
2 − π2

6
 Matrice de rotation autour l’axe ⃗y :
 
cos θ 0 − sin θ
Ry (θ) =  0 1 0 (10)
 

sin θ 0 cos θ

Pour S5 , on va appliquer une rotation autour l’axe ⃗y : θ = − π2 .

Pour les poutres étudies , la matrice de rotation est :
 
[r] 0 0 0
 
 0 [r] 0 0 
[Re ] = 
  (11)
 0 0 [r] 0 

0 0 0 [r]

La matrice [r] est Ry (θ) ou Rz (θ) (dépend de la poutre).

3.4 Assemblage

Convertir la matrice de rigidité locale en système de base global :

[Ke ] = [Re ]T [Keloc ] [Re ] (12)

Assembler les matrices Ke en une matrice K qui peut caractériser l’ensemble du système :
 
Ksj (nume , nume ) = [Ke ] (13)

Où, nume est un vecteur représentant la position de déplacement de chaque élément .

Enfin, les poutres ensemble :
X 7
 
[K] = Kg1sj (14)
j=1

h i
 
Kglsj numsj , numsj = Ksj (15)
On va adopter la même approche avec la matrice de masse.

4 Méthode d’analyse de la réponse harmonique

D’après le cours VIBRATIONS DES STRUCTURES (B. Bergeot & V. Denis), on va étudier
ce problème sous l’influence de forces harmoniques.
L’équation de mouvement :
Mû¨ + Cû + Kû = Fejωt (16)
La solution stationnaire : û = Ûejωt qui est introduite dans l’équation du mouvement.

−ω 2 M + jωC + K Û = F


(17)

La matrice rigidité dynamique :

D = −ω 2 M + jωC + K (18)

7
La matrice des amplitudes de vibration :

Û = D−1 F (19)

Ici, on a :

C = a.M + b.K = 0.01M

 
0
 . 
 .. 
 
 
−F0 
 
 0 
 
 . 
 .
[F] =  .  (20)
 
 0 
 
F 
 0 
 . 
 .. 
 
0

4.1 Analyse harmonique de la déformée

Pour étudier le déplacement, on choisit le nombre d’éléments N = 10 et l’intervalle de chan-

gement de fréquence f = [0,50], ∆f = 0.1 et l’amplitude de la force appliquée F0 = 1.

Figure 4 – Déplacement du S3 (au-dessus) et du S6 (en dessous)

8
 En examinant la figure, on observe l’émergence de pics qui traduisent la résonance de la sous-
structure. Par ailleurs, il est notable que l’utilisation de différents points de mesure conduit à
des fréquences générant la résonance dans une sous-structure sans affecter les autres.

Figure 5 – Déplacement du S3 (au-dessus) et du S6 (en dessous)

De plus, on peut voir qu’aux basses fréquences, notre modèle représente clairement l’état vi-
bratoire de la structure. Cependant, aux hautes fréquences (supérieure à 300 Hz), les vibrations
de la structure ne sont pas distinctement visibles. Ceci est attribuable à la finesse insuffisante
de notre maillage pour décrire de manière exhaustive les oscillations du système.

Figure 6 – Déformée de la structure en f = 7.4 Hz

9
 Figure 7 – Déformée de la structure en f = 11.9 Hz et f = 24.7 Hz

En examinant les trois images ci-dessus, on peut clairement distinguer des différences signifi-
catives entre elles. Nous en déduisons ainsi que chaque maillage n’est précise que dans une plage
spécifique de fréquences.

4.2 Mise en place du maillage

Nous sommes actuellement en phase de détermination du nombre optimal d’éléments pour

la plage de fréquences que nous avons définie, allant de la fréquence 0 à la fréquence maximale
fmax .
Le choix de ce nombre d’éléments revêt une importance capitale pour la création d’un maillage
finement optimisé. Une fois cette étape accomplie, nous serons en mesure de comparer les
résultats obtenus avec les maillages fin (intermédiaire 2), moyen (intermédiaire 1) et
grossier.

Figure 8 – Déformée de la structure en f = 10.3 Hz

L’équation de vibration :
∂4ω
EI + ρ.S.ω = 0, x ∈]0, L[ (21)
∂x4
ρS
On pose le nombre d’ondes β = ω 2 . EI On a donc :

2πC 2π 2π EI 1
λ= = =√ ( )4 (22)
ω β ω ρS

Le critère de maillage nécessaire pour être précis jusqu’une certaine fréquence de coupure fc . Il

10
faut au moins 4 élements par longeur d’onde.
r
2π EI 1
λ(fmax ) = ( )4 (23)
fmax ρS
Longeur des élements :
λ(fmax )
Le = (24)
4
Nombre d’élement : r r
e 8fmax EI 1 2π EI 1
Nopt =L ( )4 = ( )4 (25)
π ρS f ρS
Ensuite, nous continuerons à considérer le changement de fréquence de 0 à 1000 Hz. Par
conséquent, nous appliquerons fmax = 1000 Hz pour calculer le nombre optimal d’éléments
pour les éléments. Par ailleurs, à partir du nombre optimal d’éléments, on peut en déduire
les valeurs des maillages intermédiaire 2 (70 % Nopt ) , intermédiaire 1 (40 %Nopt ) et grossier
(10 %Nopt ). Nous pouvons obtenir Nopt en exécutant le code Matlab.

Sous-structure Nopt Intermédiaire 2 Intermédiaire 1 Grossier

S3 84 64 36 8
S4 84 64 36 8
S5 9 8 4 1
S6 42 32 16 4
S7 42 32 16 4

Tableau du nombre d’éléments

Figure 9 – Intermédiaire 2 (Noir), Intermédiaire 1 (Rouge), Grossier (Bleu) du S3

Pour chaque maillage distinct, il est possible de déterminer la fréquence de coupure, à partir de
laquelle la précision du maillage diminue par rapport au maillage fin. En mesurant la fréquence
de coupure fc nous pourrons ensuite déterminer le nombre minimal d’éléments par longueur
d’onde Nc associé.

λ(ωc ) 2π EI 1
Nc = = √ ( )4 (26)
Le Le 2πfc ρS
Ensuite, on va rechercher la fréquence de coupure pour chaque maillage.

11
— Maillage grossier et fin :

Figure 10 – Réponse harmonique pour maillage fin (Noir) et maillage grossier (Bleu) du S3

On a : Le = L8 = 2.5 (m) et on voir que le courbe bleu décroche à partir de fc = 17.46 Hz.
Donc, Nc = 7 (satisfaire aux conditions).
— Maillage intermédiaire 1 et fin :

Figure 11 – Réponse harmonique pour maillage fin (Noir) et intermédiaire 1 (Rouge) du S3

L
On a : Le = 36 = 0.56 (m) et on voir que le courbe rouge décroche à partir de fc = 122.04
Hz. Donc, Nc = 12 (satisfaire aux conditions).

12
 — Maillage intermédiaire 2 et fin :

Figure 12 – Réponse harmonique pour maillage fin (Noir) et intermédiaire 2 (Vert) du S3

L
On a : Le = 64 = 0.312 (m) et on voir que le courbe vert décroche à partir de fc = 266.6
Hz. Donc, Nc = 15 (satisfaire aux conditions).
A travers les 3 cas ci-dessus, nous remarquons qu’avec une fréquence plus élevée, le nombre
minimal d’éléments par longueur d’onde augmente. Cette observation est cohérente car à des
fréquences élevées, les longueurs d’onde sont plus petites, nécessitant ainsi la construction d’un
maillage plus fin.
Grâce aux résultats du nombre optimal d’éléments reçus, nous pouvons voir qu’il y a une
erreur entre eux, ce qui rend notre résultat final pas le plus parfait. Mais pour éviter de perdre
du temps et des efforts à calculer, nous aborderons une méthode présentée ci-dessous.

5 Réduction de modèle

Il s’agit d’une technique classique de réduction basée sur les modes de vibration de la structure.
Nous devons considérer le problème des valeurs propres du système :

[K]X = Ω2 [M]X (27)

Pour calculer le vecteur des déplacements, parmi tous les n vecteurs propres [Xi ]on ne choisit
que les [M] premiers vecteurs.
Le principe de décomposition modale, u(t) = [X]α(t) , on a :
M
X
u(t) ≈ u
e (t) = Xj α
ei (t) (28)
i=1

Ou, [X]
e est la matrice réduite des modes de vibration
α
e(t) le vecteur des coordonnés généralisées
L’équation de mouvement :

[M] ¨
f α̃(t) + [C] ˙
e α̃(t) + [K]
e α(t)
e = F(t)
e (29)

La solution :
α e −1 F
e = [D] e ⇒u
e (t) = [X]
e α(t)
e (30)

13
 e = −ω 2 [M]
Avec, [D] f + iω[C]
e + [K]
e

6 Conclusion

Ces travaux pratiques ont été une opportunité précieuse pour renforcer ma compréhension
de la méthode des éléments finis, particulièrement dans son application à l’analyse statique et
dynamique des structures.
En parallèle, j’ai pu observer des applications concrètes de la méthode des éléments finis.
D’autre part, le choix du maillage doit concilier les exigences géométriques de la structure avec
une optimisation du temps de calcul. De ce fait, je considère crucial de respecter rigoureusement
les critères de maillage.

14