IMC150 – RÉSISTANCE DES STRUCTURES II
COURS 8 – MÉTHODE DES ÉLÉMENTS
FINIS – 2E PARTIE
Walid Belgacem
15 novembre 2011
Mise en contexte
Treillis (1D et 2D)
Traction-compression seulement
Efforts
Étude sur la fibre neutre
Propriétés de section (aire, densité, …)
3 degrés de liberté par nœuds
3 translation
COMSOL : Treillis dans le plan
Cours 8 – Méthode des éléments finis II
Mise en contexte
Poutres (1D et 2D)
Traction, compression et flexion
Efforts et moments
6 degrés de liberté par nœuds
3 translation + 3 rotation
Étude sur la fibre neutre
Définis la section (géométrie et inertie)
COMSOL : Poutre d’Euler dans le plan
Cours 8 – Méthode des éléments finis II
MEF pour les treillis 2D
Au cours 7 : Treillis 1 D
On considère le treillis 2D dans un système d’axes
global (x , y).
On cherche le déplacement à chacun des nœuds
dans le repère global est noté (u ,v).
y 2
P
c b
1 a 3
x
Cours 8 – Méthode des éléments finis II
MEF pour les treillis 2D
On veut solutionner [K]{u}={f}, dans le système
global :
On débute en cherchant [K] et {f} dans tous les
systèmes locaux, pour chaque barre.
On suppose une barre quelconque : v2 x
v2 u2
Repère local (x,y)
u2
Le déplacement est seulement
entre u1 et u2.
y Élément « e »
v1
y u1
v1 β
u1
x Cours 8 – Méthode des éléments finis II
MEF pour les treillis 2D
On peut écrire dans le repère local :
1
u K u u f qh
2 2 R1
0
où : u
u T
u1 v1 u2 v2 f e
et qh
R2
2
1 0 1 0 0
0 0 0 0
e
K
AE
L 1 0 1 0
Rigidité seulement
0 0 0 0 en u1 et u2.
Cours 8 – Méthode des éléments finis II
MEF pour les treillis 2D
On veut exprimer les déplacements {u} du repère
local (x,y) dans le repère global (x , y ).
v2 x
v2 u2
u2 u1 u1 cos v1 sin
y v1 u1 sin v1 cos
v1 u1 u2 u2 cos v2 sin
v1 β
v2 u2 sin v2 cos
y u1
x
Cours 8 – Méthode des éléments finis II
MEF pour les treillis 2D
On en déduit l’équation suivante :
u G u
où [G] est la matrice de transformation :
cos sin 0 0
sin cos 0 0
G
0 0 cos sin
0 0 sin cos
De la même façon, on peut écrire :
T
u G
u u G u
T T T
u G
T
Cours 8 – Méthode des éléments finis II
MEF pour les treillis 2D
L’équation d’énergie dans le repère global s’écrit
donc :
1
2
e
u G K G u u G f
T T e
et on pose : K G K G
e T e
f G f
e e
Cours 8 – Méthode des éléments finis II
MEF pour les treillis 2D
Donc :
1
e
u K u u f
2
e
u K u f 0
e e
Selon les conditions
de stationnarité
K u f
(p. 2-54)
e e
Cours 8 – Méthode des éléments finis II
MEF pour les treillis 2D
La matrice de raideur dans le système global
devient :
K G K e G
e T
T
c s 0 0 1 0 1 0 c s 0 0
s c 0 0 0 0 s c 0 0
e
K
AE
L 0 0 c s 1
0 0
0 1 0 0 0 c s
0 0 s c 0 0 0 0 0 0 s c
Rappel : Associativité
A x (B x C) = (A x B) x C
Cours 8 – Méthode des éléments finis II
MEF pour les treillis 2D
Une fois multipliée :
c2 cs c 2 cs
2
e
K
AE cs s 2
cs s
L c 2 cs c 2 cs
2 2
cs s cs s
Valide que pour les treillis (traction/compression)
De la même façon, c(qh 2 R1 )
s (qh 2 R )
f c(qh 2 R )
e
on obtient { } :
f e 1
2
s (qh 2 R2 )
Cours 8 – Méthode des éléments finis II
MEF pour les treillis 2D
Exemple 1 p. 5-5
Faire l’assemblage des matrices
Dans COMSOL …
Confirmer à l’aide de Castigliano (voir cours 3).
3 P
y c b
1 a 2
x
Cours 8 – Méthode des éléments finis II
MEF pour les treillis 2D
3
1. Calcul élémentaire P
a. Élément a: 0; s 0; c 1 y c b
T
c s 0 0 1 0 1 0 c s 0 0 u 1
AE s c 0 0 0 0 0 0 s c 0 0 v1 a 2
. 1 x
L 0 0 c s 1 0 1 0 0 0 c s u 2
0 0 s c 0 0 0 0 0 0 s c v 2
1 0 1 0 0 0 u1
0
1 0 1 0 u 1
0 0 0 0 0 v1
AE 0 0 0 0 v1 AE 1 0 1 0 0 0 u 2
.
L 1 0 1 0 u 2 .
L 0 0 0 0 0 0 v 2
0 0 0 0 v 2 0 0 0 0 0 0 u3
K
0 0 0 0 0 0 v3 a
Cours 8 – Méthode des éléments finis II
MEF pour les treillis 2D
3
1. Calcul élémentaire P
b. Élément a: 120o ; s 0.86; c 0.5 y
c b
T
c s 0 0 1 0 1 0 c s 0 0 u 2
AE s c 0 0 0 0 0 0 s c 0 0 v 2 a 2
. 1 x
L 0 0 c s 1 0 1 0 0 0 c s u 3
0 0 s c 0 0 0 0 0 0 s c v 3
0 0 0 0 0 0 u1
0
0 0 0 0 0 v1
0.25 0.43 0.25 0.43 u 2
AE 0 0 0.25 0.43 0.25 0.43 u 2
AE 0.43 0.75 0.43 0.75 v 2 .
. L 0 0 0.43 0.75 0.43 0.75 v 2
L 0.25 0.43 0.25 0.43 u 3
0 0 0.25 0.43 0.25 0.43 u 3
0.43 0.75 0.75 0.75 v 3
0 0 0.43 0.75 0.75 0.75 v 3
K b
Cours 8 – Méthode des éléments finis II
MEF pour les treillis 2D
3
1. Calcul élémentaire P
c. Élément c: 60o ; s 0.86; c 0.5 y c b
T
c s 0 0 1 0 1 0 c s 0 0 u 1
AE s c 0 0 0 0 0 0 s c 0 0 v1 a 2
. 1 x
L 0 0 c s 1 0 1 0 0 0 c s u 3
0 0 s c 0 0 0 0 0 0 s c v 3
0.25 0.43 0 0 0.25 0.43 u1
0.43 0.75
0 0 0.43 0.75 v1
0.25 0.43 0.25 0.43 u1
AE 0 0 0 0 0 0 u2
AE 0.43 0.75 0.43 0.75 v1 .
. L 0 0 0 0 0 0 v 2
L 0.25 0.43 0.25 0.43 u 3
0.25 0.43 0 0 0.25 0.43 u 3
0.43 0.75 0.75 0.75 v 3
0.43 0.75 0 0 0.75 0.75 v 3
K c
Cours 8 – Méthode des éléments finis II
MEF pour les treillis 2D
3
2. Assemblage P
u1 R1x y c b
v1 R1 y
a b c u2
K K K .
v
2y
0
R
1 a 2
2
u3 P x
v 3 0
1.25 0.43 1 0 0.25 0.43 u1 R1x
0.43 0.75 0 0 0 .43 0 .75 R
v1 1 y
1 0 1.25 0.43 0.25 0.43 u 2 0
.
0 0 0.43 0.75 0.43 0 .75 v
2y
2 R
0.25 0.43 0.25 0.43 0.5 0 u3 P
0.43 0.75 0.43 0.75 0 1.5 v 3 0
Cours 8 – Méthode des éléments finis II
MEF pour les treillis 2D
3
3. Imposition des C.L. P
y c b
u1 0 v1 0 v 2 0
1 a 2
x
1.25 0.43 1 0 0.25 0.43 u1 R1x
0.43 0.75 0 0 0 .43 0 .75 R
v1 1 y
1 0 1.25 0.43 0.25 0.43 u 2 0
.
0 0 0.43 0.75 0.43 0 .75 v
2y
2 R
0.25 0.43 0.25 0.43 0.5 0 u3 P
0.43 0.75 0.43 0.75 0 1.5 v 3 0
Cours 8 – Méthode des éléments finis II
MEF pour les treillis 2D
3
3. Imposition des C.L. P
y c b
1 a 2
x
4. Imposition des Résolution
Cours 8 – Méthode des éléments finis II
MEF pour les treillis 2D
3
5. Variables dérivées P
y c b
1 a 2
x
1.25 0.43 1 0 0.25 0.43 u1 R1x
0.43 0.75 0 0 0 .43 0 .75 R
v1 1 y R1x 10000 N
1 0 1.25 0.43 0.25 0.43 u 2 0
. R1 y 8660.3 N
0 0 0.43 0.75 0.43 0 .75 v
2y
2 R R2 y 8660.3N
0.25 0.43 0.25 0.43 0.5 0 u3 P
0.43 0.75 0.43 0.75 0 1.5 v 3 0
Cours 8 – Méthode des éléments finis II
MEF pour les poutres 1D et 2D
Sur une poutre, il peut y avoir de la flexion pure
ainsi que des efforts transverses qui eux, engendre
des efforts tranchants.
On peut négliger l’effort tranchant lorsque :
h
0.10
L
où h est la hauteur de la poutre
L est la longueur de la poutre
C’est l’hypothèse d’Euler-Bernouilli
Cours 8 – Méthode des éléments finis II
Poutre 1D
On utilise la même démarche vue au chapitre 4, p.4-3 à
4-14 et on obtient, pour des efforts en flexion
seulement :
h3
EI e
K u
f e
Cours 8 – Méthode des éléments finis II
Poutre 1D
où la matrice élémentaire de rigidité, des forces et des
déplacements sont données par :
12 6h 12 6h
2
2
e EI 6h 4h 6h 2h
K 3
h 12 6h 12 6h
2 2
6h 2h 6h 4h
6h P1 v1
h2 M
f
e qo
1
12 6h P2
1
u
h 2 M 2 v2
2
Cours 8 – Méthode des éléments finis II
Poutre 1D
Exemple 1 de la plage 6-7
12 6h 12 6h v1 6h P1
6h 4h 2 6h 2h 2 h2 M
EI 1 q 1
h 3 12 6h 12 6h v2 12 6h P2
2 2
6h 2h
2
6h 4h 2 h M2
Cours 8 – Méthode des éléments finis II
Poutre 1D
Cours 8 – Méthode des éléments finis II
Poutre 2D
Exemple de la page 6-9, poutre en 2D
Dans ce problème, il y a de la compression et de la
flexion ; les deux matrices peuvent être combinées :
0 0 0 0 0 0 1 0 0 1 0 0
0 12 6h 0 12 6h 0 0 0 0 0 0
EI 0 6h 4h 0 6h 2h EA 0 0 0 0 0 0
2 2
e
K 3
h 0 0 0 0 0
0 h 1 0 0 1 0 0
0 12 6h 0 12 6h 0 0 0 0 0 0
2
0 6h 2h 0 6h 4h 2 0 0 0 0 0 0
Ah 2 0 0 Ah 2 0 0
0 12 I 6 Ih 0 12 I 6 Ih
E 0 6 Ih 4 Ih 2 0 6 Ih 2 Ih 2
e
K 3
h Ah 2
0 Ah 2 Ah 2 0 0
0 12 I 6 Ih 0 12 I 6 Ih
0 6 Ih 2 Ih 2 0 2
6 Ih 4 Ih Cours 8 – Méthode des éléments finis II
