DIMENSIONNEMENT DES STRUCTURES

APPLICATION A LA MODELISATION DES STRUCTURES PAR ELEMENTS FINIS

ex-ef-souplesse-poutre-encastree.doc/version 01/11/2010/JG

DETERMINATION DE LA MATRICE DE SOUPLESSE D’UNE STRUCTURE

1 OBJECTIFS
Calcul de la matrice de souplesse à partir d’un exemple
Soit une structure poutre encastrée à une extrémité et sollicitée par deux forces F1 et F2 ayant pour caractéristique
géométrique un moment quadratique autour de l'axe z Iz, noté ici I et pour caractéristique de matériau un module de Young E.
y

G x

x 1 2
F1 F2
x1
x2
2 TRAVAIL DEMANDE
On souhaite déterminer les déplacements verticaux (flèches) au droit des points de chargements 1 et 2 et les mettre
sous la forme de la relation de comportement déplacements-forces {d} = [α ]• {F}
Méthodologie
Pour chaque chargement, calculer les déplacements aux points 1 et 2 et superposer les effets de ces deux chargements cumulés
afin de calculer les flèches totales, mettre ensuite ces résultats sous le forme matricielle {d} = [α ]• {F}
y

G x

x 1 2
F1
x1
x2

G x

x 1 2
F2
x1
x2

1
ELEMENTS DE CORRECTION
Premier cas de chargement : F1 = F1 y à l'abscisse x1
y
v21
G θz1
v11 x

x 1 2
F1
x1
x2
Eléments de réduction du torseur de cohésion au centre géométrique G d'une section courante d'abscisse 0≤x≤x1:
 
{Coh }G = T y = F1 y 
Mf z = F1 ( x1 − x ) z G

d2v Mf z
La relation de comportement au moment de flexion s'écrit : =
2 EI
dx
où v(x) est le déplacement suivant y de tout point de la ligne moyenne de la poutre.
d2v
Dans la zone 0 ≤ x ≤ x1 la déformée v(x) vérifie donc l'équation différentielle : EI = Mf z = F1 ( x 1 − x )
dx 2
donc par intégrations successives (les constantes d'intégration disparaissant du fait de l'encastrement en x=0) :
dv x2 x3 x2 x2  x
EI = − F1 + F1 x 1 x et ensuite EIv( x ) = −F1 + F1 x 1 = F1  x1 − 
dx 2 6 2 2  3
- déplacements à l'abscisse x1 :
x 13 dv( x 1 ) x 2
une flèche : v( x 1 ) = F1 et une rotation de section : θ z ( x 1 ) ≅ = F1 1
3EI dx 2EI
- déplacement à l'abscisse x2 > x1 s'écrit alors, en remarquant que la portion de poutre x1 ≤ x ≤ x2 reste rectiligne :
dv( x 1 ) x 3 x 2
v( x 2 ) = v( x 1 ) + ( x 2 − x 1 ) ⇒ v( x 2 ) = F1 1 + F1 1 ( x 2 − x 1 )
dx 3EI 2EI
x 2 x 
soit : v( x 2 ) = F1 1  x 2 − 1 
2EI  3 
Avec la notation indicielle les déplacements obtenus pour l’application de F1 s’écrivent :
x 3
au point (1) : v11 = F1 1
3EI
x 2 x 
au point (2) : v 21 = F1 1  x 2 − 1 
2EI  3 
Second cas de chargement : F2 = F2 y à l'abscisse x2 :
θ z2

y
v22
G
v12 x

x 1 2
F2
x1
x2

Considérons maintenant l'effort F2 à l'abscisse x2. En tenant compte que x 2 > x 1 , on peut utiliser l'expression
3 2 2
x x x  x
EIv( x ) = −F1 + F1 x 1 = F1  x 1 −  dans laquelle on remplace F1 par F2 et x1 par x2.
6 2 2  3
Avec la notation indicielle les déplacements obtenus pour l’application de F2 s’écrivent :
x 2 x 
au point 1 : v12 = F2 1  x 2 − 1 
2EI  3 

2
 x 23
au point 2 : v 22 = F2
3EI
Superposition des effets des deux chargements : F1 = F1 y à l'abscisse x1 et F2 = F2 y à l'abscisse x2 :
- déplacements totaux :
y

v22
v2
v11 v12

v21
v1
x

1 2
x1 F1 F2
x2

x 13 x 2 x 
v1 = v11 + v12 = F1 + F2 1  x 2 − 1 
3EI 2EI  3 
x 12  x  x 3
v 2 = v 21 + v 22 = F1  x 2 − 1  + F2 2
2EI  3  3EI
 x 13 x12  x 
v1    x 2 − 1  F1 
  3   
soit sous forme matricielle :   =  2
3EI 2EI 
⇔ {d} = [α ]• {F}
  
•
x  3
v   1  x − x 1  x2 
 F2 
 2  2 3 
 2EI   3EI 
on peut remarquer la symétrie des coefficients α 12 = α 21 .