Méthode des éléments finis: Corrigé du TD2

Ibrahim ALAME

14/02/2024

Problème à résoudre
Nous avons
dN
= p(x); 0≤x≤L équation d’équilibre
dx
du
σ = Eε =⇒ N = EA 0≤x≤L loi de comportement
dx
N (0) = −f1 ; N (L) = f2 conditions aux limites
2

 −EA d u = p(x) 0 ≤ x ≤ L

équiliblre + loi de comportement
=⇒ dx2
 EA du (0) = −f0 ; EA du (L) = fL

conditions aux limites de Neumann
dx dx

Formulation variationnelle
On cherche une solution u ∈ V = H 1 (Ω). Multiplions l’équation −EAu′′ = p par une fonction test v ∈ V :

d2 u
−EA ·v =p·v
Z L dx2 Z L
2
d u
− EA 2 · v(x) dx = p · v(x)dx
0 dx  0
Z L L Z L
du dv du
EA · dx − EA · v(x) = p · v(x)dx
dx dx dx
Z L0 Z L 0 0
du dv
EA · dx = p · v(x)dx + fL · v(L) + f0 · v(0)
0 dx dx 0

D’où le problème variationnelle :


Trouver u ∈ V vérifiant
(Pv )
a(u, v) = ℓ(v) ∀v ∈ V

où Z L
du dv
a(u, v) = EA · dx
0 dx dx
et Z L
ℓ(v) = p · v(x)dx + fL · v(L) + f0 · v(0)
0

Matrice élémentaire
On approche l’espace V par P1 = {ax + b, (a, b) ∈ R2 } et on considère un maillage de la barre [0, L] réduit
à un élément fini unique.
Les deux fonctions de base φi , i = 1, 2 définies par φi (aj ) = uij sont données pour tout x ∈ [0, L] par
x x
φ1 (x) = 1 − , φ2 (x) =
h h
Nous avons 2 2
h  Z h
−1
Z
dφ1 EA
a(φ1 , φ1 ) = EA dx = EA dx =
0 dx 0 L h
Z h 2 Z L  2
dφ2 1 EA
a(φ2 , φ2 ) = EA dx = EA dx =
0 dx 0 h h

1
 h    h   
−1
Z Z
dφ2 dφ1 1 EA
a(φ1 , φ2 ) = a(φ2 , φ1 ) = EA dx = EA dx = −
0 dx dx 0 h h h
Nous avons aussi ℓ(φ1 ) = f1 et ℓ(φ2 ) = f2 . D’où le système élémentaire
    
EA 1 −1 u1 f1
=
h −1 1 u2 f2

où f1 et f2 sont les efforts appliqués à la barre aux deux extrémités.

Barre à section variable

Soit (u0 , u1 , u2 , u3 , u4 ) les déplacements aux 5 nœuds. La matrice élémentaire de l’élément ei , i = 0, 3 est
donnée par     
EAi 1 −1 ui fi
=
hi −1 1 ui+1 fi+1
Soit explicitement pour chacune des 4 barres :
 EA0
− EA
   
h0 h0
0
0 0 0 u0 f0
 − 0EA EA
  u1  
0 0 0  f1
0    
 h0 h0 

 0 0   u
0 0 0  2  =
 f2 

 0 0 0 0 0   u 3
  f3 
0 0 0 0 0 u 4 f4
    
0 0 0 0 0 u0 f0
EA1

 0 h1 − EA
h1
1
0 0 
 u1  
  f1 


 0 − EA
h1
1 EA1
h1 0 0 
 u2 =
  f2 

 0 0 0 0 0  u3   f3 
0 0 0 0 0 u4 f4
    
0 0 0 0 0 u0 f0
 0 0 0 0 0  u1   f1 
EA2
− EA
    

 0 0 h2 h2
2
0 
 u2 =
  f2 

 0 0 − EA
h2
2 EA2
h2 0  u3   f3 
0 0 0 0 0 u4 f4
    
0 0 0 0 0 u0 f0

 0 0 0 0 0 
 u1  
  f1 

 0 0 0 0 0  u2 = f2 
EA3
− EA
    
 0 0 0 h3 h3
3  u3   f3 
0 0 0 − EA
h3
3 EA3
h3
u4 f4
La matrice d’assemblage s’obtient en sommant les 4 matrices élémentaires :
 EA
− EA

0 0
0 0 0
  
h0 h0 u0 = 0 f0
 − EA0 EA0 + EA1 − EA 1
0 0 u1 f1

h0 h0 h1 h1
    
EA1 EA1 EA2 EA2
  
0 − h1 + − 0 u2 = f2
  
 h1 h2 h2    
EA2 EA2 EA3 EA3   u3 f3
− h2 −

 0 0 h2 + h3 h3
   
0 0 0 EA3
− h3 EA3 u4 f4
h3

On u0 = 0, La première ligne nous donne la réaction à l’origine f0 = − EA

h0 u1 . Le système n’a que 4 inconnues,
0

on obtient le système d’ordre 4 en supprimant la première ligne et la première colonne :

 EA EA1

h0 + h1
0
− EA 1
0 0
  
h1 u1 f1

 − EA
h1
1 EA1
h1 + h2
EA2
− EA
h2
2
0 
 u2  =  f2 
  
EA2 EA2 EA3 EA3  
0 − h2 h2 + h3 − h3  u3 f3 
  

EA3 EA3 u4 f4
0 0 − h3 h3

2
Soit numériquement (E = 2 × 1011 Pa) :
1939 1210
 
π −
   
π 0 0 u1 20

 9 9 
   
 1210 8045 625    
 − π π − π 0   u2   0 
9 18 2
106   = 103 
    
 
 625 2125 125  
    
 0 − π π − π   u3   20 

 2 6 3 






 125 125 
u 10
0 0 − π π 4
3 3
La solution en (m) est :  
1 1 8237 1 1147237 1 1441267 1
, , ,
1620 π 9801000 π 1225125000 π 1225125000 π
soit en (mm) :
[0.1964875840, 0.2675154098, 0.2980731589, 0.3744675314]
Le déplacement à l’extrémité est alors u4 = 0.37mm.
La réaction f0 est donc f0 = − EA
h0 u1 = 50kN.
0

Poutre en flexion
On cherche une solution u ∈ V = H 2 (Ω). Multiplions l’équation EIu′′′′ = p par une fonction test v ∈ V :

d4 u
·v =p·v
EIz
Z L dx4 Z L
4
d u
EIz 4 · v(x) dx = p · v(x)dx
dx
Z L 03 0 L Z L
d3 u

du dv
−EIz ·
3 dx
dx + EI z · v(x) = p · v(x)dx
0 dx dx3 0 0
Z L 2 L  L Z L
d u d2 v d2 u dv d3 u

EIz ·
2 dx2
dx − EI z · + EIz · v(x) = p · v(x)dx
0 dx dx2 dx 0 dx3 0 0
L Z L
d2 u d2 v d2 u dv d2 u dv d3 u d3 u
Z
EIz · dx − EIz 2 · + EIz 2 · + EIz 3 v(L) − EIz 3 v(0) = p · v(x)dx
0 dx2 dx2 dx dx x=L dx dx x=0 dx x=L dx x=0 0
Z L 2 L
d u d2 v
Z
dv dv
EIz 2
· 2 dx − ML (L) − M0 (0) − FL v(L) − F0 v(0) = p · v(x)dx
Z0 L dx dx
L
dx dx 0
d2 u d2 v
Z
dv dv
EIz 2
· 2 dx = p · v(x)dx + F0 v(0) + FL v(L) + M0 (0) + ML (L)
0 dx dx 0 dx dx

D’où le problème variationnelle :


Trouver u ∈ V vérifiant
(Pv )
a(u, v) = ℓ(v) ∀v ∈ V

où
L
d2 u d2 v
Z
a(u, v) = EIz · dx
0 dx2 dx2
et Z L
dv dv
ℓ(v) = p · v(x)dx + F0 v(0) + FL v(L) + M0 (0) + ML (L)
0 dx dx
En notant Φi (x) les fonctions de base associées aux valeurs nodales de la fonction vi et Ψi (x) les fonctions
dv
de base associées aux valeurs nodales de la dérivée ( dx )i , on écrit :
n n
X X dv
v h (x) = vi Φi (x) + ( )i Ψi (x)
i=1 i=1
dx

Sur un élément ek = [xk , xk+1 ], cette approximation s’écrit :

dv dv
v h (x) = vk Φk (x) + ( )k Ψk (x) + vk+1 Φk+1 (x) + ( )k+1 Ψk+1 (x)
dx dx

3
 Φi et Ψi vérifient les conditions :
Φi (xj ) = uij Φ′i (xj ) = 0

Ψi (xj ) = 0 Ψ′i (xj ) = uij
D’après le cours :

Φi (x) = [1 − 2(x − xi )λ′i (xi )] λ2i (x)


Ψi (x) = (x − xi )λ2i (x)
x x x
Sur l’élément de référence x0 = 0 et x1 = 1, λ0 = 1 − L et λ1 = L D’où en posant L =ξ :


 Φ0 (x) = (1 + 2ξ)(1 − ξ)2 = 1 − 3ξ 2 + 2ξ 3
Φ1 (x) = (3 − 2ξ)ξ 2 = 3ξ 2 − 2ξ 3


 Ψ0 (x) = Lξ(1 − ξ)2 = L(ξ − 2ξ 2 + ξ 3 )
Ψ1 (x) = L(ξ − 1)ξ 2 = L(−ξ 2 + ξ 3 )

 ′
Φ (x) = −6ξ + 6ξ 2
 0′


Φ1 (x) = 6ξ − 6ξ 2
Ψ′ (x) = L(1 − 4ξ + 3ξ 2 )
 0′


Ψ1 (x) = L(−2ξ + 3ξ 2 )
 ′′
Φ (x) = −6 + 12ξ
 0′′


Φ1 (x) = 6 − 12ξ
 Ψ′′0 (x) = L(−4 + 6ξ)
 ′′

Ψ1 (x) = L(−2 + 6ξ)
Calcul de la matrice de rigidité :

L 2 1 2 1
d2 Φ1 d2 Φ1
Z  Z  Z
EI EI 2 EI
a11 = a(Φ0 , Φ0 ) = EI dx = 3 dξ = 3 (−6 + 12ξ) dξ = 12
0 dx2 L 0 dξ 2 L 0 L3
L 2 1 2
d2 Ψ0 d2 Ψ0 EI 1 2
Z Z Z
EI 2 EI
a22 = a(Ψ0 , Ψ0 ) = EI 2
dx = 3 2
dξ = 3 L (−4 + 6ξ) dξ = 4
0 dx L 0 dξ L 0 L
Z L 2 1
d Φ0 d2 Ψ0
Z
EI EI
a12 = a(Φ0 , Ψ0 ) = EI 2 2
dx = 3 (−6 + 12ξ) L(−4 + 6ξ)dξ = 6 2
0 dx dx L 0 L
Les autres coefficients se calculent de la même façon, d’où la matrice élémentaire :
 
12 6L −12 6L
EI  2 2 
 6L 4L −6L 2L 
L3  −12 −6L 12 −6L 
6L 2L2 −6L 4L2
dv dv
1. Si la poutre est bi-encastrée alors v(0) = v(L) = 0 et dx (0) = dx (L) = 0 donc
Z L Z L
dv dv
ℓ(v) = p · v(x)dx + F0 v(0) + FL v(L) + M0 (0) + ML (L) =⇒ ℓ(v) = p · v(x)dx
0 dx dx 0

Le second membre élémentaire est alors

  RL   R1   pL

p(x)Φ0 (x)dx pL 0 1 − 3ξ 2 + 2ξ 3 dξ

ℓ(Φ0 ) R 0L 2
pL2
R1
 ℓ(Ψ0 )  
  R0L p(x)Ψ0 (x)dx pL2 0 (ξ − 2ξ 2 + ξ 3 )dξ
    
12
 ℓ(Φ1 )  =  = =
     
R1 2 pL
pL 0 3ξ − 2ξ 3 dξ

 0 p(x)Φ1 (x)dx     2 
ℓ(Ψ1 ) RL R1 −pL2
0
p(x)Ψ1 (x)dx pL2 0 (−ξ 2 + ξ 3 )dξ 12

2. Si la poutre est encastrée à l’origine alors v(0) = 0 et si elle est soumise uniquement à une force ponctuelle
P⃗ à l’autre extrémité alors p = 0, ML = 0 donc
ℓ(v) = P · v(L)
Le second membre élémentaire est alors
     
ℓ(Φ0 ) P Φ0 (L) 0
 ℓ(Ψ0 )   P Ψ0 (L)   0 
 ℓ(Φ1 )  = 
   = 
P Φ1 (L)   P 
ℓ(Ψ1 ) P Ψ1 (L) 0

4
Poutre en flexion
matrice élémentaire
 12I 6I 12I 6I 
−
 L3 L2 L3 L2 
 
 
 6I 4I 6I 2I 

2
− 

 L L L2 L 
E
 

 
 12I 6I 12I 6I 
 − − − 2 
L3 L2 L3 L 
 

 
6I 2I 6I 4I
 
−
L2 L L2 L
Assemblage
12I0 6I0 12I0 6I0
 
−

 L30 L20 L30 L20 

 
6I0 4I0 6I0 2I0
 
−
    
F0 
 L20 L0 L20 L0

 u0 = 0
    
    
 M0  

 θ = 0
 0 
  
 12I0 6I0 12I0 12I1 6I0 6I1 12I1 6I1
    
 − 3 − + 3 − + 2 −
   
L20 L30 L20 L31 L21




  L0 L1 L1  

 0  
   u1
 
  
   6I0 2I0 6I0 6I1 4I0 4I1 6I1 2I1  
− + 2 + −
    
 0  L20 L0 L20 L1 L0 L1 L21 L1   θ1
  
   
    
 =E


 
    
  
 0   12I1 6I1 12I1 12I2 6I1 6I2 12I2 6I2   u2
 
   − − + 3 − + 2 − 
L31 L21 3
L1 L2 L12 L2 L32 L22
    
    
 0  
   θ2
 
6I1 2I1 6I1 6I2 4I1 4I2 6I2 2I2
    
− −
  
 + 2 + 
L21 2 L22
  





 L1 L1 L2 L1 L2 L2 
 

 P  
   u3
 
    
    
0
 12I2 6I2 12I2 6I2 
θ3
 − − − 

 L32 L22 L32 L22 

 
 
 6I2 2I2 6I2 4I2 
−
L22 L2 L22 L2

12I0 12I1 6I0 6I1 12I1 6I1

 
+ 3 − + 2 −
 L30 L1 L20 L1 L31 L21 
 
    
0 6I0 6I1 4I0 4I1 6I1 2I1 u1
 
 − 2 + 2 + − 2
 




  L 0 L1 L0 L1 L1 L1  


 0 




 θ1 

    
    
   12I1 6I1 12I1 12I2 6I1 6I2 12I2 6I2  
− 3 − + 3 − + 2 −
    
0 L1 L21 L31 L2 L21 L2 L32 L22 u2 
   
    
 =E


 
  
 0   6I1 2I1 6I1 6I2 4I1 4I2 6I2 2I2 
θ2 
   − + 2 + −  
L21 L1 L12 L2 L1 L2 L22 L2
    
    
    
    
P u3 
   
 
12I2 6I2 12I2 6I2
    
− − −
    
0

 L32 L22 L32 L22 
 θ3
 
 
 6I2 2I2 6I2 4I2 
−
L22 L2 L22 L2

5
Application numérique

103515625π
− 4921875π − 15625000π
    
0 3 2 3 468750π 0 0 u1

 0  
  − 4921875π
2 446875π −468750π 28125π 0 0  θ1 
− 15625000π 78125000π
− 781250π − 31250000π 625000π
 

 0  
= 3 −468750π 9 3 9 3

 u2 


 0  
  468750π 28125π − 781250π
3
218750π
3 − 625000π
3
25000π
3

 θ2 

 −1000   0 0 − 31250000π
9 − 625000π
3
31250000π
9 − 625000π
3
 u3 
625000π 25000π
0 0 0 3 3 − 625000π
3
50000π
3
θ3

Solution :
104 64 3184 892 1442 1567
 
u=− 234375π 15625π 1171875π 46875π 234375π 46875π

soit
 
u≃− 0.0001412447068 0.001303797293 0.0008648522046 0.006057224925 0.001958412184 0.01064088728

1442 1442
u3 = − ≃ −0.001958412184 m ≃ −2 mm , θ3 = − ≃ −0.01064088728 rd ≃ −0.6◦
234375π 234375π