Treillis plans à nœuds articulés 3

Introduction

Un treillis est un ensemble de poutres droites (éléments) reliées entre elles par des rotules (nœuds).
Les liaisons extérieures sont des rotules et des appuis simples. Les charges sont des forces portées par
les rotules, des gradients thermiques et des déplacements d’appui. La force intérieure dans une
section droite se réduit à l’effort normal.

Le treillis est plan si :

– Le plan {O; x, y} est un plan de symétrie pour toutes les sections droites.
– Les forces appliquées sont situées dans le plan {O; x, y}.

On suppose que les déplacements sont petits.

1 Matrices élémentaires
Soit (i → j) un élément de treillis plan de section droite constante (figure 1).

Figure 1 – Élément i → j

L est la longueur de l’élément et A l’aire de sa section droite.

(xi , yi ) et (xj , yj ) sont les coordonnées des nœuds de l’élément.

Le vecteur unitaire ~n porté par l’axe de la poutre est défini par :

½ ¾ ½ ¾ ½ ¾
nx 1 xj − xi cos θ
= = , L2 = (xj − xi )2 + (yj − yi )2 (1.1)
ny L yj − yi sin θ

où θ est l’angle que fait ~n avec l’axe x.

E et α sont respectivement le module de Young et le coefficient de dilatation du matériau.

L’élément est soumis à un effort normal N (positif : traction, négatif : compression) et à une variation
de température ∆T constante.
4 Méthode des éléments finis

(ui , vi) et (uj , vj ) sont les déplacements nodaux (figure 2).

Figure 2 – Déplacements élémentaires

Les efforts aux extrémités de l’élément sont (figure 3) :

− N ~n en i , N ~n en j (1.2)

Figure 3 – Efforts élémentaires

Différentions la relation :
L2 = (xj − xi )2 + (yj − yi )2 (1.3)
Il vient :
2 L dL = 2 (xj − xi ) (dxj − dxi ) + 2 (yj − yi ) (dyj − dyi ) (1.4)
d’où l’expression de l’allongement unitaire suivant ~n :
µ ¶
dL 1 (xj − xi ) (yj − yi )
εn = = (dxj − dxi ) + (dyj − dyi ) (1.5)
L L L L
soit :
1
εn = ( nx (uj − ui ) + ny (vj − vi ) ) (1.6)
L
Cet allongement unitaire est dû à l’effort normal (loi de Hooke) et à la variation de température :
N
εn = + α ∆T (1.7)
EA
L’effort normal s’écrit en fonction des déplacements nodaux :

N = EA (εn − α ∆T )
EA (1.8)
= ( nx (uj − ui ) + ny (vj − vi ) ) − EA α ∆T
L
soit :  
 ui 
 
EA £ ¤  vi 
N= −nx −ny nx ny − EA α ∆T (1.9)
L 
uj 
 
vj
Treillis plans à nœuds articulés 5

On en déduit :
{fnod } = [ k ] {u} − {fth } (1.10a)
avec :    

 −nx 
 
 ui 
    
−ny vi
{fnod } = N , {u} = (1.10b)

 n   u 
 x    j
 
ny vj
   2 

 −nx 
 nx nx ny −n2x −nx ny
EA −ny  £ ¤ EA 
 n2y −nx ny −n2y  
[k] = −nx −ny nx ny =  2  (1.10c)
L 
 nx 
 L n x nx n y
  2
ny sym. ny
 

 −nx 

 
−ny
{fth } = EA α ∆T (1.10d)

 n 
 x  
ny
{fnod } est le vecteur force nodal (N).

{u} est le vecteur déplacement élémentaire (m).

[ k ] est la matrice de rigidité (N/m).

{fth } est le vecteur force équivalent au gradient thermique (N).

Remarque 1 : la matrice de rigidité peut se mettre sous la forme :

· ¸ · ¸
[k̄] −[k̄] EA n2x nx ny
[k] = avec [k̄] = (1.11)
−[k̄] [k̄] L nx ny n2y

Remarque 2 : la matrice de rigidité est égale à :

 
nx 0 · ¸· ¸
ny 0  EA 1 −1 nx ny 0 0
[k] =   (1.12)
0 nx  L −1 1 0 0 nx ny
0 ny
ou :
    
nx −ny 0 0 1 0 −1 0 nx ny 0 0
 ny nx 0 0  EA 0 0 0 0 −ny nx 0 0
[k] = 
0
    (1.13)
0 nx −ny  L −1 0 1 0  0 0 nx ny 
0 0 ny nx 0 0 0 0 0 0 −ny nx

Remarque 3 : l’énergie de déformation est égale à (à un coefficient près indépendant des déplacements
et de leurs dérivées) :

1
Edef = EA ε2n L − EA εn α ∆T L
2 (1.14)
1
= {u}T [ k ] {u} − {u}T {fth }
2

Le travail des forces extérieures se réduit au travail des forces nodales :

Wext = {u}T {fnod } (1.15)

6 Méthode des éléments finis

L’énergie potentielle est égale à :

Epot = Edef − Wext (1.16)

La matrice de rigidité est la matrice hessienne (ou matrice de Hess) de l’énergie de déformation par
rapport aux déplacements nodaux (programme tre mat) :

∂ 2 Edef
kij = (= kji ) (1.17)
∂ui ∂uj

Le vecteur des efforts aux nœuds est le gradient de l’énergie de déformation par rapport aux dépla-
cements nodaux :
∂Edef
fnod,i = (1.18)
∂ui

2 Exemple 1 : treillis soumis à une force nodale

2.1 Énoncé
Le treillis plan à noeuds articulés représenté sur la figure 4 est com-
posé de trois poutres de même nature et de même section droite.

Soient E le module de Young du matériau et A l’aire des sections

droites.

Le noeud 1 est articulé et le nœud 3 repose sur un appui simple dont

la normale est horizontale.

Le noeud 2 porte une charge de composantes (0, P ).

Application numérique : on donne :

A = 100 mm2 , L = 0.2 m , E = 200000 MPa , P = −10000 N Figure 4 – Exemple 1

2.2 Partition des degrés de liberté

Effectuons une partition des degrés de liberté en déplacements connus et inconnus ([1], [12]) :
 

 u2 
    
 v2 
u2  u1  ½ ¾ 
 


{UL } v3
{UL } = v2 , {US } = v1 d’où {U } = =
    {US } u1 
 
v3 u3 
 


 
v 1 
 
u3

On en déduit la localisation des degrés de liberté dans les matrices globales :

 

u1 → 0 

 v → 0 


 1 

 
u2 → 1
{DDL} =

 v2 → 2 

 

u3 → 0 

 
v3 → 3
Treillis plans à nœuds articulés 7

2.3 Étude élémentaire

– coordonnées nodales :
nœud x y
1 0 L
2 L 0
3 0 −L
– élément 1 → 2 :
√ 1 1
caractéristiques : L 2 , A , E , nx = √ , ny = − √
2 2
   

 u1 → 0 1 −1 −1 1
  EA 
v1 → 0  −1 1 1 −1
{ddl1−2 } = , [k1−2 ] = √ 

 u → 1
 2 2L −1 1 1 −1
 2 
v2 → 2 1 −1 −1 1
– élément 3 → 1 :
caractéristiques : 2 L , A , E , nx = 0 , ny = 1
   

 u3 → 0  0 0 0 0
  EA 
v3 → 3 0 1 0 −1

{ddl3−1 } = , [k3−1 ] =

 u → 0  2 L 0 0 0 0
 1 
v1 → 0 0 −1 0 1
– élément 3 → 2 :
√ 1 1
caractéristiques : L 2 , A , E , nx = √ , ny = √
2 2
   

 u3 → 0  1 1 −1 −1
  
v3 → 3 EA  1 1 −1 −1
{ddl3−2 } = , [k3−2 ] = √ 

 u → 1
 2 2L −1 −1 1 1
 2 
v2 → 2 −1 −1 1 1

2.4 Assemblage et calcul des déplacements inconnus

Les déplacements inconnus sont les solutions de l’équation [KLL ]{UL } = {Fnod,L } :
    
2 0 −1 u2   0 
EA 
√ 0 2 −1√  v2 = P
2 2 L −1 −1 1 + 2  v   0 
3

d’où (programme exemple 1) :

PL PL √
u2 = = −0.050 mm , v2 = (1 + 2 2) = −0.191 mm
2 EA 2 EA
PL
v3 = = −0.100 mm
EA

2.5 Efforts normaux dans les éléments

Ils sont calculés à l’aide de la formule (1.8) :
µ ¶ √
EA 1 1 P 2
N1−2 = √ √ u2 − √ v2 = − = 7071 N
L 2 2 2 2
EA P
N3−1 = ( −v3 ) = − = 5000 N
2L 2
µ ¶ √
EA 1 1 P 2
N3−2 = √ √ u2 + √ ( v2 − v3 ) = = −7071 N
L 2 2 2 2
8 Méthode des éléments finis

2.6 Actions de liaison

Les actions de liaisons sont calculées à partir des efforts normaux (équation 1.2) :

– nœud 1 :
F~1 = −N1−2 ~n1−2 + N3−1 ~n3−1 d’où

1 P 1
F1x = − √ N1−2 = = −5000 N , F1y = √ N1−2 + N3−1 = −P = 10000 N
2 2 2
– nœud 3 :
F~3 = −N3−1 ~n3−1 − N3−2 ~n3−2 d’où

1 P
F3x = − √ N3−2 = − = 5000 N
2 2
Remarque : l’équilibre de la structure est vérifié :

F1x + F2x + F3x = 0 , F1y + F2y + F3y = 0 , −2 LF1x − LF2x + LF2y = 0

3 Exemple 2 : treillis soumis à une force nodale

3.1 Énoncé
Le treillis plan représenté sur la figure 5 est composé de trois
poutres de même section.

Soient E le module de Young du matériau et A l’aire des sec-

tions droites.

Le nœud 1 est articulé et le nœud 2 repose sur un appui simple

dont la normale est horizontale.

Le nœud 3 porte une charge d’intensité (P, 3 P, 0).

Figure 5 – Exemple 2
Application numérique : on donne :

A = 100 cm2 , E = 200000 MPa , L = 0.7 m , P = −120 kN

3.2 Partition des degrés de liberté

Effectuons une partition des degrés de liberté en déplacements connus et inconnus ([1], [12]) :
 

 v2 
    u 
 
¾  3
 v2  u1  ½   
{UL } v3
{UL } = u3 , {US } = v1 d’où {U } = =
    {US } u1 
 
v3 u2 
 

 v1 

  
u2
Treillis plans à nœuds articulés 9

On en déduit la localisation des degrés de liberté dans les matrices globales :

 

 u1 → 0

 v → 0 


 1 

 
u2 → 0
{DDL} =

 v2 → 1 


 

 u3 → 2

 
v3 → 3

3.3 Étude élémentaire

– coordonnées nodales :
nœud x y
1 0 0
2 0 L
3 L 0
– élément 1 → 2 :
caractéristiques : L , A , E , nx = 0 , ny = 1
   

 u1 → 0  0 0 0 0
  EA  
v1 → 0 0 1 0 −1
{ddl1−2 } = , [k1−2 ] =

 u → 0 

L 0 0 0 0
 2 
v2 → 1 0 −1 0 1
– élément 1 → 3 :
caractéristiques : L , A , E , nx = 1 , ny = 0
   

 u 1 → 0
 1 0 −1 0
  
v1 → 0 EA  0 0 0 0

{ddl1−3 } = , [k1−3 ] =

 u → 2
 L −1 0 1 0
 3 
v3 → 3 0 0 0 0
– élément 2 → 3 :
√ 1 −1
caractéristiques : L 2 , A , E , nx = √ , ny = √
2 2
   

 u2 → 0  1 −1 −1 1
  EA 
v2 → 1 −1 1 1 −1

{ddl2−3 } = , [k2−3 ] = √ 

 u → 2  2 2 L −1 1 1 −1
 3 
v3 → 3 1 −1 −1 1

3.4 Assemblage et calcul des déplacements inconnus

Les déplacements inconnus sont les solutions de l’équation [KLL ]{UL } = {Fnod,L } :
 √    
1+2 2 1 √ −1  v2   0 
EA 
√ 1 1 + 2 2 −1 u3 = P
2 2L    
−1 −1 1 v3 3P

d’où (programme exemple 2) :

3PL
v2 =
= −0.126 mm
EA
√
4PL (7 + 6 2) P L
u3 = = −0.168 mm , v3 = = −0.650 mm
EA EA
10 Méthode des éléments finis

3.5 Efforts normaux dans les éléments

Ils sont calculés à l’aide de la formule (1.8) :
EA
N1−2 = v2 = 3 P = −360 kN
L
EA
N1−3 = u3 = 4 P = −480 kN
L
µ ¶
EA 1 1 √
N2−3 = √ √ u3 − √ ( v3 − v2 ) = −3 2 P = 509 kN
L 2 2 2

3.6 Actions de liaison

Les actions de liaisons sont calculées à partir des efforts normaux (équation 1.2) :

– nœud 1 :
F~1 = −N1−2 ~n1−2 − N1−3 ~n1−3 d’où
F1x = −N1−3 = −4 P = 480 kN , F1y = −N1−2 = −3 P = 360 kN
– nœud 2 :
F~2 = N1−2 ~n1−2 − N2−3 ~n2−3 d’où
1
F2x = − √ N2−3 = 3 P = −360 kN
2
Remarque : l’équilibre de la structure est vérifié :

F1x + F2x + F3x = 0 , F1y + F2y + F3y = 0 , LF3y − LF2x = 0

4 Exemple 3 : treillis soumis à une variation de température

4.1 Énoncé
Le treillis plan à noeuds articulés représenté sur la figure 6 est composé de trois poutres de même
matériau et de même section droite.

Figure 6 – Exemple 3

Soient E et α respectivement le module de Young et le coefficient de dilatation du matériau.

Soit A l’aire des sections droites.

Treillis plans à nœuds articulés 11

Les nœuds 1, 2 et 4 sont liés à l’extérieur par une rotule.

Premier cas de charge : la structure est soumise à une variation de température ∆T .

Deuxième cas de charge : la poutre (2 − 3) est soumise à une variation de température ∆T .

Application numérique : on donne :

A = 100 mm2 , L = 0.1 m , E = 200000 MPa , α = 10−5 K−1 , ∆T = 100 K

4.2 Partition des degrés de liberté

Effectuons une partition des degrés de liberté en déplacements connus et inconnus ([1], [12]) :
 
  
 u3 
u1  
 v 


  
 3 


 v1 
 
 u 


   1 
½ ¾    ½ ¾   
u3 u2 {UL } v1
{UL } = , {US } = d’où {U } = =
v3  v2 
  {US } u2 
 

 
 
 


 
u4  
 
v2 
  
 
v4 
 u4 
  
v4

On en déduit la localisation des degrés de liberté dans les matrices globales :

 

 u1 → 0
v → 0
 


 1 


 u → 0

 2
 

 
v2 → 0
{DDL} =

 u3 → 1 

 

 v3 → 2 


 

 u4 → 0 

 
v4 → 0

4.3 Étude élémentaire

– coordonnées nodales :
nœud x y
1 0 0
2 L 0
3 L L
4 2L L
– élément 1 → 3 :
√ 1 1
caractéristiques : L 2 , A , E , α , nx = √ , ny = √
2 2
     

 u1 → 0 1 1 −1 −1 
 −1
  EA    
v1 → 0 1 1 −1 −1  {fth,1−3 } = √1 EA α ∆T −1
{ddl1−3 } = [k1−3 ] = √  

 u → 1 
 2 2L −1 −1 1 1 2 1
 
 3   
v3 → 2 −1 −1 1 1 1
– élément 2 → 3 :
caractéristiques : L , A , E , α , nx = 0 , ny = 1
12 Méthode des éléments finis
     

u 2 → 0 
 0 0 0 0 
 0
  EA    
v2 → 0 0 1 0 −1
 −1
{ddl2−3 } = [k2−3 ] = {fth,2−3 } = EA α ∆T

u → 1  L 0 0 0 0 
 0
 3    
v3 → 2 0 −1 0 1 1
– élément 3 → 4 :
caractéristiques : L , A , E , α , nx = 1 , ny = 0
     

u3 → 1  1 0 −1 0 
−1
     
v3 → 2 EA  0 0 0 0  {fth,3−4 } = EA α ∆T 0
{ddl3−4 } = [k3−4 ] =

u → 0 
 L −1 0 1 0 
 1
 4    
v4 → 0 0 0 0 0 0

4.4 Premier cas de charge

4.4.1 Assemblage et calcul des déplacements inconnus
Les déplacements inconnus sont les solutions de l’équation [KLL ]{UL } = {Fth,L } :
· √ ¸½ ¾ ½√ ¾
EA 1 + 2 2 1√ u3 1 2 − 2
√ = EA α ∆T √
2 2L 1 1 + 2 2 v3 2 2+2
d’où (programme exemple 3) :
√ √
u3 = ( 2 − 2) L α ∆T = −0.0586 mm , v3 = 2 L α ∆T = 0.1414 mm

4.4.2 Efforts normaux dans les éléments

Ils sont calculés à l’aide de la formule (1.8) :
µ ¶
EA 1 1 √
N1−3 = √ √ u3 + √ v3 − EA α ∆T = ( 2 − 2) EA α ∆T = −11716 N
L 2 2 2
EA √
N2−3 = v3 − EA α ∆T = ( 2 − 1) EA α ∆T = 8284 N
L
EA √
N3−4 = (−u3 ) − EA α ∆T = (1 − 2) EA α ∆T = −8284 N
L

4.4.3 Actions de liaison

Les actions de liaisons sont calculées à partir des efforts normaux (équation 1.2) :

– nœud 1 :
F~1 = −N1−3 ~n1−3d’où :
1 √
F1x = − √ N1−3 = ( 2 − 1) EA α ∆T = 8284 N , F1y = F1x = 8284 N
2
– nœud 2 :
F~2 = −N2−3 ~n2−3 d’où :
√
F2x = 0 , F2y = −N2−3 = (1 − 2) EA α ∆T = −8284 N
– nœud 4 :
F~4 = N3−4 ~n3−4 d’où :
√
F4x = N3−4 = (1 − 2) EA α ∆T = −8284 N , F4y = 0
Remarque : l’équilibre de la structure est vérifié :
F1x + F2x + F4x = 0 , F1y + F2y + F4y = 0 , L F2y − L F4x = 0
Treillis plans à nœuds articulés 13

4.5 Deuxième cas de charge

4.5.1 Assemblage et calcul des déplacements inconnus
Les déplacements inconnus sont les solutions de l’équation [KLL ]{UL } = {FL } :
· √ ¸½ ¾ ½ ¾
EA 1 + 2 2 1√ u3 0
√ = EA α ∆T
2 2L 1 1 + 2 2 v3 1

d’où (programme exemple 3) :

√ √
1− 2 3− 2
u3 = L α ∆T = −0.0207 mm , v3 = L α ∆T = 0.0793 mm
2 2

4.5.2 Efforts normaux dans les éléments

Ils sont calculés à l’aide de la formule (1.8) :
µ ¶√
EA 1 1 2− 2
N1−3 = √ √ u3 + √ v3
= EA α ∆T = 5858 N
L 2 2 2 2
√
EA 1− 2
N2−3 = v3 − EA α ∆T = EA α ∆T = −4142 N
L 2
√
EA 2−1
N3−4 = (−u3 ) = EA α ∆T = 4142 N
L 2

4.5.3 Actions de liaison

Les actions de liaisons sont calculées à partir des efforts normaux (équation 1.2) :

– nœud 1 :
F~1 = −N1−3 ~n1−3 d’où :
√
1 1− 2
F1x = − √ N1−3 = EA α ∆T = −4142 N , F1y = F1x = −4142 N
2 2
– nœud 2 :
F~2 = −N2−3 ~n2−3 d’où :
√
2−1
F2x = 0 , F2y = −N2−3 = EA α ∆T = 4142 N
2
– nœud 4 :
F~4 = N3−4 ~n3−4 d’où :
√
2−1
F4x = N3−4 = EA α ∆T = 4142 N , F4y = 0
2
Remarque : l’équilibre de la structure est vérifié :

F1x + F2x + F4x = 0 , F1y + F2y + F4y = 0 , L F2y − L F4x = 0