REPUBLIQUE DU BENIN

************
MINISTERE DE L’ENSEIGNEMENT
SUPERIEURE ET DE LA FORMATION
TECHNIQUE
********************
ECOLE POLYTECHNIQUE D’ABOMEY-
CALAVI (EPAC)
*****************

Matière : Modélisation des systèmes et simulation

Filière : Génie Civil (4ème année)
GROUPE 1

THEME : RESOLUTION D’UN PROBLEME EN DEFORMATION

PLANE : CAS D’UN BARRAGE HYDRAULIQUE

Membres du Groupe :
1 ABAI Larissa 13 FAYOMI Shadrak
2 ABLET Benito 14 FONGNIKIN Thierry
3 AGOSSOU Emmanuel 15 GBENOU Ulrich
4 AHOKPE Othniel 16 GBETISANON Malchade
5 AKPI DIOUF 17 IBRAHIM A. Youssoufa
6 ALLOHOUTADE Ferdinand 18 LEDJOU Lucien
7 AMADOU Saratou 19 MIKAILA Moubarak
8 AVOSSE-GUEHOU Pygas 20 SAHOSSI Audry
9 ATIMANDJEGANGNI Wilfried 21 SOYITO Prosper
10 DEGAN Belvida 22 WEINSOU Juste
11 DELE Josias 23 ZOUNTCHEME Spéros
12 DOSSOU Chance

Sous la supervision de :
Emmanuel OLODO, Ing PhD Professeur Titulaire des Universités
CAMES
Assistant : Ing DATCHOSSA Abbas

P a g e 1 | 39
2022-2023
Table des matières
Liste des figures..................................................................................................... ii

Liste des tableaux ................................................................................................. iii

INTRODUCTION ................................................................................................. 1

I. SUJET ............................................................................................................. 2

II. Hypothèse........................................................................................................ 3

III. Discrétisation de maillage ............................................................................ 4

IV. Rappel des conditions de déformations planes d’un élément triangulaire à

03 nœuds................................................................................................................ 6

V. Matrice des déformations élémentaires B ....................................................... 7

VI. Matrice de rigidité élémentaire Re .............................................................. 9

VII. La matrice de rigidité élémentaire 𝑹𝒏′ ..................................................... 10

VIII. La matrice de rigidité globale .................................................................... 12

IX. La matrice des forces F .............................................................................. 12

X. La matrice des déplacements F ..................................................................... 12

XI. Les contraintes et déformations ................................................................. 13

Conclusion ........................................................................................................... 14

Annexes ............................................................................................................... 15

P a g e i | 39
 Modélisation des systèmes et simulation

Liste des figures

Figure 1:Barrage et schéma statique du problème ................................................ 3
Figure 2 Repère global de l’élément ..................................................................... 3
Figure 3:Discrétisation du domaine en 100 éléments ........................................... 4
Figure 4:Notation des nœuds dans le repère global On remarque aisément que
notre élément contient 66 nœuds. .......................................................................... 5
Figure 5: Représentation des forces agissantes sur le système ............................. 5
Figure 6:Elément de type 1 ................................................................................... 6
Figure 7:Elément de type 2 ................................................................................... 6

P a g e ii | 39

EPAC-GC4-2022-2023 GROUPE 1
 Modélisation des systèmes et simulation

Liste des tableaux

P a g e iii | 39

EPAC-GC4-2022-2023 GROUPE 1
 Modélisation des systèmes et simulation

INTRODUCTION
La méthode des éléments finis (MEF) est une méthode numérique couramment
utilisée pour résoudre des problèmes d'ingénierie et de physique. Elle est basée
sur le principe de discrétisation du domaine de calcul en éléments finis, sur
lesquels une solution approchée est déterminée. La MEF est largement utilisée
dans de nombreux domaines, tels que la mécanique des structures, la mécanique
des fluides, l'électromagnétisme, la thermique, etc.
La MEF consiste à subdiviser le domaine de calcul en un certain nombre
d'éléments finis, qui sont des polygones ou des polyèdres simples, tels que des
triangles, des quadrilatères, des tétraèdres ou des hexaèdres. Sur chaque élément
fini, une solution approchée est déterminée à l'aide d'une fonction d'interpolation,
qui est souvent une combinaison linéaire de fonctions de base. La solution globale
est obtenue en assemblant les solutions locales de chaque élément fini.

Grâce au calcul, même simplifié, l'ingénieur peut tester plusieurs configurations

pour optimiser le comportement d'un modèle à une prestation donnée. Cela évite
de multiplier les prototypes et les essais tests réels, les supports physiques ne
servent plus à chercher une solution, ils permettent de la valider. Le calcul de
conception dans les bureaux d'études, c'est sans doute le plus répandu car grâce
aux outils de calcul simplifié dont disposent les logiciels de CAO modernes, la
simulation numérique fait partie des outils de conception pour obtenir un
comportement défini à priori qui détermine le dimensionnement, donc le dessin,
des pièces mécaniques

Il s’agira de calculer les contraintes et les déformations au sein d’un barrage de

forme triangulaire.

Pour a résolution du problème ci-dessous nous nous servirons du logiciel Matlab.

P a g e 1 | 39

EPAC-GC4-2022-2023 GROUPE 1
 Modélisation des systèmes et simulation

I. SUJET
On considère un barrage constitué d'un matériau élastique isotrope et posé sur un
massif de sol semi infini (voir figures ci-dessous). On suppose que le barrage est
encastré à sa base AB dans un massif infiniment rigide.
En adoptant un maillage constitué d'éléments triangulaires à 3 nœuds et en
résolvant le problème en déformations planes, déterminer les déplacements, les
déformations et les contraintes dans le barrage pour les sollicitations suivantes :
➢ Le côté AC est soumis à une pression constante Pm=90 kg/m2 ;

➢ La masse volumique du matériau du barrage 𝜌m=2500 kg/m3 .

Dans tout le problème, on prendra E = 2000 MPa ; ν = 0,2.

P a g e 2 | 39

EPAC-GC4-2022-2023 GROUPE 1
 Modélisation des systèmes et simulation

Figure 1:Barrage et schéma statique du problème

II. Hypothèse
➢ Repère global de notre élément (barrage)

Figure 2 Repère global de l’élément

➢ Degré de liberté
Ici, nous sommes en présence d’un élément de volume tridimensionnel. A cet
élément est appliquée une pression constante suivant l’axe (ox) et une force de
volume suivant l’axe (oy).
Sachant très bien que les forces créent les sollicitations, qui à leur tour,
engendrent les déplacements qui conduisent aux déformations et aux contraintes,
nous n’aurons pas de déplacement suivant l’axe (oz). De ce fait, deux
déplacements (un déplacement horizontal notée u et un déplacement vertical
notée v) seront observés sur notre élément de structure.
Il s’ensuit qu’à chaque nœud nous aurons deux degrés de liberté et le
problème sera résolu dans le plan (oxy) (avec une épaisseur d’une unité)
puisqu’on observe la même chose suivant l’axe (oz).

P a g e 3 | 39

EPAC-GC4-2022-2023 GROUPE 1
 Modélisation des systèmes et simulation

III. Discrétisation de maillage

Nous allons discrétiser la surface ABC en élément fini triangulaire à trois
(03) nœuds. Un élément fini considéré, aura une base de 2 m ; 1 m de hauteur et
une épaisseur d’une unité.

Ce qui nous conduit au maillage ci- après :

Figure 3:Discrétisation du domaine en 100 éléments

P a g e 4 | 39

EPAC-GC4-2022-2023 GROUPE 1
 Modélisation des systèmes et simulation

Figure 4: Représentation des forces agissantes sur le système

Figure 5:Notation des nœuds dans le repère global On remarque aisément que
notre élément contient 66 nœuds.
➢ Différents types d’éléments constituants le domaine.
On distingue deux types d’éléments :
▪ Les éléments triangulaires de type 1 dont l’angle droit est vers le
Haut

P a g e 5 | 39

EPAC-GC4-2022-2023 GROUPE 1
 Modélisation des systèmes et simulation

Figure 6:Elément de type 1

▪ Les éléments triangulaires de type 2 dont l’angle droit est vers le bas

Figure 7:Elément de type 2

IV. Rappel des conditions de déformations planes

d’un élément triangulaire à 03 nœuds
Matrice des constantes élastiques (Rigidité D)

P a g e 6 | 39

EPAC-GC4-2022-2023 GROUPE 1
 Modélisation des systèmes et simulation

Dans ce problème, nous somme en déformation plane donc utilisons la formule

L’hypothèse élaborée par KIRCHOFF permet d’obtenir l’expression suivante
pour la rigidité D de :
𝛎
𝟏 𝟎
𝟏−𝛎
𝐄 𝛎
𝐃= × 𝟏 𝟎
𝟏 − 𝟐𝛎 𝟏 − 𝛎
𝟏 − 𝟐𝛎
𝟎 𝟎
[ 𝟐(𝟏 − 𝛎)]

Avec 𝐄 = 𝟐𝟎𝟎𝟎 𝐌𝐏𝐚 et 𝛎 = 𝟎, 𝟐

1
1 0
4
2000 1
D= × 1 0
2 4
1−
5 3
[0 0
8]

2000 × 5 1 8 2 0
D= × × [ 2 8 0]
3 8
0 0 3

2500 8 2 0
D= × [2 8 0]
6
0 0 3

𝟏𝟐𝟓𝟎 𝟖 𝟐 𝟎
𝐃= × [𝟐 𝟖 𝟎]
𝟑
𝟎 𝟎 𝟑

V. Matrice des déformations élémentaires B

Elément triangulaire 1 :
La matrice des déformations élémentaires B a pour expression :
0 0
𝑋 (2) 𝑌 (0)
0 1
P a g e 7 | 39

EPAC-GC4-2022-2023 GROUPE 1
 Modélisation des systèmes et simulation

Calcul de l’air
1
𝐴2 = (𝑋21 . 𝑌31 + 𝑋31 𝑌21 )
2
1
𝐴2 = (2 × 1 + 0 × 1)
2
𝐴2 = 1
Calcul de B2

1 𝑌23 0 𝑌31 0 𝑌12 0

𝐵1 = [ 0 𝑋32 0 𝑋13 0 𝑋21 ]
2𝐴 𝑋 𝑌23 𝑋13 𝑌31 𝑋21 𝑌12
32

1 −1 0 1 0 0 0
𝐵1 = [ 0 −2 0 0 0 2]
2
−2 −1 0 1 2 0

Elément triangulaire 2
2 1
𝑋 (0) 𝑌 (1)
2 0
Calcul de l’air
1
𝐴2 = (𝑋21 . 𝑌31 + 𝑋31 𝑌21 )
2
1
𝐴2 = (2 × 1 + 0 × 1)
2
𝐴2 = 1
Calcul de B2

1 𝑌23 0 𝑌31 0 𝑌12 0

𝐵2 = [ 0 𝑋32 0 𝑋13 0 𝑋21 ]
2𝐴 𝑋 𝑌23 𝑋13 𝑌31 𝑋21 𝑌12
32

1 1 0 −1 0 0 0
𝐵2 = [0 2 0 0 0 −2]
2
2 1 0 −1 −2 0

P a g e 8 | 39

EPAC-GC4-2022-2023 GROUPE 1
 Modélisation des systèmes et simulation

VI. Matrice de rigidité élémentaire Re

La matrice de rigidité élémentaire est donnée par : 𝑅𝑒 = ∫ 𝐵𝑡 𝐷𝐵𝑑Ω𝑒

Elément triangulaire 1 :
On a :
𝑅1 = 𝐵1𝑡 . 𝐷. 𝐵1 . 𝐴𝑒
−1 0 −2
0 −2 −1
625 1 8 2 0 −1 0 1 0 0 0
0 0
𝑅1 = × [2 8 0] × [ 0 −2 0 0 0 2]
6 0 0 1
0 0 3 −2 −1 0 1 2 0
0 0 2
[0 2 0]

20 10 −8 −6 −12 −4
10 35 −4 −3 −6 −32
625 −8 −4 8 0 0 4
𝑅1 =
6 −6 −3 0 3 6 0
−12 −6 0 6 16 0
[ −4 −32 4 0 0 32 ]

Elément triangulaire 2 :
On a 𝑅2 =𝐵2𝑡 D𝐵2 A
D’où
1
𝑅2 = ∗
2
1 0 2
0 2 1
1 0.25 0 1 0 −1 0 0 0
1 2000 – 1 0 0
[0.25 1 0 ] 𝑥 [0 2 0 0 0 −2]
2 0.6 0 0 –1
0 0 0.375 2 1 0 −1 −2 0
0 0 –2
[ 0 –2 0 ]

P a g e 9 | 39

EPAC-GC4-2022-2023 GROUPE 1
 Modélisation des systèmes et simulation

20 10 −8 −6 −12 −4
10 35 −4 −3 −6 −32
625 −8 −4 8 0 0 4
𝑅2 =
6 −6 −3 0 3 6 0
−12 −6 0 6 16 0
[ −4 −32 4 0 0 32 ]

VII. La matrice de rigidité élémentaire 𝑹′𝒏

Elle concerne chaque élément et sa formule est la suivante :
𝑅𝑒′ = 𝑇𝑒𝑡 𝑅𝑒 𝑇𝑒
Elément triangulaire 1 :
On a : 𝐑′𝟏 = tT R1 T
Avec :

20 10 −8 −6 −12 −4
10 35 −4 −3 −6 −32
625 −8 −4 8 0 0 4
• 𝑅1 =
6 −6 −3 0 3 6 0
−12 −6 0 6 16 0
[ −4 −32 4 0 0 32 ]
1 0 0 0 0 0 0 0
0 1 0 0 0 0 0 0
0 0 0 0 0 0 1 0
• 𝑇1 = (la matrice de localisation)
0 0 0 0 0 0 0 1
0 0 1 0 0 0 0 0
[0 0 0 1 0 0 0 0]

1 0 0 0 0 0
0 1 0 0 0 0
0 0 0 0 1 0
0 0 0 0 0 1
• 𝑇1𝑡 =
0 0 0 0 0 0
0 0 0 0 0 0
0 0 1 0 0 0
[0 0 0 1 0 0]

P a g e 10 | 39

EPAC-GC4-2022-2023 GROUPE 1
 Modélisation des systèmes et simulation

𝐑′𝟏 =
20 10 −8 −6 −12 −4
10 35 −4 −3 −6 −32 1 0 0 0 0 0 0 0
−12 −6 0 6 12 0 0 1 0 0 0 0 0 0
625 −4 −32 4 0 0 32 0 0 0 0 0 0 1 0
𝑥 𝑥 °
6 0 0 0 0 0 0 0 0 0 0 0 0 0 1
0 0 0 0 0 0 0 0 1 0 0 0 0 0
−8 −4 8 0 0 4 [0 0 0 1 0 0 0 0]
[ −6 −3 0 3 6 0 ]

20 10 −12 −4 0 0 −8 −6
10 35 −6 −32 0 0 −4 −3
−12 −6 12 0 0 0 0 6
625 −4 −32 0 32 0 0 4 0
𝐑′𝟏 = 𝑥
6 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0
−8 −4 0 4 0 0 8 0
[ −6 −3 6 0 0 0 0 3]

Elément triangulaire 2 :
On a : R’2 = tT2R2T2

0 0 0 0 1 0 0 0
0 0 0 0 0 1 0 0
0 0 1 0 0 0 0 0
𝑇2 =
0 0 0 1 0 0 0 0
0 0 0 0 0 0 1 0
[0 0 0 0 0 0 0 1]

P a g e 11 | 39

EPAC-GC4-2022-2023 GROUPE 1
 Modélisation des systèmes et simulation

0 0 0 0 0 0
0 0 0 0 0 0 20 10 −8 −6 −12 −4
0 0 1 0 0 0 10 35 −4 −3 −6 −32
625 0 0 0 1 0 0 −8 −4 8 0 0 4
𝑅′ 2 =
6 1 0 0 0 0 0 −6 −3 0 3 6 0
0 1 0 0 0 0 −12 −6 0 6 16 0
0 0 0 0 1 0 [ −4 −32 4 0 0 32 ]
[0 0 0 0 0 1 ]
0 0 0 0 1 0 0 0
0 0 0 0 0 1 0 0
0 0 1 0 0 0 0 0
0 0 0 1 0 0 0 0
0 0 0 0 0 0 1 0
[0 0 0 0 0 0 0 1]

0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0
0 0 8 0 −8 −4 0 −4
625 0 0 0 3 −6 −3 6 0
𝑅′2 =
6 0 0 −8 −6 20 10 −12 −4
0 0 −4 −3 10 35 −6 −32
0 0 0 6 −12 −6 12 0
[0 0 4 0 −4 −32 0 32 ]

VIII. La matrice de rigidité globale

La matrice de rigidité globale se trouve en faisant la somme des toutes les matrice
de rigidité élémentaire globale : 𝑅 = ∑100 ′
𝑖=1 𝑅𝑒 nous sommes servi d’un
programme matlab pour calculer cette matrice.
IX. La matrice des forces F

X. La matrice des déplacements F

Les déplacements se déterminent en résolvant l’équation : RU= F

P a g e 12 | 39

EPAC-GC4-2022-2023 GROUPE 1
 Modélisation des systèmes et simulation

XI. Les contraintes et déformations

Les contraintes quant à elles se déterminent comme suit : 𝜎 = 𝜀 ∗ 𝐸

P a g e 13 | 39

EPAC-GC4-2022-2023 GROUPE 1
 Modélisation des systèmes et simulation

Conclusion
La MEF présente de nombreux avantages, tels que sa flexibilité pour traiter des
géométries complexes, sa capacité à gérer des matériaux non homogènes et
anisotropes, sa précision adaptable en fonction de la taille de la maille, et sa
capacité à résoudre des problèmes non linéaires. Toutefois, la MEF peut
également présenter des inconvénients, tels que la nécessité de définir une maille
fine pour obtenir une solution précise, des erreurs d'arrondi numériques, et des
problèmes liés aux conditions aux limites.
La méthode des éléments finis est une méthode numérique puissante pour
résoudre une grande variété de problèmes d'ingénierie et de physique. Elle est
largement utilisée dans l'industrie et la recherche, et constitue une compétence
importante pour les ingénieurs et les scientifiques dans de nombreux domaines.

P a g e 14 | 39

EPAC-GC4-2022-2023 GROUPE 1
 Modélisation des systèmes et simulation

Annexes

P a g e 15 | 39

EPAC-GC4-2022-2023 GROUPE 1
 Modélisation des systèmes et simulation

ANNEXE 1 : programme Matlab

% Propriétés de l'élément
E=2000; % module de young E
v=0.2;% coefficient de poisson v

% les coordonnées à chaque noeud de l'élément triangulaire 1

x=[0;2;0];
y=[0;0;1];

%la matrice de rigidité D pour une déformation plane

D=(E/(1-2*v))*[1 v/(1-v) 0;v/(1-v) 1 0;0 0 (1-2*v)/(2*(1-v))];

% l'aire du triangle élémentaire

A1=((x(2)-x(1))*(y(3)-y(1))-((x(3)-x(1))*(y(2)-y(1))))/2;
% matrice des déformations
B1=(1/(2*A1))*[y(2)-y(3),0,y(3)-y(1),0,y(1)-y(2),0;0,x(3)-x(2),0,x(1)-
x(3),0,x(2)-x(1);x(3)-x(2),y(2)-y(3),x(1)-x(3),y(3)-y(1),x(2)-x(1),y(1)-
y(2)];

% matrice de rigidité élémentaire dans le repère local

Re1=B1'*D*B1*A1;

% les coordonnées à chaque noeud de l'élément triangulaire 2

x=[2;0;2];
y=[1;1;0];

%la matrice de rigidité D pour une déformation plane

D=(E/(1-2*v))*[1 v/(1-v) 0;v/(1-v) 1 0;0 0 (1-2*v)/(2*(1-v))];

% l'aire du triangle élémentaire

A2=((x(2)-x(1))*(y(3)-y(1))-((x(3)-x(1))*(y(2)-y(1))))/2;
% matrice des déformations
B2=(1/(2*A2))*[y(2)-y(3),0,y(3)-y(1),0,y(1)-y(2),0;0,x(3)-x(2),0,x(1)-
x(3),0,x(2)-x(1);x(3)-x(2),y(2)-y(3),x(1)-x(3),y(3)-y(1),x(2)-x(1),y(1)-
y(2)];

% matrice de rigidité élémentaire dans le repère local

Re2=B2'*D*B2*A2;

%après execution, on note que Re1=Re2. nous utiliserons Re1 pour la suite.

%nombre de noeuds après maillage

nombre_noeuds=66;

%nombre d'élément après maillage

nombre_elements=100;

P a g e 16 | 39

EPAC-GC4-2022-2023 GROUPE 1
 Modélisation des systèmes et simulation

% matrice renseignant sur la position des noeuds des éléments dans le

% repere global
X=[1 2 12;13 12 2;2 3 13;14 13 3;3 4 14;15 14 4;4 5 15;16 15 5;5 6 16;17 16
6;6 7 17;18 17 7;7 8 18;19 18 8;8 9 19;20 19 9;9 10 20;21 20 10;10 11 21;12
13 22;23 22 13;13 14 23;24 23 14;14 15 24;25 24 15;15 16 25;26 25 16;16 17
26;27 26 17;17 18 27;28 27 18;18 19 28;29 28 19;19 20 29;30 29 20;20 21
30;22 23 31;32 31 23;23 24 32;33 32 24;24 25 33;34 33 25;25 26 34;35 34
26;26 27 35;36 35 27;27 28 36;37 36 28;28 29 37;38 37 29;29 30 38;31 32
39;40 39 32;32 33 40;41 40 33;33 34 41;42 41 34;34 35 42;43 42 35;35 36
43;44 43 36;36 37 44;45 44 37;37 38 45;39 40 46;47 46 40;40 41 47;48 47
41;41 42 48;49 48 42;42 43 49;50 49 43;43 44 50;51 50 44;44 45 51;46 47 52;
53 52 47; 47 48 53;54 53 48;48 49 54;55 54 49;49 50 55;56 55 50;50 51 56;52
53 57;58 57 53;53 54 58;59 58 54;54 55 59;60 59 55;55 56 60;57 58 61;62 61
58;58 59 62;63 62 59;59 60 63;61 62 64;65 64 62;62 63 65;64 65 66];

%initialisation de la matrice R_global

R_global=zeros(2*nombre_noeuds,2*nombre_noeuds);

%initialisation du vecteur force associé au poids volumique

FE=zeros(2*nombre_noeuds,1);

%force volumique élémentaire

rho=0.025; %en MN/m^3
f=(1/3)*[0;-rho;0;-rho;0;-rho];

%initialisation de la matrice de localisation T des éléments

T=zeros(6,2*nombre_noeuds);

for i=1:100
T=zeros(6,2*nombre_noeuds);
for j=1:3
p=X(i,j);
T(2*j-1:2*j,2*p-1:2*p)=[1 0;0 1];
end
R_global=R_global+transpose(T)*Re1*T; %matrice de rigidité globale
FE=FE+transpose(T)*f;%force volumique globale
end

%définition de la charge P
P=0.0009; %MN/m^2

%indices des noeuds où la charge est appliquée

Y=[1;12;22;31;39;46;52;57;61;64;66];

%initialisation du vecteur de charge

FR=zeros(2*nombre_noeuds,1);

%application de la charge sur les noeuds sélectionnés

for j=1:length(Y)
i=Y(j);
FR(2*i-1)=0.5*P;
end

f_global=FE+FR; %force globale appliquée au barrage

%initialisation du vecteurs contenant les déplacements nodaux

U=zeros(2*nombre_noeuds,1);

%détermination des déplacements nodaux

P a g e 17 | 39

EPAC-GC4-2022-2023 GROUPE 1
 Modélisation des systèmes et simulation

H=R_global(23:132,23:132);
U(23:132,1)=H\f_global(23:132,1);

%initialisation de la matrice de deplacement par élément

Ue=zeros(6,1);
Tabl=zeros(6,nombre_elements);

for i=1:100
%remplissage de la matrice Ue
a=X(i,1);
b=X(i,2);
c=X(i,3);
Ue(1,1)=U(2*a-1,1);
Ue(2,1)=U(2*a,1);
Ue(3,1)=U(2*b-1,1);
Ue(4,1)=U(2*b,1);
Ue(5,1)=U(2*c-1,1);
Ue(6,1)=U(2*c,1);
Tabl(:,i)=Ue;
end

%pour l'ensemble des triangles de type 1

tab1=Tabl(:,[1:2:19,20:2:36,37:2:51,52:2:64,65:2:75,76:2:84,85:2:91,92:2:96
,97:2:99,100]);
W=zeros(3,55); Z=zeros(3,55);
for j=1:55
epsilon_j=B1*tab1(:,j);
W(:,j)=epsilon_j;
sigma_j=D*epsilon_j;
Z(:,j)=sigma_j;
end
disp(W); %les deplacements
disp(Z); %les contraintes

%pour l'ensemble des triangles de type 2

tab2=Tabl(:,[2:2:18,21:2:35,38:2:50,53:2:63,66:2:74,77:2:83,86:2:90,93:2:95
,98]);
w=zeros(3,45); z=zeros(3,45);
for k=1:45
epsilon_k=B2*tab2(:,k);
w(:,k)=epsilon_k;
sigma_k=D*epsilon_k;
z(:,k)=sigma_k;
end
disp(w); %les deplacements
disp(z); %les contraintes

P a g e 18 | 39

EPAC-GC4-2022-2023 GROUPE 1
 Modélisation des systèmes et simulation

Annexe 2 : matrice de force

0.0004

-0.0083

-0.0250

-0.0250

-0.0250

-0.0250

-0.0250

-0.0250

-0.0250

-0.0250

-0.0250

-0.0083

0.0004

-0.0250

-0.0500

P a g e 19 | 39

EPAC-GC4-2022-2023 GROUPE 1
 Modélisation des systèmes et simulation

-0.0500

-0.0500

-0.0500

-0.0500

-0.0500

-0.0500

-0.0500

-0.0250

0.0004

-0.0250

-0.0500

-0.0500

-0.0500

-0.0500

-0.0500

0
P a g e 20 | 39

EPAC-GC4-2022-2023 GROUPE 1
 Modélisation des systèmes et simulation

-0.0500

-0.0500

-0.0250

0.0004

-0.0250

-0.0500

-0.0500

-0.0500

-0.0500

-0.0500

-0.0500

-0.0250

0.0004

-0.0250

-0.0500

-0.0500

0
P a g e 21 | 39

EPAC-GC4-2022-2023 GROUPE 1
 Modélisation des systèmes et simulation

-0.0500

-0.0500

-0.0500

-0.0250

0.0004

-0.0250

-0.0500

-0.0500

-0.0500

-0.0500

-0.0250

0.0004

-0.0250

-0.0500

-0.0500

-0.0500

0
P a g e 22 | 39

EPAC-GC4-2022-2023 GROUPE 1
 Modélisation des systèmes et simulation

-0.0250

0.0004

-0.0250

-0.0500

-0.0500

-0.0250

0.0004

-0.0250

-0.0500

-0.0250

0.0004

-0.0250

-0.0250

0.0004

-0.0083

P a g e 23 | 39

EPAC-GC4-2022-2023 GROUPE 1
 Modélisation des systèmes et simulation

Annexe 3 : matrice des déplacements

1.0e-03 *

-0.0307

-0.0737

-0.0130

P a g e 24 | 39

EPAC-GC4-2022-2023 GROUPE 1
 Modélisation des systèmes et simulation

-0.0635

-0.0044

-0.0560

0.0002

-0.0485

0.0025

-0.0410

0.0036

-0.0335

0.0039

-0.0263

0.0036

-0.0191

0.0026

-0.0118

0.0007

-0.0041

-0.0478

-0.1383

-0.0258

-0.1213

-0.0125

-0.1059

-0.0047

-0.0904

-0.0005

-0.0749

0.0015
P a g e 25 | 39

EPAC-GC4-2022-2023 GROUPE 1
 Modélisation des systèmes et simulation

-0.0598

0.0021

-0.0450

0.0015

-0.0305

-0.0004

-0.0157

-0.0601

-0.1945

-0.0383

-0.1721

-0.0234

-0.1490

-0.0140

-0.1254

-0.0086

-0.1017

-0.0060

-0.0786

-0.0052

-0.0563

-0.0059

-0.0343

-0.0713

-0.2428

-0.0516

-0.2155

-0.0371
P a g e 26 | 39

EPAC-GC4-2022-2023 GROUPE 1
 Modélisation des systèmes et simulation

-0.1850

-0.0275

-0.1531

-0.0218

-0.1210

-0.0188

-0.0899

-0.0176

-0.0600

-0.0836

-0.2838

-0.0668

-0.2511

-0.0540

-0.2133

-0.0453

-0.1732

-0.0399

-0.1328

-0.0368

-0.0937

-0.0982

-0.3174

-0.0848

-0.2789

-0.0746

-0.2338

-0.0677
P a g e 27 | 39

EPAC-GC4-2022-2023 GROUPE 1
 Modélisation des systèmes et simulation

-0.1854

-0.0634

-0.1368

-0.1161

-0.3436

-0.1063

-0.2988

-0.0994

-0.2462

-0.0953

-0.1896

-0.1375

-0.3622

-0.1318

-0.3108

-0.1291

-0.2505

-0.1628

-0.3729

-0.1613

-0.3148

-0.1915

-0.3756

P a g e 28 | 39

EPAC-GC4-2022-2023 GROUPE 1
 Modélisation des systèmes et simulation

Annexe 4 : matrice des déformations

- Déformation pour les triangles élémentaires de type 1

P a g e 1 | 39

EPAC-GC4-2022-2023 GROUPE 1
 Modélisation des systèmes et simulation

- déformations pour les triangles élémentaires de type 2

P a g e 2 | 39

EPAC-GC4-2022-2023 GROUPE 1
 Modélisation des systèmes et simulation

P a g e 3 | 39

EPAC-GC4-2022-2023 GROUPE 1
 Modélisation des systèmes et simulation

P a g e 4 | 39

EPAC-GC4-2022-2023 GROUPE 1
 Modélisation des systèmes et simulation

Annexe 5 : matrice des contraintes

- Contraintes des triangles élémentaires de type 1

-0.0614 -0.0529 -0.0466 -0.0404 -0.0342 -0.0280 -0.0219 -0.0159 -0.0099 -0.0034 -0.0243

-0.2455 -0.2117 -0.1866 -0.1618 -0.1367 -0.1118 -0.0876 -0.0637 -0.0394 -0.0138 -0.2081

-0.0384 -0.0162 -0.0055 0.0002 0.0032 0.0045 0.0049 0.0045 0.0033 0.0009 -0.0150

-0.0339 -0.0340 -0.0309 -0.0265 -0.0214 -0.0162 -0.0111 -0.0065 -0.0101 -0.0202 -0.0231

-0.1890 -0.1645 -0.1387 -0.1127 -0.0873 -0.0626 -0.0385 -0.0138 -0.1780 -0.1639 -0.1407

-0.0113 -0.0055 -0.0014 0.0009 0.0019 0.0022 0.0019 0.0010 -0.0047 -0.0060 -0.0040

-0.0220 -0.0189 -0.0148 -0.0103 -0.0064 -0.0040 -0.0112 -0.0143 -0.0141 -0.0117 -0.0081

-0.1148 -0.0883 -0.0625 -0.0377 -0.0134 -0.1521 -0.1382 -0.1158 -0.0900 -0.0634 -0.0374

-0.0019 -0.0007 -0.0002 -0.0000 -0.0000 -0.0001 -0.0022 -0.0023 -0.0020 -0.0020 -0.0020

-0.0044 -0.0013 -0.0054 -0.0076 -0.0072 -0.0048 -0.0013 -0.0000 -0.0018 -0.0026 -0.0013

-0.0127 -0.1283 -0.1126 -0.0904 -0.0646 -0.0379 -0.0122 -0.1051 -0.0873 -0.0647 -0.0387

-0.0019 0.0018 0.0001 -0.0012 -0.0022 -0.0032 -0.0038 0.0022 0.0011 -0.0007 -0.0028

P a g e 5 | 39

EPAC-GC4-2022-2023 GROUPE 1
 Modélisation des systèmes et simulation

0.0019 0.0006 0.0005 0.0010 0.0037 0.0007 0.0015 0.0033 0.0006 0.0012 0.0003

-0.0120 -0.0818 -0.0621 -0.0388 -0.0122 -0.0579 -0.0371 -0.0126 -0.0334 -0.0122 -0.0083

-0.0049 0.0018 0.0013 -0.0009 -0.0041 0.0012 0.0010 -0.0017 0.0006 0.0008 0.0005

- Contraintes des triangles élémentaires de type 2

-0.0233 -0.0324 -0.0329 -0.0302 -0.0261 -0.0214 -0.0164 -0.0114 -0.0067 -0.0115 -0.0194

-0.2043 -0.1830 -0.1599 -0.1357 -0.1114 -0.0875 -0.0638 -0.0398 -0.0146 -0.1834 -0.1608

-0.0099 -0.0008 0.0048 0.0079 0.0092 0.0094 0.0090 0.0078 0.0057 -0.0054 -0.0005

-0.0220 -0.0212 -0.0185 -0.0147 -0.0105 -0.0065 -0.0061 -0.0111 -0.0135 -0.0133 -0.0112

-0.1365 -0.1114 -0.0865 -0.0623 -0.0383 -0.0138 -0.1604 -0.1377 -0.1126 -0.0869 -0.0616

0.0036 0.0059 0.0069 0.0070 0.0065 0.0054 -0.0017 0.0008 0.0032 0.0046 0.0050

-0.0080 -0.0044 -0.0033 -0.0057 -0.0071 -0.0066 -0.0044 -0.0012 -0.0017 -0.0022 -0.0022

-0.0371 -0.0129 -0.1362 -0.1137 -0.0882 -0.0621 -0.0365 -0.0119 -0.1117 -0.0890 -0.0634

0.0049 0.0044 0.0005 0.0019 0.0031 0.0036 0.0035 0.0031 0.0015 0.0025 0.0029

-0.0009 0.0020 -0.0008 -0.0000 0.0012 0.0039 -0.0004 0.0011 0.0033 -0.0004 0.0010

-0.0369 -0.0113 -0.0870 -0.0640 -0.0381 -0.0115 -0.0623 -0.0387 -0.0123 -0.0376 -0.0131

0.0025 0.0019 0.0016 0.0024 0.0022 0.0011 0.0011 0.0018 0.0008 0.0003 0.0006

P a g e 6 | 39

EPAC-GC4-2022-2023 GROUPE 1
 Modélisation des systèmes et simulation

-0.0008

-0.0127

-0.0005

P a g e 7 | 39

EPAC-GC4-2022-2023 GROUPE 1