Exo8 Corige

Les treillis par la MEF
Exercice 8 : Modélisation EF d'une colonne sous son poids propre

Objectifs : Notion d'erreur de discrétisation, et analyse des résultats EF.
Étude de convergence en affinant le maillage.
x 1- Établir l'équation matricielle d'un modèle à un élément fini
Analyser les résultats aux nœuds (déplacement & efforts)
Tracer les résultats sur l'élément (déplacement & efforts)
g ℓ = 6h 2- Construire le modèle utilisant deux éléments finis identiques
Analyser les résultats aux nœuds (déplacement & efforts)
Tracer les résultats sur les éléments (déplacement & efforts)
3- Modèle à 3 éléments, pour affiner le maillage dans la zone la plus contrainte
nous utilisons 3 éléments de longueur h, 2h, et 3h.
Déduire des calculs précédents l'équation matricielle du modèle.
Tracer les résultats sur les éléments (déplacement & efforts)
Pour améliorer la solution éléments finis nous avons augmenté le nombre d'éléments et densifié le
maillage dans la zone la plus chargée. Cette méthode dite « h convergence » demande en général un
nombre élevé d'éléments finis.
La figure suivante présente les résultats d’un modèle éléments finis en contraintes planes. Pour quantifier l’erreur relative à
cette discrétisation, la démarche est identique à celle de cet exercice, elle est basée sur l’analyse de la discontinuité du
champ des contraintes entre deux éléments adjacents.
en MPa
Dans cette section le diagramme des contraintes est le suivant

σ VM solution cherchée
Discontinuité
145
solution éléments finis
constante par morceau L’erreur est beaucoup trop importante.
83
Ce modèle n’est pas satisfaisant, il faut
62
affiner le maillage
Rappel : la solution analytique est donnée par

 x  mg  x 
L’effort normal : N = − mg 1 −  Le déplacement : u = − x 1 − 
 6h  ES  12h 
21
1- Modèle à 1 élément fini
Modèle à 2 degrés de liberté {U } = {u1 u2 }

T
xo
2
1 − 1
[K ] = [K1 ] = ES 
g 6h − 1 1 
(1) ρ gS
 6h 
  mg 1
1
{Fd } = − ρ gS  2  = −  
 6h  2 1
X1  2 
 X1 
La condition de déplacement imposé u1 = 0 ==> {Fi } =   effort de liaison inconnu
0
ES  1 −1  0  mg 1  X1 
D’où l'équation matricielle à résoudre :  −1 1  u  = − 2 1 +  0 
6h   2   
La solution est :
3mgh mg
Champ de déplacement : u2 = − soit une approximation u ( x ) = − x
ES 2 ES
Effort à l’encastrement : X1 = mg
ES mg
L’état de contrainte sur l'élément est une constante : N = u2 = −
6h 2
La solution éléments finis est exacte aux nœuds
Elle donne une approximation sur l'élément.
Courbes représentant le champ de déplacement et le champ de contrainte.
Champ de déplacement Champ de contrainte

Le saut de contrainte entre l'information nodale et l'information élémentaire
permet de dire que cette modélisation n’est pas satisfaisante
22
2- Modèle à 2 éléments finis : deux éléments de longueur identique
Ce modèle à 3 degrés de liberté {U } = {u1 u2 u3 }

T
xo
3
ES  1 −1 mg 1
pour chaque élément : [ K e ] =  −1 1  et {Fd } = −   (2)
3h   4 1 g 6h 2 ρ gS
L’équation matricielle obtenue après assemblage est : (1)
 1 −1 0   0  1   X1  1
ES     mg    
−1 2 −1 u2  = − 2  +  0 
3h    4     X
1
 0 −1 1  u3  1   0 
D’où la solution
9mgh 3mgh
Déplacements nodaux : u2 = − et u3 = −
4 ES ES
Effort à l’encastrement : X1 = mg
Post-traitement
 ES 3mg
 N1 = 3h u2 = − 4
État de contrainte sur les éléments : 
 N = ES ( u − u ) = − mg
 2 3h 3 2 4
Isolons l'élément N°1 pour calculer l'effort au nœud interne N°2
ES  1 −1  0  mg 1  X1 
 −1 1  u  = − 4 1 +  X  X 21 = − mg / 2
3h   2    21 
Reportons ces résultats sur les courbes de la solution analytique
Nous retrouvons les informations exactes (déplacement et efforts) aux nœuds. Une meilleure
approximation des champs de déplacements et de contrainte sur les éléments, l’erreur sur la
contrainte maximale est maintenant de 25%. Cette modélisation n’est toujours pas satisfaisante.
23
3- Modèle à 3 éléments xo
L’équation matricielle obtenue après assemblage de ce modèle à 4 degré de 4
liberté est : (3)
 6 −6 0 0   0  1   X1  g 6h 3 ρ gS
 −6 9 −3 0  u      (2)
  2  = − mg 3 +  0 
ES  2
(1)
6h  0 −3 5 −2  u3  12 5  0  1
   3  0 
 0 0 −2 2  u4  X
1
Sur {U } = {u1 u2 u4 }
T
u3
D’où la solution
11mgh 3mgh 3mgh
• Déplacements nodaux : u2 = − , u3 = − et u4 = −
12 ES 4 ES ES
• Effort à l’encastrement : X1 = mg
 ES 11mg
 N1 = h u2 = − 12

 ES 2mg
• État de contrainte sur les éléments :  N 2 = ( u3 − u2 ) = −
 2h 3
 ES mg
 N3 = 3h ( u4 − u3 ) = − 4

Reportons ces résultats sur les courbes de la solution analytique
L’approximation des champs de déplacements et contraintes sur les éléments est meilleur, erreur de
8% sur la contrainte maximale, on note cependant que la convergence est lente.
Pour améliorer la solution éléments finis nous avons augmenté le nombre d'éléments et densifié le
maillage dans la zone la plus chargée. Cette méthode dite « h convergence » demande en général
un nombre élevé d'éléments finis. Il est possible d'utiliser des éléments de degré 2 qui autorisent
une variation linéaire de la contrainte sur l'élément.
24

Exo8 Corige

Transféré par

Droits d'auteur :

Formats disponibles

Exo8 Corige

Transféré par

Informations du document

Titre original

Copyright

Formats disponibles

Partager ce document

Partager ou intégrer le document

Options de partage

Avez-vous trouvé ce document utile ?

Ce contenu est-il inapproprié ?

Droits d'auteur :

Formats disponibles

Exo8 Corige

Transféré par

Droits d'auteur :

Formats disponibles

Les treillis par la MEF

Exercice 8 : Modélisation EF d'une colonne sous son poids propre

Dans cette section le diagramme des contraintes est le suivant

Rappel : la solution analytique est donnée par

1- Modèle à 1 élément fini

Modèle à 2 degrés de liberté {U } = {u1 u2 }

Champ de déplacement Champ de contrainte

2- Modèle à 2 éléments finis : deux éléments de longueur identique

Ce modèle à 3 degrés de liberté {U } = {u1 u2 u3 }

Champ de déplacement Champ de contrainte

Champ de déplacement Champ de contrainte

Vous aimerez peut-être aussi