Université Abou Bekr-Belkaid Tlemcen

Département Génie Civile

Option : Structure

Thème :
Rapport sur un Algorithme d’un Bâtiment

Nom / Prénom :
- Belkhiri Ismail Abdelbasset
- Hamdache Abdelkader
 Méthodes des Eléments finis
Portique dans l’espace
1-Discrétisation
On a une structure (portique) 3D de 3travée selon « x » et 5 travée
selon « y » et 2 selon « z » avec la longueur de travée selon « x »
poutre principale est 4m et selon « z » poutre secondaires de 4m
et selon « y » les poteaux 3m.

Noeuds X Y Z
1 0 0 0
2 4 0 0
3 8 0 0
4 12 0 0
5 0 3 0
6 4 3 0
7 8 3 0
8 12 3 0
9 0 6 0
10 4 6 0
11 8 6 0
12 12 6 0
13 0 9 0
14 4 9 0
15 8 9 0
16 12 9 0
17 0 12 0
18 4 12 0
19 8 12 0
20 12 12 0
21 0 15 0
22 4 15 0
23 8 15 0
24 12 15 0
25 0 0 4
26 4 0 4
27 8 0 4
28 12 0 4
29 0 3 4
30 4 3 4
31 8 3 4
32 12 3 4
33 0 6 4
34 4 6 4
35 8 6 4
36 12 6 4
37 0 9 4
38 4 9 4
39 8 9 4
40 12 9 4
41 0 12 4
42 4 12 4
43 8 12 4
44 12 12 4
45 0 15 4
46 4 15 4
47 8 15 4
48 12 15 4
49 0 0 8
50 4 0 8
51 8 0 8
52 12 0 8
53 0 3 8
54 4 3 8
55 8 3 8
56 12 3 8
57 0 6 8
58 4 6 8
59 8 6 8
60 12 6 8
61 0 9 8
62 4 9 8
63 8 9 8
64 12 9 8
65 0 12 8
66 4 12 8
67 8 12 8
68 12 12 8
69 0 15 8
70 4 15 8
71 8 15 8
72 12 15 8

Les éléments :
L’élements Origine Extrémité Type
1 1 25 1
2 2 26 1
3 3 28 1
4 4 29 1
5 5 30 1
6 6 31 1
7 7 32 1
8 8 33 1
9 9 34 1
10 10 35 1
11 11 36 1
12 12 37 1
13 13 38 1
14 14 39 1
15 15 40 1
16 16 41 1
17 17 42 1
18 18 43 1
19 19 44 1
20 20 45 1
21 21 46 1
22 22 47 1
23 23 48 1
24 24 49 1
25 25 50 1
26 26 51 1
28 28 52 1
29 29 53 1
30 30 54 1
31 31 55 1
32 32 56 1
33 33 57 1
34 34 58 1
35 35 59 1
36 36 60 1
37 37 61 1
38 38 62 1
39 39 63 1
40 40 64 1
41 41 65 1
42 42 66 1
43 43 67 1
44 44 68 1
45 45 69 1
46 46 70 1
47 47 71 1
48 48 72 1
49 25 26 2
50 26 27 2
51 27 28 2
52 29 30 2
53 30 31 2
54 31 32 2
55 33 34 2
56 34 35 2
57 35 36 2
58 37 38 2
59 38 39 2
60 39 40 2
61 41 42 2
62 42 43 2
63 43 44 2
64 45 46 2
65 46 47 2
66 47 48 2
67 49 50 2
68 50 51 2
69 51 52 2
70 53 54 2
71 54 55 2
72 55 56 2
73 57 58 2
74 58 59 2
75 59 60 2
76 61 62 2
77 62 63 2
78 63 64 2
79 65 66 2
80 66 67 2
81 67 68 2
82 69 70 2
83 70 71 2
84 71 72 2
85 25 29 2
86 29 33 2
87 33 37 3
88 37 41 3
89 41 45 3
90 26 30 3
91 30 34 3
92 34 38 3
93 38 42 3
94 42 46 3
95 27 31 3
96 31 35 3
97 35 39 3
98 39 43 3
99 43 47 3
100 28 32 3
101 32 36 3
102 36 40 3
103 40 44 3
104 44 48 3
105 49 53 3
106 53 57 3
107 57 61 3
108 61 65 3
109 65 69 3
110 50 54 3
111 54 58 3
112 58 62 3
113 62 66 3
114 66 70 3
115 51 55 3
116 55 59 3
117 59 63 3
118 63 67 3
119 67 71 3
120 52 56 3
121 56 60 3
122 60 64 3
123 64 68 3
124 68 72 3
125 73 76 3
126 77 81 3
127 81 85 3
 128 85 89 3
129 89 93 3
130 93 97 3
131 78 82 3
132 82 86 3
133 86 90 3
134 90 94 3
135 94 98 3
136 83 83 3
137 83 87 3
138 87 91 3
139 91 95 3
140 95 99 3
141 80 84 3
142 84 88 3
143 88 92 3
144 92 96 3
145 96 100 3

2-Matrice de rigidité locale

L’élément dans cette structure dans l’espace à deux nœuds et
donc la matrice de rigidité locale élémentaire de dimension
(6*02)=12 soit (12*12) elle est obtenue par superposition des
types de déplacements qui composent l’élément soit les
déplacements normaux, la torsion autour de ox, la flexion dans le
plan oxy et la flexion dans le plan oxz, Les matrices ainsi utilisées
dans la superposition

3-Transformation :
Le passage du système local au système global s’appelle la
transformation, en usant des vecteurs unitaires des 6 axes (3
locaux ox, oy et oz et 3 globaux OX, OY,et OZ)
4-Matrice de rigidité élémentaire globale :
Elle est toujours obtenue par l’équation [ K ( e ) ]=¿
¿La matrice de transformation
[ K ( e ) ]La matrice de rigidité locale
On obtient la matrice [ K ( e ) ] de dimension 12*12 de rigidité
élémentaire globale OXYZ

L’algorithme de la méthode des éléments pour ce système sur

MATLAB :