REPUBLIQUE DU BENIN

¤¤¤¤¤¤
UNIVERSITE D’ABOMEY CALAVI (UAC)
¤¤¤¤¤
ECOLE POLYTECHNIQUE D’ABOMEY
CALAVI (EPAC)
¤¤¤¤¤
SECTEUR INDUSTRIEL
¤¤¤¤¤¤
DEPARTEMENT DU GENIE CIVIL
¤¤¤¤¤

METHODE DES ELEMENTS FINIS (MEF)

APPLICATION
CHAPITRE I : RESOLUTION D’UN PROBLEME D’ELASTCITE

Document rédigé par :

Génie Civil 4ème Année (GC4)

Sous la direction de :
Prof Rostand MOUTOU PITTI

Année Universitaire 2023-2024

 APPLICAITON CHAPITRE 1 : RESOLUTION D’UN PROBLEME D’ELASTCITE

5-1-Solution de la théorie des poutre.

Expression de la flèche v ( x ) et du moment fléchissant M ( x) .

• Degré d’hyperstaticité
D = Co − 3 avec Co le nombre d’inconnues et 3 représente le nombre d’équations
de la statique.
Or Co = 3 donc D = 3 − 3 d’où D = 0
Par conséquent la poutre est isostatique.
• Les réactions aux appuis
À l’équilibre statique, on a :
 F =0
 F = 0  x

   Fx = 0
 M F /  = 0 
 M  = 0

F x = 0  H1 = 0
F y = 0  R1 + R2 = −Q

L
M 1 = 0  Q 2 + R L = 0 2

Q
 R2 = −
2
Q qL
Or la poutre étant symétriquement chargée, R1 = R2 = − =−
2 2

GC4 EPAC 1
 APPLICAITON CHAPITRE 1 : RESOLUTION D’UN PROBLEME D’ELASTCITE

Détermination du moment fléchissant

 x
 M o = 0  − R x − qx  2  − M ( x) = 0
1

L x2
 M ( x) = q x−q
2 2
x
M ( x) = q ( L − x)
2

Détermination de l’effort tranchant

F y = 0  R1 + Q − T ( x) = 0
 T ( x) = Q + R1

L
T ( x) = qx − q
2

Détermination de la flèche
EIv( x) = − M ( x)
x2 L
EIv( x) = q −q x
2 2
 x2 L 
  EIv( x) =   q − q x  dx
 2 2 
x3 L
 EIv( x) = Co + q − q x2
6 4
 x 3
L 
  EIv( x) =   Co + q − q x 2 dx
 6 4 
x4 L
 EIv( x) = C1 + C0 x + q − q x3
24 12

Conditions aux limites

v(0) = 0  C1 = 0
 x = 0; v( x) = 0 
À ; L4 qL 3
 x = L; v( x) = 0 v( L) = 0  C0 L + q −  L
 24 12

C1 = 0 C0 = 0
 
 qL3 qL3   qL3
 0
C = −  1
C =
 12 24  24

1  x2 x3 L3 
v( x) =  q − qL + q x
EI  24 12 24 

GC4 EPAC 2
 APPLICAITON CHAPITRE 1 : RESOLUTION D’UN PROBLEME D’ELASTCITE

v( x) =
q
24 EI
( x 4 − 2 Lx3 + L3 x )
5-2- Application du PTV
2-1- Déterminons les expressions des TV des efforts extérieurs et intérieurs
en fonction du déplacement virtuel

Formulation variationnelle (dérivée du PTV)

 wi +  we = 0  u CA

➢ Déterminons les travaux intérieurs

− wi =  J =   ij ij dv


 J =   dv or  x =  x E (Loi de Hooke pour un matériau homogène, isotrope et

v

élastique)
 J =   x E x dv (1)
v

 2v
v( x) = (a)
x 2
M
v( x ) = − (b)
EIz
M
= y (c) (Flexion simple)
Iz

 2v M
(a)=(b)  2
=− (d)
dx EIz
M 
De ©, on a : = (e)
Iz y

 2v  
(e) dans (d) donne : 2
=− or x =  x (d’après la loi de Hooke) 
dx Ey E

 2v x  2v
= −   = − y
x 2 x 2
x
y

 2 v
De même,  x = − y
x 2

GC4 EPAC 3
 APPLICAITON CHAPITRE 1 : RESOLUTION D’UN PROBLEME D’ELASTCITE

 2v  ² v 
L’équation (1) devient :  J =  − y  E   − y dV
x ²  x ² 

 2 v  ² v
 J =  Ey ² dV or  dV =  [Link] or x  0; L
x ² x ²

Donc
 2 v  ² v
 J = E  y ² dS  dx or  y²dS = Iz donc
x ² x ²

L  2 v  ² v
 J = EIz  dx
0 x ² x ²

➢ Déterminons les travaux extérieurs

 we =   ui dV + 
V S
f i 
f i
 ui dS  u CA

Hypothèse : Les forces volumiques sont négligeables.

Donc  we = 
S
f  ui dS
i

Avec dS = dxdz
ui étant le déplacement transversal des points de la poutre. On a u = v ( x )

 we =  
S
f  v( x)dxdz or f dz = f s z = q 
s


L L
 we =  q v( x)dx   we = q   v( x )dx
0 0

Donc  wi +  we = 0   we = − wi =  J
L L  2 v  ² v
 q   v( x)dx = EIz  dx  u CA (2)
0 0 x ² x ²

2-2-a- Expression approchée de la flèche sous la forme v( x) =  x ² +  x + 

Conditions aux limites

 x = 0; v( x) = 0 v(0) = 0   = 0  = 0
À  
 x = L; v( x) = 0 v( L) = 0   L² +  L = 0   = − L

On a donc
v( x) =  x ² −  Lx  v( x) =  ( x ² − Lx) (3)

GC4 EPAC 4
 APPLICAITON CHAPITRE 1 : RESOLUTION D’UN PROBLEME D’ELASTCITE

(3) dans (2) donne

² ²
q   ( x ² − Lx ) dx = EIz   ( x ² − Lx ) •  ( x ² − Lx ) dx
L L

0 0 x ² x ²
² ²
Or  ( x ² − Lx ) = 2 et  ( x ² − Lx ) = 2
x ² x ²
L
1 1  L
 q  x 3 − Lx ²  = EIz  4 dx
3 2 0 0

1 1 
 q  L3 − L3  = EIz  4 x 
L

3 2  0

L3
 −q = 4 EIz L
6
2
L
 −q = 4 EIz
6
qL ²
D’où  = −
24 EIz
L²
α dans (3) donne : v( x) = −q ( x ² − Lx )
24 EIz
L² x
Par conséquent v( x) = −q ( x − L)
24 EIz

Calcul d’erreur Err sur Vmax

L
À x= ; v = vmax
2

• Selon la théorie des poutres on a :

q  L  
4 3
L 3 L
vRdm =   − 2 L   + L    
24 EIz  2  2  2  

q L L  q  L4 − 4 L4 + 8L4 
4 4 4
L
vRdm =  − +  =  
24 EIz  16 4 2  24 EIz  16 

5qL4
D’où vRdm =
384 EIz

• Selon PTV, on a :
L
qL ²  
L  2 L− L
vmax = v  =  
2 24 EIz  2

GC4 EPAC 5
 APPLICAITON CHAPITRE 1 : RESOLUTION D’UN PROBLEME D’ELASTCITE

L qL4
d’où v   =
 2  96 EIz

L
vRdm − v  
On a donc Err = 2
vRdm

L
vRdm − v  
Err = 2
vRdm

qL4  5 1  1
 − 
EIz  384 96  384 1
Err = = = F
5qL4 5 5
384 EIz 384

Err = 0, 2 soit Err = 20 %

Calcul d’erreur Err’ sur Mmax

• Selon la théorie des poutres, on a :
x
M ( x) = q ( L − x)
2
M max = M L ( 2 ) = q L4  L − L2  = q L8
• Selon PTV
 ²v( x)
M ( x) = − EIz
x ²
v( x) qL ²
Or = ( L − 2x)
x 24 EIz
 ²v( x) 2qL ²
D’où =−
x ² 24 EIz
qL ²
Donc M ( x) =
12 EIz

L
M Rdm − M  
Err = 2
M Rdm

GC4 EPAC 6
 APPLICAITON CHAPITRE 1 : RESOLUTION D’UN PROBLEME D’ELASTCITE

L² L² 1 1 1
q −q −
Err = 8 12  Er  = 8 12 = 24
r
L² 1 1
q
8 8 8
1
Err = 0,3333 soit Err = 33,33 %
3

Conclusion
Une erreur de 20% sur la flèche maximale (vmax) et de 33,33% sur le moment
fléchissant maximal (Mmax) or l’erreur admissible en Génie Civil est de 2%.
Par conséquent les résultats ne sont pas précis.
Pour améliorer l’approximation, on peut soit augmenter le degré de
l’interpolation, soit garder le même degré d’interpolation mais subdiviser le
domaine d’étude. Dans notre cas, nous allons subdiviser le domaine d’étude.
2-2-b-
On divise la poutre en quatre éléments de même longueur. Sur chaque él
ément i recherchons la flèche sous la forme approchée de a), c’èst-à dire
v( x) =  i x ² + i x +  i
Considérons maintenant le schéma ci-après

vi ( x) =  i x ² + i x +  i

Premier élément : x  0;  3

L
 4

v1 ( x) = 1 x ² + 1 x +  1 (*)

Deuxième élément : x   ; 
L L
4 2

v2 ( x) =  2 x ² +  2 x +  2 (**)

GC4 EPAC 7
 APPLICAITON CHAPITRE 1 : RESOLUTION D’UN PROBLEME D’ELASTCITE

Conditions aux limites

 x = 0; v1 ( x) = 0 a

 x = L ; v1 ( x) = v2 ( x)  v1  L  = v2  L  b
 4    
4 4

 L L L c
À  x = ; v1 ( x) = v2 ( x)  v1   = v2  
 4 4 4
 L d
 x = ; v2 ( x) = 0
 4
De (a), v1 (0) = 0   1 = 0
2
L L  L
v1   = 1   + 1   et
4 4 4
2
L L L
v2   =  2   +  2   +  2
4 4 4
2 2

De (b), v1   = v2    1   + 1   =  2   +  2   +  2
L L L L L L
4 4 4 4 4 4

v1 ( x) = 21 x + 1 et v2 ( x) = 2 2 x + 2 donc

L  L L  L
v1   = 21   + 1 et v2   = 2 2   + 2
4 4 4 4

De (c), v1   = v2    21   + 1 = 2 2   + 2

L L L L
4 4 4 4
L
 (1 −  2 ) =  2 − 1
2

Aussi De (d),
L  L L L
v2   = 2 2   + 2 or v2   = 0  2 2   +  2 = 0
2 2 2 2

  2 = − L 2

GC4 EPAC 8
 APPLICAITON CHAPITRE 1 : RESOLUTION D’UN PROBLEME D’ELASTCITE

 L2 L
16 (1 −  2 ) + 4 ( 1 −  2 ) =  2 (a)

Soit  (1 −  2 ) =  2 − 1 (b)
L
2 (c )
  2 = − 2 L



L2 L L
(b) dans (a) donne (1 −  2 ) −  (1 −  2 ) =  2
16 4 2
2
L2 2L
(1 −  2 ) − (1 −  2 ) =  2
16 4

L2
Soit − (1 −  2 ) =  2 (d)
16
L L
(s) dans (b) donne : (1 −  2 ) = − 2 L − 1  1 = − (1 +  2 ) (e)
2 2
L
(e) dans (*) donne : v1 ( x) = 1 x ² − ( 2 − 1 ) x
2
2
L
(c) dans (**) donne : v2 ( x) =  2 x² −  2 Lx − (1 −  2 )
16

Selon PTV
D’après la question 5-5-a , on a :
L  2 v( x)  ² v( x)
 J = EIz  dx  u CA
0 x ² x ²
L
 we = q   v( x )dx
0

Premier élément : x  0; 

L
 4
L
 2 v1 ( x)  ² v1 ( x)
 J1 = EIz  4 dx
0 x ² x ²
L
Or v1 ( x) = 21 x − (1 −  2 )
2

v1 ( x) = 21

GC4 EPAC 9
 APPLICAITON CHAPITRE 1 : RESOLUTION D’UN PROBLEME D’ELASTCITE

L
Donc  J1 = EIz  21   (21 )dx
0
4

L
= 4 EIz11  x 04

 J1 = EIz11 L

De même,
L
 we1 = q  4  v1 ( x)dx
0

 
L
L
= q  4  1 x 2 − (1 +  2 ) x  dx
0
 2 
L
1 L 4
= q  1 x3 − (1 +  2 ) x 2 
3 4 0

 1  L 3 L L 
2

= q  1   − (1 +  2 )   
 3  4  4  4  

 L
3
L 
3
= q 1 −  (1 +  2 ) 
 192 64 

 L 3 3L
3

= q 1 − (1 +  2 )
192 192 

 2L 3 3L
3

= q − 1 −  2 
 192 192 
3
L
 we1 = −q  21 + 3 2 
192

Deuxième élément : x   ; 
L L
4 2
L
 2v2 ( x)  ² v2 ( x)
 J 2 = EIz L2 dx
4 x² x²

Or v2 ( x) = 2 2 x −  2 L

Donc v2 ( x) = 2 2

GC4 EPAC 10
 APPLICAITON CHAPITRE 1 : RESOLUTION D’UN PROBLEME D’ELASTCITE

L
 J 2 = EIz  2 2   (2 2 )dx
2
L
4

L
= EIz L2 4 2   2 dx
4

L
= EIz  4 2   2 x L2
4

L L
= 4 EIz 2   2  − 
2 4

Soit  J 2 = EIz 2   2 L

De même,
L
 we 2 = q   v2 ( x)dx
2
L
4

L
1 L L
2
2
= q   2 x3 −  2 x 2 − (1 −  2 ) x 
3 2 16 L
4

1  L  3  L  3  L  L  2  L  2  L 2 L L 
= q   2   −    −  2   −    − (1 −  2 )  −  
 3  2   4   2  2   4   16  2 4  
 7 L3 3L3 L3 
= q  2 −  2 − (1 −  2 ) 
192 32 64 

 8L3 L3  qL3
= q −  2 − 1  d’où  we 2 = − (8 2 + 31 )
 192 64  192

Or
 J =  J1 +  J 2

= EIz11 L + EIz 2 2 L

 J = EIzL(11 +  2 2 )

De même ,
 we =  we1 +  we 2
3
L qL3
 we = −q  21 + 3 2  − (8 2 + 31 )
192 192

GC4 EPAC 11
 APPLICAITON CHAPITRE 1 : RESOLUTION D’UN PROBLEME D’ELASTCITE

3
L
 we = −q 51 + 11 2  ∀ 1
192

D’après PTV, on a :
 wi +  we = 0   we = − wi =  J

L 3
 −q 5 + 11  = EIzL(  +   )
192  1 2 1 1 2 2

Par identification, on a :
5qL ²
1 = −
192 EIz
5qL ²
 2 = − 2 L alors  2 = L
192 EIz

(1 +  2 ) donc 1 = −  −
L L 5qL² 11qL² 
1 = − − 
2 2  192 EIz 192 EIz 

16qL ² qL ²
1 = − par conséquent, 1 = −
384 EIz 24 EIz
L²
2 = ( 2 − 1 )
16

L²  11qL² 5qL² 
2 = − + 
16  192 EIz 192 EIz 

−6qL4 qL4
2 = Soit  2 = −
16 192  EIz 512 EIz
1 = 0

Alors on a :
5qL ² qL3 qL ²
v1 ( x) = − x² + x  v1 ( x) = ( −5 x ² + 8 x )
192 EIz 24 EIz 192 EIz

11qL ² 11qL3 qL4

Et v2 ( x) = − x² + x−
192 EIz 192 EIz 512 EIz
qL ²
v2 ( x) = − ( −88 x ² + 88 Lx − 3L ²)
1536 EIz

GC4 EPAC 12
 APPLICAITON CHAPITRE 1 : RESOLUTION D’UN PROBLEME D’ELASTCITE

Calcul d’erreur Err sur vmax

L
La flèche est maximale à x =
2

L L L
Or  ; : Deuxième élément
2  4 2 
2
L qL² L L
v2   = − (−88   + 88L   − 3L²)
2 1536 EIz 2 2

 L  L  19qL
4
qL²
v2   = (19L²)  v2   =
 2  1536EIz  2  1536 EIz

L 5qL4 19qL4

vRdm − v   −
Err =  2  = 384 EIz 1536 EIz
vRdm 5qL4
384 EIz
1
Err = = 0, 05 Soit Err =5%
20

Calcul d’erreur Err’ dur Mmax

• Selon la théorie des poutres, on a :
x
M ( x) = q ( L − x)
2

L L L
M max = M   = q  L − 
2 4 2

L²
M max = q
8

• Selon PTV
 ²v2 ( x)
M ( x) = − EIz
x ²
 ² v2 ( x) 176qL ²
=− donc
x ² 1536 EIz

qL ² Alors M  L  = 11 qL²
11
M ( x) =
96  2  96

GC4 EPAC 13
 APPLICAITON CHAPITRE 1 : RESOLUTION D’UN PROBLEME D’ELASTCITE

L
M Rdm − M  
Err ' = 2
M Rdm

qL ² 11 1 11 1
− qL ² −
Err ' = 8 96 = 8 96 = 96
qL ² 1 1
8 8 8
8
Err ' = = 0, 0833 Soit Err’=8,33%
96

Conclusion
On constate que les erreurs sur vmax et Mmax lorsqu’on a subdivisé la poutre en
quatre parties égales sont quatre fois plus petites que lorsqu’on a étudié la poutre
comme un seul élément.
Cela nous permet de conclure que plus on discrétise le domaine, plus les calculs
sont précis.

i x 
c- Trouvons v(x) sous la forme : v( x) = i sin 
n


i =1  L 

Pennons i=1
x  L  v ( x )  ² v ( x )
2
v( x) = 1 sin   on a :  J = EIz  dx
 L  0 x ² x ²

v( x)   
 = 1 cos  x 
x L L 

v²( x) ²  
 = −1 sin  x  D’où
x² L² L 

L  ²     ²   
 J = EIz   −1 sin  x  •   −1 sin  x    dx
0
 L² L   L²  L 

L 4   
 J = EIz   1 • 1 sin 2  x   dx
 L 
4
0
L

4 L   1 − cos(2 )
 J = EIz 1 • 1  sin 2  x  dx Or sin 2  = Alors
L 
4
L 0 2

GC4 EPAC 14
 APPLICAITON CHAPITRE 1 : RESOLUTION D’UN PROBLEME D’ELASTCITE

 
1 − cos  x 
 4
L
 L  dx
 J = EIz 1 • 1 
L4 0 2
L
EIz 41 • 1  L  2  
J =  x− sin  x 
 2  L  0
4
2L

EIz 41 • 1

J = ( L − 0)
2 L4

EIz 41 • 1

J =
2 L3

De même,
L
 we = q   v( x )dx
0

L x 
 we = q  1 sin   dx
0
 L 
L
 L   x 
 we = q1  − cos   
   L 0

 L 
 we = q1  − ( cos ( ) − cos ( 0 ) )
  
L
 we = −q1  −1 − 1

L
 we = 2q1


Selon le PTV, on a :
 wi +  we = 0   we = − wi =  J

L EIz 41 • 1

 2q1 =
 2 L3

L4
 4q1 = EIz 41 • 1


4qL4
 1 =
EIz 5

GC4 EPAC 15
 APPLICAITON CHAPITRE 1 : RESOLUTION D’UN PROBLEME D’ELASTCITE

4qL4 x 
Donc v( x) = sin  
EIz 5
 L 

Calcul d’erreur Err sur vmax

 L  4qL  L 
4
v  = sin    Selon PTV
 2  EIz L 2
5

 L  4qL  
4
v  = sin  
 2  EIz 2
5

 L  4qL
4
v  =
 2  EIz
5

L 5qL4 4qL4

vRdm − v   − 5
Err =  2  = 384 EIz  EIz
vRdm 5qL4
384 EIz

5 4
− 5
Err = 384 
5
384 EIz

Err = −0, 00386 Soit Err = 0, 4 %

Calcul d’erreur Err’ dur Mmax

 ²v2 ( x)
M ( x) = − EIz Or
x ²

 ²v2 ( x) qL4  2  
= −4 5  2 sin  x 
x ²  EIz L L 

 ²v2 ( x) qL2  
= −4 3 sin  x  Donc
x ²  EIz  L 

qL2  
M ( x) = 4 sin  x  Alors
3
L 

L qL2  L
M   = 4 3 sin   
2  L 2

GC4 EPAC 16
 APPLICAITON CHAPITRE 1 : RESOLUTION D’UN PROBLEME D’ELASTCITE

L  
2
qL
M   = 4 3 sin  
2  2

L
2
qL
M =4 3 D’où
2 

L
M Rdm − M  
Err ' = 2
M Rdm

qL2 4qL2
− 3
Err ' = 8 
2
qL
8

Err ' = −0, 00801 0, 008 Soit Err ' = 0,8 %

Conclusion
Ce modèle basé sur la fonction trigonométrique est beaucoup plus précis que celui
basé sur les fonctions polynômiales

GC4 EPAC 17