Dpartement de gnie civil

Laboratoire de mcanique des sols

COLE POLYTECHNIQUE
FDRALE DE LAUSANNE

Cycle postgrade:
Gologie Applique l'Ingnierie et
l'Environnement
Michel Dysli

B3-5: Introduction aux

lments finis
Documents distribus

2me dition, juin 1997

Liste des
Cycle postgrade
: Gologie Applique l'Ingnierie et l'Environnement
documents
Module B3-5: Introduction aux lments finis
M. Dysli
mise jour: 1997-06-10
Document
No
1
2
2a
3
4
4a
5
5a
6
7
8
8a
9
9a
10
11
11a
12
13
14
15
16
Annexe

Titre
Buts du cours
Bibliographie
Scalaire, vecteur, tenseur
Equations de contraintes et dformations (cours mcanique des sols B2-2)
Exemple d'un tenseur: le tenseur des contraintes
Les quations et leur solution
Solution quations de contraintes et dformations (cours mcanique des sols B2-2)
Diffrences finies
Calcul dun coulement souterrain avec un tableur
Elments finis
Autres lments finis
Exemples de discrtisation en lments finis
Exemples de discrtisation en lments finis (suite)
Fonctions de transformation et d'interpolation
Fonctions de transformation et d'interpolation (suite)
Matrices d'lasticit
Relations entre les noeuds et l'intrieur de l'lment (1re partie)
Relations entre les noeuds et l'intrieur de l'lment (suite)
Dtermination de la matrice dformation - dplacement
Application numrique
Application numrique (suite)
Principe de la rsolution par la mthode des dplacements
FEM, BEM, KEM
Rappel du calcul matriciel

Cycle postgrade
: Gologie Applique l'Ingnierie et l'Environnement
Module B3-5: Introduction aux lments finis
M. Dysli

Buts
du cours

Buts du cours "Introduction aux lments finis" (8 heures):

Comprendre les principes de la mthode

des lments finis applique au calcul
des contraintes et dformations et
celui des coulements souterrains.
Avoir les bases ncessaires pour s'en
servir dans des cas simples.
Et non pas :
Etre capable d'crire des programmes sur ordinateur relatifs
cette mthode... Il faudrait pour cela un cours de 100 200
heures!

Mthode: approche physique par des cas particuliers

Cycle postgrade
: Gologie Applique l'Ingnierie et l'Environnement
Module B3-5: Introduction aux lments finis
M. Dysli

Document
No 1
ajout 1997-06-24

Bibliographie
1. Acquisitions conseilles
BATHE K.-J. (1996). Finite element procedures . Prentice-Hall.
BRITTO A. M., GUNN M. J. (1987). Critical State Soil Mechanics via Finite Elements. Ellis
Horwood Limited, Halsted Press: a division of John Wiley & Sons,
DHATT G., TOUZOT G. (1984). Une prsentation de la mthode des lments finis. Maloine SA
Paris.
PRAT M. (1995). La modlisation des ouvrages, Herms, Paris.
ZIENKIEWICZ O. C. (1979). La mthode des lments finis. 3e d., McGraw Hill.
ZIENKIEWICZ O. C. (1977). The Finite Element Method. 3rd ed., McGraw Hill.

2. Autres
CHEN W. F., MIZUNO E. (1990). Nonlinear Analysis in Soil Mechanics. Elsevier
DESAI C. S., SIRIWARDANE H. J. (1984). Constitutive laws for engineering materials. PrenticeHall.

Des ditions plus rcentes pourraient exister.

Document
No 2

Cycle postgrade
: Gologie Applique l'Ingnierie et l'Environnement
Module B3-5: Introduction aux lments finis
M. Dysli

Scalaire, vecteur, tenseur

En physique, lorsqu'un simple nombre suffit dfinir une grandeur (temprature, densit, module
lastique, par ex.), ce nombre est un scalaire.
Si la grandeur a non seulement une dimension mais aussi une direction, on l'appellera alors vecteur.
Une force, une acclration, un dplacement sont des vecteurs. Dans un systme d'axes cartsien 3
dimensions, il faut 3 composantes au moins pour dfinir un vecteur, soit par exemple deux angles et
une intensit, ou 3 projections sur les axes du systme choisi.
Une modification du systme de rfrence induit une modification des 3 composantes dfinissant le
vecteur. Il est cependant possible de driver de ce vecteur une quantit scalaire qui est indpendante
de sa direction. Par exemple, si l'on considre un vecteur dplacement de composantes d1, d 2, d3, la
grandeur de ce vecteur: d = d12 + d22 + d32 est indpendante de sa direction; c'est l'invariant du
vecteur.
La notion de tenseur est plus abstraite. On pourrait lui donner la signification physique d'une reprsentation d'un champ: champ de vitesses, champ de contraintes par exemple. Dans la dfinition des
composantes d'un tel champ, nous devons nous rfrer deux fois des directions, par exemple, pour
un champ de contrainte: premirement l'orientation du solide lmentaire et deuximement
l'orientation des contraintes proprement dites sur les faces de ce solide. Il faudrait donc 3 . 3 = 9
composantes pour dfinir un tel tenseur. Cependant un tenseur cartsien est symtrique et seules 6
composantes suffisent le dfinir compltement.
Si pour un vecteur il est possible de driver un invariant, un tenseur donne trois invariants dans un
systme 3 dimensions.
Type

Scalaire

Vecteur

Tenseur

poids
volumique

dplacement

contrainte

Notations

d1

d2
d
3

Notations indicielles

di

ij

Nombre de
composantes dans un
systme de coordonnes 3 dimensions

30

31 = 3

32 = 9

Valeurs
indpendantes

9 en gnral

Invariants

Ordre de la matrice
Exemple

11 12

21 22
31 32

13
23
33

6 avec symtrie: ij = ji
1

Document
No 02a

Cycle postgrade
: Gologie Applique l'Ingnierie et l'Environnement
Module B3-5: Introduction aux lments finis
M. Dysli

Les quations des contraintes et dformations

Notation indicielle (d'Einstein)

Nbre
Nbre
quations inconnues

Notation xyz

'xx yz zx u
+
+
+ + Fx = 0
x
y
z
x

1. Equations d'quilibre
i = 1,2,3 = direction
'ij u
des normales
+
+ Fj = 0
aux faces
xi xi
j = 1,2,3 = direction
des axes
= ' + u
u = pression interstitielle
Fj = force par unit de volume
En gnral, F x (ou F 1) = F y (ou F 2) = 0
sauf si forces d'inertie (sisme par ex.)
Fz (ou F3) = -.g car sol = milieu pesant

qui s'crit habituellement :

3
z ou 3

'x yz zx u
z
+
+
+ + Fx = 0
x
y
z x
zx zy
xy 'y zy u
yz
+
+ + Fy = 0
+
y
xz
x
y
z y
x xy yx y ou 2
zx yz 'z u
+
+ + Fz = 0
+
x
y
z z
x ou 1

2. Relations dformation-dplacement (quations gomtriques)

ux ux ux u
x = 1
+
=
=
ui u j
i = 1,2,3
2
ij = 1
+
x
x
x x
j
=
1,2,3
2 x j xi
ux = u ; uy = v ; uz = w
ui,j = vecteur des dplacements
w
v
z =
y =
z
Formulation de Lagrange : la position d'un
y
point est dcrite en fonction de ses coordon
u v
u v
nes initiales (avant dformation).
xy = xy = 1
+
xy =
+
y x
2 2 y x
v w
Par opposition la formulation d'Euler o la
yz =
+
z y
position d'un point est dcrite en fonction de
ses coordonnes aprs dformation.
w u
+
zx =
x z

x
y

ij = z
xy
yz
zx

i = direction des ,
j = direction des

loi linaire
loi non linaire

sij est la matrice du matriau et son inverse

sij-1 la matrice d'lasticit . sij est symtrique,
elle a donc, pour le cas le plus gnral d'un
matriau anisotrope : 36 - 15 = 21 termes.
Exemples (lastiques linaires) :

fil

x = s 11 x = 1/E x
z

z
x = s31 3 = -/E z
y = s32 3 = -/E z
z = s33 3 = 1/E z
z

x
z

voir 2.2.3 du
cours polycopi

x
y

ij = z
xy
yz
zx

sij =

15

il manque 6 quations

(lois constitutives, quations d'tat, quations de dformabilit)

R(ij, ij) = 0 ou ij = f(ij) = sij . ij

3
3
ij = ji

u, v, w
x , y ,
z
xy, yz,
zx

Total intermdiaire

3. Relations contrainte-dformation

sij = cte =>

sij = f(, , t) =>

15

15

s11 s12 s13 s14 s15 s16

s21
s31
etc.
s41
s51
s61

Dans le cas de la mcanique des sols,

les des exemples seraient ngatifs
(compression = +)

Total final

Document
No 3

Cycle postgrade
: Gologie Applique l'Ingnierie et l'Environnement
Module B3-5: Introduction aux lments finis
M. Dysli

Exemple d'un tenseur: le tenseur des contraintes

1, 2, 3 = axes des contraintes
principales
= 0 sur les faces
z

z
z
zx
xz
x

xy

zy
yz
y

'
z' z

z'x'

x'z'

yx

x'

z
3
3

z'y'

y'
'
y'z
y'

x'y'

2
y

y'x'

x'
1

tenseur des contraintes :

x
yx
zx

xy
y
zy

tenseur des contraintes

principales :
1 0
0
2
0
0
3
0
0

xz
yz
z

xy = yx
avec: xz = zx
zy = yz

Invariants du tenseur des contraintes (sans dmonstration)

I1 = x + y + z

I1 = 1 + 2 + 3

I2 = x y + y z + z x xy2 yz2 zx2

I2 = 1 2 + 2 3 + 3 1
I3 = 1 2 3

I3 = x y z + 2 xy yz zx x yz2 y zx2 z xy2

' = - u
x' xy
'
' y'
yx
zx
' zy
'

xz'
x
yz' = yx
z'
zx

xy
y
zy

' = contrainte effective

= contrainte totale
u

xz
u 0 0
yz - 0 u 0
z
0 0 u

pression interstitielle =
tenseur sphrique

= pression interstitielle

1'
0
0

0
2'
0

1
0
0 = 0
3'
0

0
2
0

0
u 0 0
0 - 0 u 0
3
0 0 u

Cycle postgrade
: Gologie Applique l'Ingnierie et l'Environnement
Module B3-5: Introduction aux lments finis
M. Dysli

Document
No 4

Les quations et leur solution

Ecoulement souterrain
1 h
h = c
v t
avec

h = charge hydraulique
t = temps
cv = coefficient de

Contraintes-dformations
' u
+
+ Fx = 0
quilibre des forces
x x
avec

consolidation

'=
u=
Fx =
x=

Conduction thermique
Champs lectriques
etc.

contraintes effectives
pression interstitielle
force volumique
axe cartsien quelconque

vi vj relation dformation +
2 xj xi dplacement
avec vi,j = vecteur des dplacements

ij = 1

ij =

= f( ')

dformations

relation contraintes - dformations

..
M .v + K .v = R - F + U
avec

M=
K=
v=
R=
F=
U=

quilibre des dplacements

(document No 15)

matrice des masses.

matrice de rigidit.
vecteur des dplacements.
vecteur des forces extrieures.
vecteur des forces intrieures
vecteur des pressions interstitielles.

Cas gomtriquement simple : solutions analytiques

Cas plus complexe: diffrences finies ou lments finis

Excavation avec paroi moule
WS

puits
filtrant

couche 1

WS

couche 2

Diffrences finies : document No 5

Elments finis : document No 6

Cycle postgrade
: Gologie Applique l'Ingnierie et l'Environnement
Module B3-5: Introduction aux lments finis
M. Dysli

Document
No 4a

Solution des quations de contraintes et dformations

Charges simples pouvant conduire une solution analytique :

Semi-infini
lastique

Solutions

Abaques

analytiques

et

formules

Recordon Ed. (1980). Abaques du cours polycopi de mcanique des sols de l'EPFL.
Giroud J.-P. (1975). Tables pour le calcul des fondations, Tomes 1 et 2. Dunod, Paris.
Poulos H. G., Davis E. H. (1974). Elastic solutions for soil and rock mechanics. John Wyley & Sons.

Cycle postgrade
: Gologie Applique l'Ingnierie et l'Environnement
Module B3-5: Introduction aux lments finis
M. Dysli

Document
No 5

Diffrences finies
Ecoulement souterrain

Contraintes-dformations
WS

puits
filtrant

couche 1

WS

couche 2

y
Fi+1,j - Fi-1,j
F
=
x
2

j+2

Fi,j+1 - Fi,j-1
F
=
y
2

j+1

2 F
x2
2 F
y2

i-2

i-1

i+1

i+2

Fi+1,j - 2Fi,j + Fi-1,j

22

j-1

Fi,j+1 - 2Fi,j + Fi,j-1

22

j-2

F = contraintes, dformations, dplacements, tempratures, pressions, moments, charge hydraulique, etc.

Exemple: coulement souterrain
z

palplanche
h = 24

24

h0 = h2/4 + h3/2 + h 4/4

22
h = 18

h = 24

h6

savoir:
hi+1,j - 2hi,j + hi-1,j

h = 18

h5

2
2
Equation pour rgime permanent: h + h = 0
2
x
y2
h = charge hydraulique = z + u / w avec u = pression interstitielle

hi,j =

hi,j+1 - 2hi,j + hi,j-1

2 2

hi-1,j + hi,j-1 + hi+1,j + h i,j+1

4

=0

ou h 0 =

h1 + h 2 + h 3 + h 4
4
h4

Cas particuliers aux limites :

h1

h0

impermable
h2

h0 = h1/4 + h3/4 + h4/2

h0 = h1/2 + h2/4 + h4/4

h0 = h 1/4 + h 2/4 + h3/4 + h 5/8 + h 6/8

Excellent moyen = tableur

h3

Calcul d'un coulement souterrain avec un tableur

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18

z
24
22
20
18
16
14
12
10
8
6
4
2
0

0
24
24
24
24
24
24
24
24
24
24
24
24
24

2
24
23.92
23.83
23.74
23.66
23.59
23.54
23.50
23.48
23.47
23.46
23.46
23.46

4
24
23.83
23.66
23.48
23.31
23.16
23.05
22.98
22.94
22.92
22.91
22.91
22.91

6
24
23.75
23.49
23.21
22.94
22.70
22.53
22.44
22.39
22.37
22.36
22.36
22.36

8
24
23.67
23.32
22.94
22.53
22.15
21.94
21.85
21.82
21.81
21.81
21.82
21.82

10
24
23.61
23.19
22.70
22.10
21.42
21.23
21.20
21.22
21.24
21.26
21.28
21.29

=(E9+F10+G9+F8)/4

H
x
12
24
23.59
23.14
22.57
21.75
20.21
20.37
20.51
20.61
20.68
20.72
20.75
20.76

14

16

18

20

22

24

18
18.36
18.81
19.53
19.86
20.04
20.14
20.20
20.23
20.24

18
18.32
18.68
19.08
19.36
19.54
19.65
19.71
19.74
19.75

18
18.25
18.51
18.76
18.96
19.10
19.20
19.25
19.28
19.29

18
18.17
18.33
18.49
18.62
18.71
18.78
18.82
18.84
18.85

18
18.08
18.16
18.24
18.30
18.35
18.39
18.41
18.42
18.42

18
18
18
18
18
18
18
18
18
18

Cycle postgrade
: Gologie Applique l'Ingnierie et l'Environnement
Module B3-5: Introduction aux lments finis
M. Dysli

Document
No 6

Elments finis
Ecoulement souterrain

Contraintes-dformations
WS

puits
filtrant

couche 1

WS

couche 2

WS

uy,Fy
hm

hj
hn

hi

N qM-N h
p
ho

hk

hi = charge hydraulique
qi = dbit

j
n

o
hl

qO

M
qM

qO-N

Elments
rectangulaires

ux,Fx

x
l

O
ui = dplacements
Fi = forces

Document
No 7

Cycle postgrade
: Gologie Applique l'Ingnierie et l'Environnement
Module B3-5: Introduction aux lments finis
M. Dysli

ajout 1997-06-24

Autres lments finis

exemples parmi de nombreux autres

barres

seul l'effort de traction / compression est transmis

poutre

une unit

2D

Efforts dans le plan

3D

joint
cinmatique

corps 1

Elment de
cisaillement

= c + tg
Intersection ?

joints

corps 2

Document
No 8

Cycle postgrade
: Gologie Applique l'Ingnierie et l'Environnement
Module B3-5: Introduction aux lments finis
M. Dysli

Exemples de discrtisation en lments finis

z

paroi
moule

axe de symtrie

couche 1

seulement dplacement
vertical

couche 2

2D
couche 3

Rseau initial
couche 4

y
aucun dplacement possible

naissance du buton
aprs tape excavation 2
e1

couche 1

couche 2

e = tapes d'excavation

e2
e4
e5

e3

Rseau aprs plusieurs

pas d'excavation

couche 3

couche 4

De par la symtrie seul 1/4

de l'excavation est modlise

3D
Elment
isoparamtrique
21 noeuds

Cycle postgrade
: Gologie Applique l'Ingnierie et l'Environnement
Module B3-5: Introduction aux lments finis
M. Dysli

Document
No 8a

Exemples de discrtisation en lments finis (suite)

Cycle postgrade
: Gologie Applique l'Ingnierie et l'Environnement
Module B3-5: Introduction aux lments finis
M. Dysli

Document
No 9
modification 1997-06-24

Fonctions de transformation et d'interpolation (1re partie)

Approche trs simple n'utilisant que les forces, les dplacements, les contraintes et les
dformations. Les mmes principes s'appliquent par analogie aux coulements souterrains et
tous les autres champs.
1) Compatibilit des dplacements entre les cts
de deux lments voisins

mauvais

2) Fonction de transformation gomtrique:

Transformation des coordonnes relles en coordonnes de rfrence
Exemple pour un lment triangulaire

uy

ux
xk

xj
xi

Elment rel
Coordonnes locales

h() = fonction de transformation gomtrique

Par exemple, pour l'lment triangulaire :

0, 1

x( , ) = h 1 ( , ) xi + h 2 ( , ) x j + h 3 ( , ) xk
1, 0

0, 0
Elment de rfrence
Coordonnes naturelles

xi

x( , ) = {h } xj

xk

De la mme
manire :
yi

y( , ) = {h } yj

yk

Document
No 9a

Cycle postgrade
: Gologie Applique l'Ingnierie et l'Environnement
Module B3-5: Introduction aux lments finis
M. Dysli

modification 1997-06-24

Fonctions de transformation et d'interpolation (suite)

Approche trs simple n'utilisant que les forces, les dplacements, les contraintes et les
dformations. Les mmes principes s'appliquent par analogie aux coulements souterrains et
tous les autres champs.
3) Fonctions d'interpolation:
Passage des valeurs aux noeuds (variables nodales) celles en tout point d'un lment
Exemple pour un lment triangulaire
ui,j,k = variables nodales = variables attaches aux noeuds de l'lment rel comme des
dplacements, des forces ou une pression interstitielle.
Passage des valeurs aux noeuds (variables nodales u) celles en tout point x,y ou , du triangle : u
i
ui

u(
,
)
=
h
(
,
)
h
(
,
)
h
(
,
)

{1
}uj
2
3


u( x, y) = {h1( x, y) h 2 ( x, y) h 3 ( x, y)} u j

h(,) = fonction d'interpolation sur uk
u
k
l'lment de rfrence
h(x) = fonction d'interpolation sur l'lment rel
Pour le triangle (interpolation linaire) :
Pour le triangle (interpolation linaire) :
1
h 1 ( , ) = 1
h 1 ( x, y) =
(y y j )(x j x) (xk x j )(y j y)
2A k
1
h 2 ( x, y) =
h 2 ( , ) =
[( y yk )( xk x) ( xi xk ) (yk y)]
2A i
1
h 3 ( x, y) =
h 3( , ) =
(y y )( xi x) (x j xi )( yi y)
2A j i

2 A = (xk x j)(yi y j ) (xi x j )(yk y j)

Elment:
isoparamtrique
sub-paramtrique
super-paramtrique

h() = h()
ordre h() < ordre h()
ordre h() > ordre h()

4) Jacobien
Toutes les expressions des fonctions d'interpolation qui impliquent des drives de u (variable nodale) en
x, y sont transformes en drives de et grace la matrice de transformation dite matrice jacobienne
ou Jacobien J :
u = quelque chose
u x

=
u x

y u
x

y u
y

{}= [ J ] { x}

{ x} = [ J] 1 {}
Problme lors de l'inversion si, par
exemple, le dterminant de J est nulle

Document
No 10

Cycle postgrade
: Gologie Applique l'Ingnierie et l'Environnement
Module B3-5: Introduction aux lments finis
M. Dysli

ajout 1997-06-24

Matrices d'lasticit
Forme gnrale de la relation - pour un corps lastique:
Premire solution: dcomposition du tenseur des contraintes en en tenseur sphrique m et en un
tenseur dviatorique o. Les tenseurs correspondant des dformations sont m et o.
L'indice o de et est la lettre o et non pas le chiffre 0.
lasticit linaire:
= 3K .
m

lasticit non linaire :

= 3K .

o = 2G . o

avec :
K = coefficient de compressibilit
[kPa par ex.]
G = module de glissement ou de
cisaillement [kPa[

o = 2G . o

Deuxime solution: usage de vecteurs de contraintes et de dformations

x
y
z
xy
yz
zx

x
y
z
=
xy
yz
zx

ij = sij . ij

s11 s12 s13 s14 s15 s16

s21
s31
etc.
s41
s51
s61
et

ij = sij-1 . ij

zx
xz
x

xy

zy
yz
y
yx

Rappel: xy = xy
2xy = xy

forme la plus utilise

sij-1 est l'inverse de la matrice sij

Simplifications ( 2D )
z 0

ij =

1
E
ij =

(1- 2)
0

0
0 ij = [ D] ij
(1-)
0
D = matrice d'lasticit
2

Dformations planes (z = 0)
1-

1-
0

0
ij = [ D] ij
0
1-2
2

ij =

z
ij = r

rz

z
z

Symtrie de rvolution
z
r

rz

E
ij =
(1+) (1-2)

x
ij = y
xy

1-
E
ij =
(1+) (1-2)

0 1-
0
0
0

1-
0

0
0
ij
0
1-2
2

ij =

x
y
xy

x
ij = y
xy

x
y
xy

Contraintes planes (z = 0)

z 0

zr
rz

Cycle postgrade
: Gologie Applique l'Ingnierie et l'Environnement
Module B3-5: Introduction aux lments finis
M. Dysli

Document
No 11
modification 1997-06-24

Relations entre les noeuds et l'intrieur de l'lment (1re partie)

Approche trs simple n'utilisant que les forces, les dplacements, les contraintes et les dformations.
1)

Relation entre les dformations l'intrieur de l'lment et les

u
dplacements aux noeuds : Matrice des dformations - dplacements
x
v1
v2
v5
v
u2
u
y

(x,y)=
document
No
35
module
B22
:
1
u5

1
v

u v
8
Une matrice B diffrente pour chaque point x,y
5
v6 2

u
8
y + x
x,y
u6
8
6
(1)
{ (x,y)} = [B (x,y)] uv noeuds
v4
v3
noeuds
7
u
4
4
Dimensions (lignes . colonnes):
u3
3
v7
u, v

B
u
.6
lment triangulaire 3 noeuds
7
3
3
6
x
.

3
3

38
3.16

8
16

lment rectangulaire 4 noeuds

lment rectangulaire 8 noeuds

2)

Relation entre les contraintes l'intrieur de l'lment et les forces aux

noeuds : Matrice des forces aux noeuds - contraintes
x
Fy2
Fy1
Fy5
( x,y) = y
xy
Fx2
Fx1
Fx5
(2)
F
=
[
C(x,y)
]
(
x,
y)

{ noeuds}
{
}
1 Fy8
2
5
Fy6
Dimensions (lignes . colonnes):
x,y
Fx6
F
8
x8

F
C
6
F
y4
lment triangulaire 3 noeuds
6. 3
6
3
Fy3
7
Fx4
4
lment rectangulaire 4 noeuds
8
3
8. 3
3
y Fx3
lment rectangulaire 8 noeuds
16
3
16.3
Fx7
F
Pas utilise par la suite
y7
x

3) Points de Gauss et mthode des rsidus pondrs

Dans le cas de l'lment isoparamtrique avec h = h > 1 (par exemple interpolation cubique) on peut
utiliser, pour passer des dplacements aux noeuds aux dformations l'intrieur de l'lment, la mthode
d'intgration numrique de Gauss qui permet de calculer, avec une trs bonne approximation, les dformations et les contraintes en des points dtermins: les points de Gauss, dont les positions permettent
d'obtenir la prcision maximale.
La mthode de Gauss fait partie d'un ensemble de mthodes dnomm: mthode des rsidus pondrs.
Cette mthode gnrale d'intgration numrique qui utilise des fonctions de pondrations permet de
passer d'un systme d'quations aux drives partielles une formulation intgrale en utilisant les
fonctions de d'interpolation (polynomiale par exemple). La prcision de la formulation intgrale dpend
de la position du point considr dans l'lment. Ces formulations intgrales sont nombreuses, par
exemple:
mthode de Gauss dj cite

mthode de Ritz

= environ Gauss
mthode de Legendre

mthode de Galerkine ( la mode aujourd'hui!)

mthode de Newton-Cotes
=> diffrente des autres et moins prcises.

Dans la relation { (x,y)} = [B (x,y)] uv noeuds , la matrice B est forme de termes provenant des
noeuds
fonctions de pondrations, des fonctions d'interpolation et des fonctions de transformation
gomtrique. Un exemple de constitution de cette matrice est donn sur le document no 12.

Document
No 11a

Cycle postgrade
: Gologie Applique l'Ingnierie et l'Environnement
Module B3-5: Introduction aux lments finis
M. Dysli

modification 1997-06-24

Relations entre les noeuds et l'intrieur de l'lment (suite)

Approche trs simple n'utilisant que les forces, les dplacements, les contraintes et les dformations.

4) Relation forces - dplacements avec usage de la matrice d'lasticit D et de la matrice B

{Fnoeuds} = [K ] {uvnoeuds
}
noeuds
[K ] =

S [B ]T [ D] [B] dS

dS = det [ J] d d

(5)

[K ] = [A] d d (6)
S

[K ] =

i,j

avec [K] = matrice de rigidit

de l'lment

(3)

(4)

Dimensions de K (lignes . colonnes):

6.6
lment triangulaire 3 noeuds
lment rectangulaire 4 noeuds
8.8
16.16 lment rectangulaire 8 noeuds

sans dmonstration, intgration sur le surface de l'lment

, = coordonnes naturelles, J = Jacobien (document 9a), det = dterminant

avec [] = [] [D] [B] det [J]

i,j [A ] i,j (7) ou [A ] i,j est value chaque point i et i

[ ] i,j sont des constantes qui dpendent des valeur de i et i (intgration de Gauss)

Il reste assembler tous les lments. La mthode gnrale utilise pour cela fait l'objet du document
15. Nous allons avant dire encore quelques mots sur les mthodes d'intgration numrique, dont celle
de Gauss (document 12).

Document
No 12

Cycle postgrade
: Gologie Applique l'Ingnierie et l'Environnement
Module B3-5: Introduction aux lments finis
M. Dysli

modification 1997-06-24

Dtermination de la matrice dformation - dplacement

Deux dimension, lment 4 noeuds avec 4 points d'intgration (degr n = 2)
y, v

noeud 2

noeud 1

0, +1
4

2
4
-0,577

4
+0,577

+1, 0
dplacements
v4
u4

4
-0,577 -1, 0

3
4
-0,577

Fonctions d'interpolation:
(sans dmonstration)
Noeuds :
4
5
6
7
8
-1/2h 8
h1 = 1/4 (1 + )(1 + ) -1/2h 5
h2 = 1/4 (1 - )(1 + ) -1/2h 5 -1/2h 6
h3 = 1/4 (1 - )(1 - )
-1/2h 6 -1/2h 7
h4 = 1/4 (1 + )(1 - )
-1/2h 7 -1/2h 8
h5 = 1/2 (1 - 2)(1 + )
h6 = 1/2 (1 - )(1 - 2)
h7 = 1/2 (1 - 2)(1 - )
h8 = 1/2 (1 + )(1 - 2)
h9 = 1/2 (1 - 2)(1 - 2)

Point de Gauss

y4

noeud 4
0, -1

x4

noeud 3
ainsi:

x = 1/4 (1 +
y = 1/4 (1 +
u = 1/4 (1 +
v = 1/4 (1 +

)(1 + )x1 + 1/4 (1 )(1 + )y1 + 1/4 (1 )(1 + )u1 + 1/4 (1 )(1 + )v1 + 1/4 (1 -

x, u

z = h1z1 + h2z2 + h3z3 + ... + hizi avec z = x, y, u ou v

)(1 + )x2 + 1/4 (1 )(1 + )y2 + 1/4 (1 )(1 + )u2 + 1/4 (1 )(1 + )v2 + 1/4 (1 -

)(1 - )x3 + 1/4 (1 +

)(1 - )y3 + 1/4 (1 +
)(1 - )u3 + 1/4 (1 +
)(1 - )v3 + 1/4 (1 +

Nous savons dj que (doc. No 11 ou module B22) :

u v
x = u
y = v
xy =
+
x
y
y x

x / = 1/4 (1 + )x1 - 1/4 (1 +

x / = 1/4 (1 + )x1 + 1/4 (1 y / = 1/4 (1 + )y1 - 1/4 (1 +
savoir en drivant (1) : y / = 1/4 (1 + )y1 + 1/4 (1 u / = 1/4 (1 + )u1 - 1/4 (1 +
u / = 1/4 (1 + )u1 + 1/4 (1
v / = 1/4 (1 + )v1 - 1/4 (1 +

x
-1
v / = 1/4 (1 + )v1 + 1/4 (1
= [J ]

et que (Jacobien, document No 9a) :

x

=
x

y
x

y
y

et

x son

= [ J ]
inverse:

y

u
x 1
-1 1 + j
ainsi : = [ J ]
ij 1 + i
u 4

-(1 + j)
1 - i

0
0

0
0

-(1 - j)
-(1 - i)

0
0

)(1 - )x4
)(1 - )y4
)(1 - )u4
)(1 - )v4

1 + j
1 + i

0
0

-(1 + j)
1 - i

0
0

-(1 - j)
-(1 - i)

1 - j

0 . u

-(1 + i) 0 v

(3)

1 - j
. u
-(1 + i) v

(4)

0
0

(1)

)x2 - 1/4 (1 - )x3 + 1/4 (1 - )x4

)x2 - 1/4 (1 - )x3 - 1/4 (1 + )x4
)y2 - 1/4 (1 - )y3 + 1/4 (1 - )y4
)y2 - 1/4 (1 - )y3 - 1/4 (1 + )y4
(2)
)u2 - 1/4 (1 - )u3 + 1/4 (1 - )u4
)u2 - 1/4 (1 - )u3 - 1/4 (1 + )u4
)v2 - 1/4 (1 - )v3 + 1/4 (1 - )v4
)v2 - 1/4 (1 - )v3 - 1/4 (1 + )v4

= j

v
x 1
-1 0
= [J ]

v
4
ij 0

y = i

9
-1/4h 9
-1/4h 9
-1/4h 9
-1/4h 9
-1/4h 9
-1/4h 9
-1/4h 9
-1/4h 9

u1
v
1
u2
u v2
u = = u
v 3
v
u3
4
v4

= j

En valuant les relations (3) et (4), on peut tablir le matrice de transformation dformation - dplacement, savoir la
matrice B, en un point i, j : ij = Bij . u. Les indices i et j indiquent que la transformation est value au point i, j.
Par exemple, si les axes x et y de l'exemple 4 noeuds sont confondus avec les axes et et si l'lment est un carr de
dimension 2 de cot, le Jacobien est une matrice unitaire et ainsi:
u / x

x
1 + j
0
-(1 + j)
0
0
-(1 - j)
0
1 - j

1
ij = y = v / y
d'o : B =
0
1 + i
0
1 - i
0
-(1 + i)
0
-(1 + i)
4

xy
1
+

-(1
+

)
)
-(1
)
-(1
)
-(1
+
1
+

i
j
i
j
i
j
i 1 - j
u / y + v / x

On peut constater que les valeurs de la matrice B dpendent des coordonnes naturelles i et j

Document
No 13

Cycle postgrade
: Gologie Applique l'Ingnierie et l'Environnement
Module B3-5: Introduction aux lments finis
M. Dysli

Application numrique
y, v
noeud 2

noeud 1

4m

En appliquant les relations (1) du document No 12:

x = 3 et y = 2

x, u

x

Jacobien = J=
x

y 3 2
3 0
=
=
y 3 2 0 2

Dans ce cas particulier J est indpendant des et

noeud 3

noeud 4
6m

0,333 0
Son inverse, J-1 =
0 0,5

Avec les relations (3) et (4) du document No 12, on peut maintenant dterminer la matrice B,
par exemple pour le point de Gauss No 1 (i = -0,5774, j = -0,5774) :
u
x 1 0,333 0 1 + j 0 -(1 + j) 0 -(1 - j) 0 1 - j 0
=

0 0,5 1 + i 0 1 - i 0 -(1 - i) 0 -(1 + i) 0

u 4
y

u 1
=
v 4

v
x 1 0,333 0 0 1 + j 0 -(1 + j) 0 -(1 - j) 0 1 - j
=

v 4 0 0,5 0 1 + i 0 1 - i 0 -(1 - i) 0 -(1 + i)

y

u 1
=
v 4

et ainsi B1 = 1
4

0.14 0 -0.14
0.21 0 0.79

0
0

-0.53
-0.79

0
0

0.53
-0.21

0 u
0 v

+
0
0

0.14
0.21

0
0

-0.14
0.79

0 -0.53
0 -0.79

0 0.53 u
0 -0.21 v

0.14
0 -0.14
0 -0.53
0
0.53
0
0
0.21
0
0.79
0
-0.79 0 -0.21
0.21 0.14 0.79 -0.14 -0.79 -0.53 -0.21 0.53

Calculons enfin les dformations au point de Gauss No 1 (1) pour le dplacement des noeuds suivant
(cas trivial pour contrle) :
noeud 2

noeud 1
0,5 m

4m
1
noeud 3

noeud 4
6m

u1
0
v
0.5
1

u2
0
x
0.14
0
-0.14
0
-0.53
0
0.53
0

1 = y = B uv2 = 14 0 0.21 0 0.79 0 -0.79 0 -0.21 0.5

0

3

0.21
0.14
0.79
-0.14
-0.79
-0.53
-0.21
0.53
v3
0
xy
u
0
4

v4
0

x
0

1 = y = -0.125

xy
0

c.q.f.d. !
Suite sur document No 14

Document
No 14

Cycle postgrade
: Gologie Applique l'Ingnierie et l'Environnement
Module B3-5: Introduction aux lments finis
M. Dysli

ajout 1997-06-24

Application numrique (suite)

Il reste maintenant constituer la matrice de rigidit K de l'lment au moyen de la relation (7) du

document No 11. Les coefficients sont donns sur le tableau ci-dessous (sans dmonstration).

[K] =

i,j

i,j [A ] i,j = 1 [A] 1 + 2 [A] 2 + 3 [A] 3 + 4 [A] 4

Nbre de pts
de Gauss

Dans notre cas 1 = 2 = 3 = 4 = 1

2
4
6

No des points de Gauss

Pour calculer A, il faut connatre la matrice d'lasticit D.

En admettant E = 20'000 et = 0,3 et des contraintes planes (doc. 10) :
21'978 6'593
0
6'593 21'978 0
0
0
7'692

[D] =

292.3
159.5
317.2
188.8
159.5
425.2
246.9
1316.3
317.2
246.9
1957.9 -595.2
188.8 1316.3 -595.2 5184.0
-1091.2 -595.2 -1183.9 -704.9
-595.2 -1587.2 -921.5 -4913.1
481.7
188.8 -1091.2 1111.2
246.9
-154.3 1269.8 -1587.2

[A] 1 = [B]1T [D] [B]1 =

2407.4
595.2
-1797.8
704.9
-1091.2
-1111.2
481.7
-188.8

[A] 2 =
etc.

595.2
1165.5
921.5
576.0
-1269.8
-1587.2
-246.9
-154.3

Fx1
F
y1
Fx2
F
y2
Fx3
Fy3
F
x4
Fy4

k11

k12

8730.5
3571.4
3571.4 12699.0
-2961.3
274.7
-274.7
3784.6
-4364.7 -3571.4
-3571.4 -6348.7
-1404.5
-274.7
274.7 -10134.9

-2961.3
274.7
8730.5
-3571.4
-1404.5
-274.7
-4364.7
3571.4

704.9
576.0
-2221.6
5924.3
921.5
-4913.1
595.2
-1587.2

0,0
0,5774
0,7746
0,0000
0,8611
0,3400

2.0000
1.0000
0,5555
0,8888
0,3479
0,6521

-1091.2 -595.2 481.7

246.9
-595.2 -1587.2 188.8
-154.3
-1183.9 -921.5 -1091.2 1269.8
-704.9 -4913.1 1111.2 -1587.2
4072.9 2221.6 -1797.8 -921.5
2221.6 5924.3 -704.9
576.0
-1797.8 -704.9 2407.4 -595.2
-921.5
576.0 -595.2 1165.5
-1091.2 -1111.2 481.7
-188.8
-1269.8 -1587.2 -246.9 -154.3
-1183.9 704.9 -1091.2 595.2
921.5 -4913.1 595.2 -1587.2
1957.9
595.2
317.2
-246.9
595.2
5184.0 -188.8 1316.3
317.2
-188.8 292.3
-159.5
-246.9 1316.3 -159.5
425.2

{Fnoeuds} = [K] {uv noeuds

}
noeuds

[K] = 1.[A] 1 + 1.[A] 2 + 1.[A] 3 + 1.[A] 4

k21

-1797.8
921.5
4072.9
-2221.6
-1183.9
704.9
-1091.2
595.2

[K] =

matrice de rigidit
d'un lment

symtrique kij = kji

-274.7
3784.6
-3571.4
12699.0
274.7
-10134.9
3571.4
-6348.7

-4364.7 -3571.4 -1404.5 274.7

-3571.4 -6348.7 -274.7 -10134.9
-1404.5
-274.7 -4364.7 3571.4
274.7 -10134.9 3571.4 -6348.7
8730.5
3571.4 -2961.3 -274.7
3571.4 12699.0 274.7
3784.6
-2961.3
274.7
8730.5 -3571.4
-274.7
3784.6 -3571.4 12699.0

u1
v
1
u2
v
2
u3
v3
u
4
v4

F
x1
Fy1

i = no noeud o l'on
dsire connatre les
forces

Fy2

v1

v2

v3

v4

connus

= k11 1 + k12 2 + k13 3 + k14 4

v1
v2
v3
v4 etc. Les sous-matrices k sont
ij
j = no noeud o les
F
u1
u2
u3
u4
aussi
elles-mmes
symtriques
x2
dplacements sont

= k21 + k22 + k23 + k24

Appuis :
noeud 2

c1, c 2, c 3, c 4 = constantes

noeud 1

k14 = k41 =

c1 0
0 0

k24 = k42 =

c2 0
0 0

k34 = k43 =

c3 0
0 0

k44 =

c4 0
0 0

noeud 4
noeud 4
noeud 3

v4 = 0 d'o :

u4 = v4 = 0 d'o : k14 = k41 = k24 = k42 = k34 = k43 = k44 = 0 0

0 0

Document
No 15

Cycle postgrade
: Gologie Applique l'Ingnierie et l'Environnement
Module B3-5: Introduction aux lments finis
M. Dysli

modification 1997-06-25

Principe de la rsolution par le mthode des dplacements

Point de calcul

ne figure pas
dans l'quation

Equation d'quilibre
des dplacements
..
. v = R - F+ U
M.v +

R2

des

R1

R
k

k31 k

k11 k 21 k

k13 k12 k 22 k

{}=

k23 k

k32 k

k24 k34

k44 k

k 42 k

k43 k

k33 k

Relation entre le vecteur des forces

aux noeuds F de l'lment et celui des
dplacements v correspondants :

{F} =

eb

an

{}

voire lasto-plastique
D = loi de comportement

k14 k

D = matrice d'lasticit

k41

de

{v}

rd

e signifie "de l'lment"

rg

eu

symtrique

Ke

la

matrice K : voir document No 11a

Matrice de rigidit de l'ensemble de la structure :

arrangement des sous-matrices kij en fonction de la
numrotation globale de tous les noeuds du rseau.
Dans cette matrice les termes kij sont
en fait la somme des termes occupant
la mme position (provenant des
lments voisins).

k 11 k 12 k13 k 14
e

k 21 k 22 k23 k 24

K = k k k k
31 32 33 34
k 41 k 42 k 43 k 44

avec : M =
K=
v=
R=
F=

matrice de rigidit de l'lment

qui tablit la relation entre le vecteur des forces F et le vecteur des
dplacements v aux noeuds 1,2,3,4.
F = Ke . v
sous-matrice relative
chaque noeud, de dimension
k ij = .
dl dl, dl = degrs de libert
du noeud

matrice des masses.

matrice de rigidit.
u
vecteur des dplacements, dnomm u ou dans les prcdents documents.
v
vecteur des forces extrieures.
vecteur des forces aux noeuds du systme; forces quivalentes aux contraintes
effectives.
U = vecteur de la rsultante aux noeuds des pressions interstitielles.

Cycle postgrade
: Gologie Applique l'Ingnierie et l'Environnement
Module B3-5: Introduction aux lments finis
M. Dysli

Document
No 16

FEM, BEM, KEM

FEM = Finite Element Method

=
MEF = Mthode des Elments finis

Un seul lment
avec rsolution
sa frontire qui
sera "tu" pour
simuler l'excavation

Simulation creuse tunnel

BE(M) = Boundary Element

(Method)

Elment eau
avec rsolution
sa frontire

Elments avec rsolution

leur frontire

= (Mthode des) Elments

aux frontires
(Mthode des lments finis
avec lments spciaux o les
dplacements ne sont connus
qu' leurs "frontires")

Effet sisme sur un barrage poids

t = 0,26 sec.

t = 0,51 sec.

t = 0,75 sec.

Tous les lments sont indpendants avec joints

cinmatiques entre eux.

KEM = Kinematic Element

Method
=
Mthode des Elments
cinmatiques
(Mthode des lments finis avec
lments indpendants et joints
cinmatiques)
Mthode dnomme aussi:
DEM = Discret Element Method