Rdmtheo

RDM
Complments thoriques
Bibliographie
Yves DEBARD
Institut Universitaire de Technologie du Mans
Dpartement Gnie Mcanique et Productique
Yves DEBARD
juin 2002
Institut Universitaire de Technologie

Dpartement Gnie Mcanique et Productique
Avenue Olivier Messiaen
72085 Le Mans Cedex 9
Tel 02 43 83 34 64
Fax 02 43 83 31 49
E-mail : [email protected]
http://iut.univ-lemans.fr/ydlogi
http://iut.univ-lemans.fr/gmp/cours/rdmyd
RDM
Complments thoriques
Bibliographie
Sommaire :
Pages
Elasticit
25
Section droite
32
Calculs
43
Maillage dun domaine plan : triangulation de Delaunay
48
Bibliographie
Elasticit
Elasticit
Nous rappelons dans ce texte, les principaux rsultats de la thorie de l'lasticit. Le matriau est
homogne et isotrope. Son comportement est linaire et lastique. Les dplacements sont petits.
A - Contraintes autour dun point

I - Dfinitions
En chaque point M d'un solide, il existe des forces intrieures que l'on met en vidence en effectuant
une coupure du solide, par une surface S, en deux parties A et B.
B
n
S
dS
dF
La partie A, par exemple, est en quilibre sous l'action des forces extrieures qui lui sont directement
appliques et des forces intrieures rparties sur la coupure.
r
Considrons un point M de S. Soit dS un lment infinitsimal de la surface S, entourant M et n le
vecteur unitaire, perpendiculaire en M S et dirig vers l'extrieur de la partie A.
Nous appellerons " facette n en M" cet lment de surface.
Soit dF la force qui s'exerce sur cette facette. On appelle contrainte en M et dans la direction n , la
quantit :
r
r
r
dF
( M, n) = lim
dS0 dS
r
r
r
r
On a : ( M, n) = ( M , n) ( galit de l'action et de la raction ).
RDM - Complments thoriques
Le vecteur contrainte peut tre dcompos en sa composante suivant n et sa projection sur la facette :
r
r
r r
( M, n) = n + .
(M,n)
r
r r
On a donc = n ( M , n) .
r
est la contrainte normale et est le vecteur cisaillement.
est une valeur algbrique positive ( traction ) ou ngative ( compression ).
II - Etat de contrainte en un point

1 - Dfinitions
r r r
Les vecteurs unitaires ( x , y , z ) associs au repre orthonorm { O , x y z } dfinissent en un point
M du solide trois facettes perpendiculaires entre elles.
(M,z)
z
(M,y)
z
O
x
y
x
(M,x)
Les contraintes qui s'exercent sur chacune de ces faces sont dfinies par leurs composantes dans le
repre { O , x y z } :
r
r
r
r
r r
facette x : (M, x ) = xx x + xy y + xz z
r
r
r
r
r
r
facette y : (M, y) = yx x + yy y + yz z
r
r
r
r
r r
facette z : (M, z) = zx x + zy y + zz z
r
Remarque : sur la facette x , la contrainte normale est gale :
xz
r
r r
= x (M, x ) = xx
(M,x)
et le vecteur cisaillement est gal :

r
r
r
= xy y + xz z .
xy
xx
Elasticit
2 - Contrainte sur une facette quelconque : tenseur des contraintes

r
r
Considrons la facette n en M. Soit a, b et c les cosinus directeurs de n .
r
La contrainte sur la facette n est gale :
(M,n)
r
r
r
r
r
r
r
r
r
( M, n ) = a ( M, x ) + b ( M, y) + c ( M, z ) = ( M ) n
n
z
Sous forme matricielle :

r
[(M, n )] = [(M )] [n ]
o [(M)], appel tenseur des contraintes en M, a pour expression :

r
r
( M, x )
r
r
( M , y)
r
r
( M , z)
xx xy xz
[(M)] = yx yy yz
zx zy zz
zz
zy
zx
z
xz
y
xx
xy
yz
yy
yx
Composantes du tenseur des contraintes
3 - Rciprocit des contraintes tangentielles

On a les relations suivantes : xy = yx ,
xz = zx ,
zy = yz
4 - Etat de contrainte en un point

Dans un solide, la distribution des contraintes autour dun point M est entirement dfinie par la
donne du repre { M , x y z } et des six quantits : ( xx yy zz xy xz yz ) . Nous appellerons
cet ensemble tat de contrainte en M.
5 - Contrainte normale dans une direction quelconque

r r
r
= n (M, n ) = [n ] T [(M )] [n ] = a 2 xx + b 2 yy + c 2 zz + 2 a b xy + 2 a c xz + 2 bc yz
III - Faces et contraintes principales
Existe t-il en M
r une facette n telle que le vecteur contrainte soit
colinaire avec n ? Dans ce cas, le vecteur cisaillement est nul sur cette
r
r
facette et le vecteur contrainte ( M, n) satisfait la relation :
n
(M,n)
r
r
r
r
( M, n ) = ( M ) n = n
soit sous forme matricielle : [(M)] [ n ] = [n ] .
r
est alors valeur propre du tenseur des contraintes et n est le vecteur propre associ.
( M ) est une matrice symtrique coefficients rels. Elle a trois valeurs propres relles ( distinctes
ou confondues ). Si les trois valeurs propres sont distinctes, les vecteurs propres associs sont
perpendiculaires entre eux.
Il existe donc en M un repre orthonorm { M , X Y Z } tel que sur les facettes X,Y et Z le vecteur
cisaillement est nul.
Soit X , Y , Z les contraintes normales associes.
Remarque : les trois contraintes principales sont les solutions de lquation : dt ( [(M )] [I] ) = 0
o [I] est la matrice unit.
Dans le repre { M , X Y Z }, le tenseur des contraintes s'crit :
Y
X 0 0
[(M)] = 0 Y 0
0 0 Z
X
X
z
M
Les facettes X , Y , Z sont appeles faces principales.
Les quantits X , Y et Z sont les contraintes

principales.
Dans le repre principal { M , X Y Z }, le vecteur contrainte sur la facette n a pour expression :

X a X
r
r
(M, n ) = Y = b Y
Z c Z
o a, b et c sont les composantes de n .

De la formule prcdente et de la relation a2 + b2 + c2 = 1, on dduit :
2X 2Y 2Z
+
+
= 1.
2X 2Y 2Z
r
r
r
Quand n varie, l'extrmit du vecteur ( M, n) se dplace sur l'ellipsode d'axes ( X , Y , Z ) appel
ellipsode de Lam.
Elasticit
IV - Cercles de Mohr des contraintes

En M, prenons comme repre, les axes rprincipaux X, Y et Z. Considrons la famille de facettes
passant par la direction principale Z. Soit n ( a , b , 0 ), une de ces facettes. Sur cette facette, le vecteur
contrainte est gal :
a X
b avec a2 + b2 = 1.
Y
0
r
r
Le vecteur contrainte ( M, n) est donc situ dans le plan { M , X Y }.
(M,n)
facette Y
Y
(d,0)
-2
X
facette X
r
facette n
Soit t le vecteur unitaire, situ dans le plan { M , X Y } et faisant avec n un angle gal /2 .
r
r r r
r
r
Projetons le vecteur contrainte sur les axes n et t : ( M, n) = n + t
Si nous appelons l'angle que fait direction n avec la direction principale X, et s'crivent :
= X cos 2 + Y sin 2
= X cos sin + Y cos sin
= d + r cos( 2)
soit
= r sin( 2)
avec d = 12 ( X + Y ) et r = 12 ( X Y )
A chaque facette n , nous pouvons donc associer un point dans le repre ( , ) orthonorm. Lorsque
langle varie, ce point dcrit le cercle de rayon r et centre (0,d).
r
Remarque : si la facette n fait un angle avec la facette X, son
point reprsentatif sur le cercle de Mohr fait un angle -2 avec le
point reprsentatif de la facette X.
Les trois directions principales nous permettent de construire trois

cercles.
,
Z
r
On montre que le point reprsentatif ( , = ) d'une facette
quelconque en M reste l'intrieur du domaine limit par les trois
cercles.
Cercles de Mohr en M
V - Etats de contraintes particuliers

1 - Etat de contrainte uniaxial : traction ou compression simple
L'tat de contraintes en un point M est uniaxial si, dans le repre principal { M , X Y Z }, le tenseur
des contraintes se rduit :
0 0
[(M)] = 0 0 0
0 0 0
y
x
Cet tat de contraintes est appel tat de traction simple si est positif et tat de compression simple
si est ngatif.
2 - Etat de cisaillement simple

Soit en M un repre orthonorm { M , x y z }. L'tat de contraintes en M est un tat de cisaillement
simple par rapport aux deux directions x et y, si le tenseur des contraintes se rduit :
0 0
[(M)] = 0 0
0 0 0
Dans le repre principal { M , X Y Z }, le tenseur des contraintes est gal :

0 0
[(M)] = 0 0
0 0 0
3 - Etat de contrainte isotrope
L'tat de contraintes en un point M est isotrope si, quelque soit la facette n , on a T ( M, n ) = n .

Les trois contraintes principales sont alors gales et le tenseur des contraintes en M a pour
expression ( quelque soit le repre ) :
0 0
[(M)] = 0 0
0 0
r
Toute facette n est face principale. En M, les cercles de Mohr des
contraintes se rduisent un point.
Elasticit
4 - Etat de contrainte plan

En un point M, l'tat de contrainte est dit plan s'il existe un repre orthonorm { M , x y z } tel que le
tenseur des contraintes soit de la forme :
xx xy 0
[(M)] = xy yy 0
0
0 0
L'axe z est donc direction principale et la contrainte

principale correspondante est nulle.
VI - Equations d'quilibre
1 - Forces de volume
r
Soit f , de composantes ( fx , fy , fz ), la force par unit de volume applique au point M du solide. Les
quations d'quilibre en M s'crivent :
xx xy xz
2 u
+
+ fx = 2
+
y
z
t
x
yx yy yz
2 v
+ fy = 2
+
+
z
y
t
x
zx
2 w
zy
+ zz + f z = 2
+
z
y
t
x
o u,v et w sont les composantes du vecteur dplacement du point M et la masse volumique du

matriau.
2 - Forces de surface
r
Considrons en M, une facette n appartenant la frontire d'un solide. Soit a, b et c les cosinus
r
r
directeurs de n . Soit p de composantes ( px , py , pz ) la force par unit de surface, qui sexerce sur la
facette.
r
r
r
Les quations d'quilibre en M s'crivent ( M , n) = p soit :
xx
xy
xz
o ( xx
yy
xy
yy
yz
zz
xz a p x
yz b = p y
zz c p z
xy
xz
yz ) est l'tat de contrainte en M.
B - Dformations autour dun point

Sous l'action des forces extrieures, le solide se dforme. Il en rsulte pour tous les points du solide un
dplacement que nous supposerons petit.
I Vecteur dplacement
z
Soient { O , x y z } un repre orthonorm et M(x,y,z) un point

du solide. Au cours de la mise en charge, le point M vient en
M'.
N
O
On appelle vecteur dplacement du point M le vecteur MM ' .
Nous noterons ses composantes :

u( x, y , z)
v( x, y , z)
w ( x, y, z)
dx
Soit N un point du solide voisin de M : MN = dy . Au cours de la mise en charge N vient en N. Le
dz
u ( x + dx, y + dy, z + dz)
dplacement du point N est gal NN ' = v( x + dx , y + dy, z + dz) soit
v( x + dx , y + dy, z + dz)
u
u
u
u
dx +
dy +
dz
y
z
u ( x , y, z) x
u x
v
v
v
v
NN ' = v( x , y, z) +
dx +
dy +
dz = v +
x
y
z

v( x , y, z)
w
w
w w w
dx +
dy +
dz
y
z
x
x
En dcomposant la matrice [D] en sa partie symtrique [E ] =

antisymtrique [] =
z dx
v
dy = MM ' + [D] MN
z

w dz
u
y
v
y
w
y
1
( [D] + [D] T ) et sa partie
2
1
( [D] [D] T ) , il vient:
2
NN ' = MM ' + [] MN + [E] MN
xx xy xz
avec : [E ] = yx yy yz
zx zy zz
et
0 z y
[ ] = z
0 x
y x
0
o lon a pos :
u
xx =
x
w
1 u v
1 u w
1 w v
xz =
xy =
+
+
+
yz =
z
2 y x
2 z x
2 y z
1 w v
1 u w
1 v u
x =

y =
z =

2 y z
2 z x
2 x y
yy =
v
y
zz =
Elasticit
r
Le dplacement du point N scrit finalement : NN ' = MM '+ MN + [E ] MN
Si [] et [E] sont nuls, le dplacement de N se rduit NN' = MM' et tous les points situs au
voisinage de M subissent la mme translation.
r
Si MM ' et [E] sont nuls le dplacement de N se rduit : NN' = MN . Si de
r
plus = est petit (les drives du vecteur dplacement sont petites), tous les
points situs au voisinage de M subissent une rotation dintensit autour de laxe

r 1
= rot (MM ') dont lorigine est situe en M.
2
M
La matrice [E] qui reprsente les dformations du solide en M est appele tenseur
des dformations.
Remarque : des relations ci-dessus, on dduit :
M ' N' = MN + [D] MN = MN + [] MN + [E] MN = ( [I] + [D] ) MN = ( [I] + [] + [E ] ) MN
MN = M ' N' MN = [D] MN = [] MN + [E ] MN = ([] + [E]) MN
II Etat de dformation en un point

1 - Allongement unitaire ou dilatation
y
Considrons deux points M et N du solide voisins l'un de

l'autre. Au cours de la mise en charge, le point M vient en M' et
r
le point N en N'. Soit n le vecteur unitaire li la direction
N
M
M
MN .
r
On appelle allongement unitaire en M dans la direction n , la
quantit :
N
y
r
M ' N'MN
(MN)
( M, n ) = lim
= lim
NM
N M MN
MN
a
r
r
Soient a , b , c les composantes de n et l la longueur de la fibre MN do : MN = l n = l b = l [n ] .
c
r
2
2
2
2
Evaluons la quantit (MN ) : (MN ) = l = 2 l l = 2 l (M, n )
Dautre part : (MN 2 ) = (MN.MN) = MN.MN + MN.MN = 2 l 2 [n ] T [E] [n ]

On obtient finalement :
r
( M, n ) = [n ]T [E ][n ] = a 2 xx + b 2 yy + c 2 zz + 2 a b xy + 2 a c xz + 2 bc yz
r r
r
r
Remarque : si n = x , on obtient : (M, x ) = xx . De mme : (M, y) = yy et
r
(M, z) = zz . Les
quantits xx , yy et zz reprsentent donc respectivement lallongement unitaire en M dans les

r r r
directions x , y et z .
10
2 - Glissement ou dformation angulaire

N2
Soient deux points N1 et N2 voisins de M et tels que les

r
r
directions MN1 et MN 2 soient orthogonales. Soient n 1 et n 2
les vecteurs unitaires associs ces deux directions. Au cours
de la mise en charge les points M, N1 et N2 viennent
respectivement en M', N1' et N2'. Soit l'angle que font entre
n2
r
r
On appelle glissement en M dans les directions n1 et n 2 , la
quantit :
r r
(M, n1 , n 2 ) = lim ( )
N1 N 2
N2
eux les deux vecteurs MN '1 et MN' 2 .

O
x
N1
n1
N1
y
N2 N
a 1
a 2
r
r
Posons : MN1 = l1 MN 2 = l 2 n 1 = b1 n 2 = b 2 .
c1
c 2
a 1
r
r
On a donc : MN1 = l1 n 1 = l1 b 1 , MN 2 = l 2 n 2 = l 2
c1
a 2
b
2
c 2
Soit langle que font entre eux les vecteurs MN1 et MN 2 .

Evaluons la quantit (MN 1 . MN 2 ) :
r r
(MN 1 . MN 2 ) = (l1l 2 cos ) = (l1 l 2 ) cos l1 l 2 sin = l1l 2 ( M, n 1 , n 2 ) car = .

2
Dautre part : (MN 1 . MN 2 ) = MN 1 . MN 2 + MN 1 . MN 2 = 2 l1l 2 [n 1 ]T [E ][n 2 ] .

r r
On en dduit ( M, n 1 , n 2 ) = 2 [ n 1 ] T [E ][n 2 ] soit :
r r
( M, n 1 , n 2 ) = 2 [ xx a 1a 2 + yy b1 b 2 + zz c1c 2 + xy (a 1 b 2 + a 2 b1 ) + yz (a 2 b 3 + a 3 b 2 ) + xz (a 3 b1 + a 1 b 3 )]
r
r
r r
r
r
Remarque : si n 1 = x et n 2 = y , lexpression ci-dessus se rduit : ( M, x , y) = xy = 2 xy . De
r r
r r
mme : ( M, x , z) = xz = 2 xz et ( M, y, z) = yz = 2 yz . Les quantits xy , xz et yz reprsentent
r r r r
r r
donc respectivement le glissement en M dans les directions (x , y ) , (x , z ) et (y, z ) .
3 - Etat de dformation en un point

La distribution des dformations autour d'un point M est entirement dfinie par la donne des six
quantits : ( xx yy zz xy xz yz ) .
Nous appellerons cet ensemble tat de dformation en M.
Elasticit
11
III - Faces principales

En M, dans le repre principal { M , X Y Z }, on a : XY = 0 , XZ = 0 , ZY = 0 et le tenseur des
dformations se rduit :
X 0 0
[E ] = 0 Y 0
0 0 Z
Les quantits X , Y et Z sont appeles dformations principales.
IV - Cercles de Mohr des dformations

En M, prenons comme repre, le repre principal { M , X Y Z }. Considrons la famille de facettes
r
r
passant par la direction principale Z. Soit n ( a , b , 0 ), une facette appartenant cette famille et t le
r
vecteur unitaire, situ dans le plan { M , X Y } et faisant avec n un angle gal /2.
/2
t
/2
(d,0)
r
-2
facette n
r
A chaque facette n , nous pouvons associer deux quantits et dfinies par :
r r
r
= ( M , n) et = ( M , n, t ) .
r
Si nous appelons l'angle que fait direction n avec la direction X, on a :
= X cos 2 + Y sin 2
= 2 X cos sin + 2 Y cos sin
soit
d = 12 ( X + Y )
= d + rcos(2)
avec
2 = rsin(2)
1
r = 2 ( X Y )
r
A chaque facette n , nous pouvons associer un point ( , /2 ) dans un repre orthonorm. Lorsque
varie, ce point dcrit le cercle de rayon r et centre ( 0 , d ).
Les trois directions principales nous permettent de construire trois cercles.
12
C - Loi de comportement
I - Bases exprimentales
Considrons une barre de section droite constante et rectangulaire. L'axe X est la ligne moyenne de la
barre. Soit O un point quelconque du solide. La barre, constitue d'un matriau homogne et isotrope
est soumise suivant l'axe X un effort de traction uniformment rparti chaque extrmit.
Le repre orthonorm { O , X , Y , Z } est le repre principal en O.
Faisons crotre la force de traction partir de 0. L'exprience montre que :
1 - Loi de Hooke
Y
-X
On a : X = E X o E est le module d'lasticit longitudinal ( module d' Young ) du matriau.

La relation prcdente est vrifie si la contrainte XX ne dpasse pas une certaine valeur E appele
limite lastique en traction.
2 - Loi de Poisson
L'allongement de la barre suivant la direction x, s'accompagne d'une contraction suivant les directions
y et z :
Y = - X , Z = - X o est le coefficient de Poisson du matriau ( 0 < < 0.5 ).
II - Notations
L'tat de contrainte en un point sera reprsent par un vecteur six composantes [] dfini par :
[]T = [ xx yy zz xy xz yz ]
L'tat de dformation en un point sera reprsent par un vecteur six composantes [] dfini par :
[] T = [ xx yy zz xy xz yz ]
Elasticit
13
III - Loi de comportement

Le vecteur contrainte [] et le vecteur dformation [] sont lis par la relation :
[] = [D] [] + [ TH ]
La matrice [D] des coefficients lastiques est gale :
+ 2

+ 2
+ 2
[ D] =
0
0
0
0
0
0
0
0
0
0 0 0
0 0 0
E
= (1 + )(1 2 )
0 0 0
avec
E
0 0
=
2(1 + )
0 0
0 0
o E et sont respectivement le module d'Young et le coefficient de Poisson du matriau.

et sont les coefficients de Lam.
Remarque : Le module de glissement G est gal .
1
1

E(T To ) 1
[TH] reprsente les contraintes d'origine thermique et est gal : [ TH ] =

1 2
0
0

0
o est le coefficient de dilatation thermique et To la temprature de rfrence (temprature de
montage).
IV - Dformation volumique
Soit V, un lment de volume infiniment petit entourant le point M. Au cours de la dformation V
augmente de dV. On appelle dformation volumique la quantit :
V =
dV
(1 2 )
= xx + yy + zz =
( xx + yy + zz )
V
E
D - Energie de dformation
Soit dV = dx dy dz un lment de volume, infiniment petit, entourant le point M. Au cours de la mise
en charge l'nergie potentielle de dformation accumule dans cet lment de volume est gale :
dW = 12 []T [] dV = 12 [] T [D] [] dV
14
E - Problmes particuliers dlasticit

I - Contraintes planes
Dfinition : Un solide est en tat de contraintes planes par
rapport au plan { O , x y }, s'il existe un repre { O , x y z }, li
au solide, tel qu'en tout point M du solide, le tenseur des
contraintes soit de la forme :
xx xy 0
[(M)] = yx yy 0
0
0 0
o xx , yy et xy sont indpendants de z.
L'axe z est donc, pour tous les points du solide, direction principale et la contrainte principale
associe est nulle. On en dduit:
xz = yz = 0
zz =
( xx + yy )
E
L'tat de dformation est donc entirement dfini par les trois quantits : xx , yy et xy.
La loi de comportement scrit :
xx
E
yy =
2
xy 1
1
0
0 0 (1 ) / 2
xx
1
E(T To )
yy
1
1
xy
0

Les dformations et les contraintes ne dpendent que des dplacements u( x, y ) et v( x, y ) parallles

aux axes x et y.
Domaine d'application : L'approximation contraintes planes convient aux plaques minces sollicites
dans leur plan. Le plan { O, x y } est alors le plan moyen de la plaque.
II - Dformations planes
Dfinition : Un solide est en tat de dformations planes par rapport au plan { O , x y }, s'il existe un
repre { O x y z }, li au solide, tel qu'en tout point du solide, le champ de dplacement soit de la
forme :
u=u(x,y) , v=v(x,y) , w=0
On en dduit : zz = 0
xz = 0
yz = 0
zz = ( xx + yy )
En tout point du solide, la direction z est donc direction principale. Les dformations et les contraintes
sont indpendantes de z.
Elasticit
15
Le tenseur des contraintes est alors de la forme :

xx xy 0
[(M)] = xy yy 0
0
0 zz

xx + 2
0 xx
1
E(T To )
1
+ 2 0 yy
yy =
1 2
xy 0
0
xy
0
Domaine d'application : L'tat de dformations planes se prsente lorsqu'on a affaire un cylindre

d'axe Oz trs long satisfaisant aux conditions suivantes :
- les bases du cylindre sont fixes.
- les forces appliques au solide sont normales l'axe Oz et indpendantes de z.
III - Problme axisymtrique

Le solide considr est de rvolution. Il en va de mme du chargement et des conditions aux limites.
Soit z l'axe de rvolution. Un point du solide est repr par ses coordonnes cylindriques ( r , , z ).
La solution est axisymtrique. Chaque point du solide se dplace dans son plan mridien ( r , z ). De
plus le champ de dplacement est indpendant de la coordonne .
z
On a donc :
- dplacement radial : u = u ( r , z )
- dplacement orthoradial : v = 0
- dplacement axial : w = w ( r , z )
r
y
On en dduit les dformations :

rr =
u
r
u
r
zz =
w
z
La direction est direction principale.

Le tenseur des contraintes est gal :
rr 0 rz
[(M)] = 0 0
rz 0 zz
rz =
u w
+
z r
r = z = 0
16
rr + 2

+ 2
=
zz
+ 2

0
0
rz 0
0
0
0
rr
1

E(T To ) 1
zz
1 2 1

0
rz
IV - Flexion des plaques

Dfinitions :
Une plaque est un corps solide limit par deux faces planes
parallles et par une surface cylindrique perpendiculaire cellesci.
L'paisseur e de la plaque est la distance entre les deux faces.

Le plan quidistant des deux faces est le plan mdiant ou
surface moyenne.
y
x
Soit { O , x y z } un repre orthonorm tel que le plan { O , x y } soit le plan moyen. L'axe z est
donc normal au plan moyen.
Le plan situ z = e/2 est la peau suprieure de la plaque.
Le plan situ z = -e/2 est la peau infrieure de la plaque.
Une fibre normale est l'ensemble des points du solide situs sur une normale au plan mdiant. Elle
est caractrise par la donne de ses coordonnes ( x , y ).
Une plaque est dite mince si son paisseur est petite par rapport aux autres dimensions.
Hypothses :
La plaque est sollicite par des forces perpendiculaires au plan moyen et des couples de composantes
(mx , my , 0).
La contrainte normale zz est ngligeable par rapport aux autres composantes du tenseur des
contraintes.
Les phnomnes de membrane et de flexion sont dcoupls. Compte-tenu des conditions de
chargement, les phnomnes de membrane sont nuls.
Au cours de la dformation, une fibre normale reste droite mais ne reste pas perpendiculaire au plan
moyen.
Elasticit
17
Champ de dplacement :
Considrons une fibre normale ( x , y ). Soit Mo le point d'intersection de cette fibre avec le plan
mdiant et M(x,y,z) un point quelconque de cette fibre.
Le vecteur dplacement ( u , v , w ) en M est gal :
u ( x , y, z ) = z y ( x , y )
v ( x , y, z ) = z ( x , y )
x
w ( x , y, z) = w ( x , y)
o
- w est le dplacement du point Mo suivant z. w est la flche en Mo.
- x est la rotation de la fibre normale suivant x.
- y est la rotation de la fibre normale suivant y.
y
w
z
y
Mo
Le champ de dplacement dans le solide est donc dfini par la connaissance de ( w , x , y ) en tout
point ( x , y ) du plan moyen.
De l'expression du champ de dplacement, on dduit les dformations :
y
, yy = z x
xx = z
x
y
y x
xy = z
y
y
w
w
xz = y + x , yz = x + y

xx

yy
xy = E
1 2
xz
yz
0
0
0 (1 ) / 2
0
0
0
0
(1 ) / 2
0
0
0
0
(1 ) / 2
0
0
0
0
0
0
xx

yy
xy

xz
yz

Remarque : Avec les hypothses prcdentes, le cisaillement transverse est pris en compte. En effet
les contraintes xz et yz sont diffrentes de zro.
18
Si on nglige le cisaillement transverse ( hypothse de Kirchhoff ), on a les relations suivantes :

x =
w
,
y
y =
w
x
Les fibres normales restent alors perpendiculaires au plan moyen, au cours de la dformation.
Rsultantes et moments :
Considrons un lment de plaque infiniment petit, limit par un cylindre perpendiculaire au plan
moyen, de section droite rectangulaire et dont les faces sont parallles x ou y.
La face AB est soumise :
z
- un moment de flexion ( a Mxx ).
Qyz
Qxz
- un moment de torsion ( a Mxy ).
Mxy
- une force de cisaillement : ( a Qxz ).
La face CB est soumise :
Myx
- un moment de flexion ( b Myy ).
Mxx
- un moment de torsion ( b Myx ).
- une force de cisaillement : ( b Qyz ).

avec :
e/2
M xx =
z xx dz
e/2
M yy =
e / 2
z yx dz
e / 2
z yy dz
e/2
M xy =
e / 2
e/2
M yx =
Myy
e/2
Q xz =
xz dz
e / 2
xy dz
e / 2
e/2
Q yz =
yz dz
e / 2
On en dduit :
xx = z Mxx / I , yy = z Myy / I , xy = z Mxy / I avec I = e3 / 12.
Elasticit
19
F - Critres de limite lastique

Soit X , Y et Z les trois contraintes principales en un point M d'un solide. Nous supposerons que la
limite lastique en traction simple est gale la limite lastique en compression simple. Soit E cette
limite lastique .
Comment vrifier, dans un tat de contrainte complexe, que la limite lastique n'est pas dpasse ? On
admet que la limite lastique est atteinte lorsqu'une certaine fonction f des contraintes principales est
gale limite lastique du matriau en traction simple :
f (X , Y , Z ) = E.
Le domaine lastique en un point du solide est donc dfini par la relation :
f ( X , Y , Z ) < E.
Nous examinons dans ce chapitre plusieurs critres de limite lastique.
I - Critre de Rankine ou de la contrainte normale maximale

1 - Enonc :
La fonction f est gale : f (X , Y , Z ) = sup ( X , Y Z ).
2 - Validit :
Y
Dans un tat de traction simple, le critre est satisfaisant. Dans un tat

de cisaillement pur, le critre impose E = E o E est la limite
lastique au cisaillement pur.
3 - Etat plan de contraintes ( Z = 0 ) :
Le domaine lastique a la forme suivante dans le plan {X , Y } :
-E
-E
II - Critre de Tresca ou du cisaillement maximal

1 - Enonc :
On vrifie que le cisaillement maximal est plus petit que la moiti de la limite lastique en traction
simple :
f (X , Y , Z ) = 2 max < E.
La quantit f est appele contrainte quivalente de Tresca .
20
2 - Validit :
Y
Ce critre est satisfaisant dans un tat de traction simple, mais il

impose E = 2E dans un tat de cisaillement pur.
3 - Etat plan de contraintes ( Z = 0 ):
E
E
-E
-E
III - Critre de Von Mises

1 - Enonc :
On vrifie ( X - Y )2 + ( X - Z )2 + ( Y - Z )2 < 2 E2 .
La quantit f =
1
[( X Y ) 2 + ( X Z ) 2 + ( Y Z ) 2 ] est appele contrainte quivalente de
2
Von Mises.
2 - Validit :
Ce critre est satisfaisant dans un tat de traction simple, et il impose E = 0.57 E dans un tat de
cisaillement pur.
Y
3 - Etat plan de contraintes ( Z = 0 ):

E
E
-E
-E
Elasticit
21
G - Dpouillement des rosettes

Une rosette est un ensemble de trois jauges colles en un point M
dun solide et faisant entre elles un angle gal . Langle est gal
45 , 60 ou 120.
y
3
2
r
Soit z le vecteur unitaire normal la surface et dirig vers
l'extrieur du solide et { M, x y zr } le repre orthonorm tel que
laxe x soit colinaire avec la jauge 2 .
r
La direction z est direction principale et la contrainte principale correspondante est nulle (en
labsence de pression extrieure) : en M l'tat de contrainte est plan. Le tenseur des contraintes et le
tenseur des dformations ont pour expression:
xx xy 0
xx xy 0
[(M)] = yx yy 0 , [E (M )] = yx yy 0
0 0
0
0 0 zz
r r
r
On mesure les allongements unitaires dans les trois directions 1 , 2 et 3 .
r
Posons : 1 = (M,1) ,
r
2 = ( M , 2) ,
r
3 = (M,3) .
Les trois mesures effectues vont nous

contraintes en M.
permettre de calculer les composantes du tenseur des
Posons : c = cos () et s = sin ().

Des relations :
1 = xx c 2 + yy s 2 xy sc
2 = xx
= c 2 + s 2 + sc
xx
yy
xy
3
on dduit les composantes du tenseur des dformations :

xx = 1
1 + 3 2 2 c 2
yy =
2s 2
1
xy = 3
2sc
zz = 1 ( xx + yy )
y
3
-
1
22
puis les composantes du tenseur des contraintes :
E
( xx + yy )
xx =
1 2
= E ( + )
xx
yy 1 2 yy
E
xy
xy =
+ )
2
(
1
o E et sont les caractristiques lastiques du matriau.

On en dduit les contraintes principales :
X = d + r
= d r avec d = xx + yy et r 2 = xx yy
Y
2
2
Z = 0
+ 2xy
les dformations principales :

1
X = E ( X Y )
1
Y = ( Y X )
E
= ( + ) =
X
Y
zz
Z
E
et la position angulaire de la direction principale X par rapport laxe x :

xx yy
r cos(2) =
2
r sin(2) = xy
Elasticit
23
Section droite
Considrons une section droite dune poutre.
C
y
Soit { Y et Z } le repre central principal. L'axe X est l'axe neutre de la poutre. Le point ( Y = 0 ,
Z = 0 ) est le centre lastique de la section. On a les relations suivantes :
E (Y, Z) Y dA = 0
E (Y, Z) Z dA = 0
E (Y, Z) Y Z dA = 0
o E(Y,Z) est le module d'Young du matriau.
Les caractristiques homognises de la section sont :

- la rigidit de membrane :
< EA >=
E (Y, Z) dA
- les rigidits principales de flexion :
< EI Y >=
E (Y, Z) Z 2 dA
< EI Z >=
E (Y, Z) Y 2 dA
- la rigidit de cisaillement transversal :

< GA >=
G (Y, Z) dA
- la rigidit de torsion :
< GJ >=
G (Y, Z) Y
Z
+ Y 2 + Z 2 dA
A
Y
Z
Dans ces formule G(Y,Z)

gauchissement de torsion.
est le module de cisaillement du matriau et la fonction de
24
La force lastique sur la section droite a pour expression :
[N X
TY TZ M X M Y M Z ]
L'nergie de dformation linique est gale :

M 2Y
M 2Z
TY2
TZ2
M 2X
dW 1 N 2X
=
+
+
+
+
+
dX 2 < EA > < EI Y > < EI Z > < Gk Y A > < Gk Z A > < GJ >
o kY et kZ sont les coefficients daire cisaille.

La contrainte normale est gale :
NX
MY
MZ
Z
Y
XX = E (Y, Z)
+
< EI Z >
< EA > < EI Y >
Les contraintes de cisaillement sont gales :

M X
1
g
h
XY = G (Y, Z)
+ TZ
Z
TY
+
Y
Y < GJ > Y
< GA >
M X
1
g
h
+ Y
XZ = G (Y, Z)
+ TZ
TY
+
Z
Z < GJ > Z
< GA >
o g et h sont les fonctions de gauchissement d'effort tranchant sur Y et Z.

Le centre de cisaillement/torsion est le point de la section qui reste fixe lorsque la force lastique
sur la section se rduit un moment de torsion.
Rfrences :
[B5] [B6] [G4] [L4] [L6]
[SEC 1] ... [SEC...]
Section droite dune poutre
25

A - Caractristiques
1 - Dfinitions
Soient une section droite (A) et { y z } un repre orthonorm du plan de
la section.
A = dy dz =
A
dA
On appelle aire de la section la quantit :
O
z
Les moments statiques de la section par rapport aux axes y et z sont

dfinis par :
M y = z dA ; M z = y dA
A
On appelle moments quadratiques de la section par rapport aux axes y et z, les quantits :
I y = z 2 dA ; I z = y 2 dA
A
On appelle moment produit de la section par rapport aux axes y et z, la quantit :

I yz =
y z dA
On appelle centre de gravit de la section le point O dont les coordonnes sont dfinies par les
relations :
yO =
Mz
A
zO =
My
A
Soit une droite (D), passant par lorigine et faisant un angle par rapport l'axe Oy. Le moment
quadratique de la section par rapport (D) est gal :
I() = I y cos 2 () + I z sin 2 () 2 I yz cos() sin()
Lorsque varie, I() passe par un minimum et un maximum. Les deux directions correspondantes
sont perpendiculaires entre elles. Le moment produit est nul par rapport ces deux directions. Ces
deux directions que nous noterons Y et Z sont appeles axes centraux principaux.
26
RDM Complments thoriques
On a les relations suivantes :
Y dA =
Z dA =
IY =
Les quantits
Y Z dA = 0
et
dA
IZ =
dA
sont les moments
quadratiques principaux.
Le logiciel value les quantits suivantes :

- moment dinertie polaire ( moment de la surface par rapport au centre de gravit ) :
IO =
(Y
+ Z 2 dA = I Y + I Z
- constantes de stabilit : Y = Y (Y 2 + Z 2 ) dA
Z = Z (Y 2 + Z 2 ) dA
- modules plastiques : Wpl Y =
ZZ
dA
Wpl Z =
Y Yp dA
o les droites Y = Yp et Z = Zp sparent la surface de la section droite en deux surfaces gales.

- rayons de girations par les relations : i Y =
IY
iZ =
A
IZ
iO =
A
IO
A
La constante de torsion de Saint Venant est gale :

J=
Y Z Z Y + Y
+ Z 2 dA
o est la fonction de gauchissement de torsion. Si la section est circulaire ( pleine ou creuse ), est
nul et J se rduit IO.
Le centre de cisaillement/torsion est le point de la section qui reste fixe lorsque la force lastique
sur la section se rduit un moment de torsion. Il est dfini par :
YC =
1
IY
Z dA
et Z C =
1
IZ
Y dA
Le moment dinertie de rotation ( moment de la surface par rapport au centre de cisaillement ) est
gal :
Ir =
[(Y Y
C)
+ ( Z Z C ) 2 ] dA = (YC2 + Z C2 ) A + I O
27
La constante de gauchissement est dfinie par :
I = 2 dA
A
La constante de stabilit est dfinie par :
= (Y 2 + Z 2 ) dA
A
L'nergie de dformation linique due leffort tranchant est gale :

TY2
TZ2
+
2 Gk Y A 2 Gk Z A
o kY et kZ sont les coefficients daire cisaille. AY = kY A est laire cisaille suivant Y et AZ = kZ A

est laire cisaille suivant Z. Ces coefficients sont dfinis par :
kY =
1
IZ
gY dA ; k
1
IY
h Z dA
o g et h les fonctions de gauchissement associes aux efforts tranchants TY et TZ.
2 - Calculs
Calcul de A , My , Mz , Iy , Iz , Iyz :
Ces quantits sont calcules par intgration sur le contour extrieur

de la section par utilisation de la formule de Green qui permet la
transformation dune intgrale double en intgrale simple de contour :
M
N
dA = N dy + M dz
z
Par exemple, le calcul de Iy seffectue en posant N = 0 et M = yz2 do I y = yz 2 dz
De mme : A = y dz
M y = yz dz
M z = yz dy
I z = zy 2 dy
I yz =
1
y 2 z dz
2
28
Contribution dun segment :

Soit llment de contour dorigine 1 (y1,z1) et dextrmit 2 (y2,z2) :
y = N 1 y1 + N 2 y 2
z=N z + N z
1 1
2 2
dy = ( N 1' y1
'
dz = ( N 1 z 1
+
+
N '2 y 2 ) d
N '2 z 2 ) d
1
1+
N2 =
N1 = 2
2
avec
1
1
N 1' =
N '2 =
2
2
2
=1
1
z
y
=-1
La contribution du segment 1-2 chacune des quantits recherches est de la forme :

+1
f ( ) d
et est calcule numriquement par quadrature de Gauss.
Contribution dun cercle :

La contribution du cercle de centre C (yC,zC), de rayon R et parcouru dans le sens
trigonomtrique aux quantits recherches est gale :
A = R 2
M y = zC A M z = yC A
R2
R
+ z C2 ) A I z = (
+ y C2 ) A I yz = y C z C A
I y = (
4
4
C
z
y
Contribution dun arc de cercle :

La contribution de larc de cercle dorigine 1 (y1,z1), dextrmit 2
(y2,z2), de centre C (yC,zC) et parcouru dans le sens
trigonomtrique laire ( Iy ... ) est gale laire ( Iy ... ) du
segment circulaire (S) moins la contribution laire ( Iy ... ) du
segment 2-1.
Soient R le rayon de larc, son angle au centre et L la longueur
de sa corde.
v
b
a
1
S
2
z
C
Soient G le centre de gravit de (S), Ga et Gb les axes centraux

principaux de (S). Introduisons les axes Gu et Gv parallles
respectivement aux axes y et z.
29
On a ( [C12] [P2] [Y1] ) :

A S = 0.5R 2 ( sin )
I Sa =
d = CG =
R4
( sin cos ) d 2 A S
8
Posons : c = cos =
y1 y 2
L
s = sin =
L3
12A S
I Sb =
R4
(3 4sin + sin cos )
24
z1 z 2
L
y G = y C s d z G = z C + cd
M Sy = z G A S
S
S
M z = y G A
I Su = I Sa c 2 + I Sb s 2
S
S 2
S 2
d'o
I v = Ia s + Ib c
S
S
S
I = I cs + I cs
a
b
uv
I Sy = I Su + z G2 A S
S
S
2
S
Iz = Iv + yG A
I S = I S + y z A S
uv
G G
yz
Remarque : si langle de larc est trs petit, larc est remplac par un lment de contour quadratique.
Contribution dun lment de contour quadratique :

Soit llment de contour quadratique dorigine 1 (y1,z1), de milieu 2 (y2,z2) et dextrmit 3 (y3,z3) :
y = N 1 y1 + N 2 y 2 + N 3 y 3 dy = ( N 1' y1 + N '2 y 2 + N 3' y 3 ) d
z=N z +N z +N z
'
'
'
1 1
2 2
3 3
dz = ( N 1 z1 + N 2 z 2 + N 3 z 3 ) d
avec
(1 )
(1 + )
N2 = 1 2 N3 =
N1 =
2
2
1 + 2
1 + 2
'
'
'
N1 =
N 2 = 2 N 3 =
2
2
=1
3
=0
2
z
y
1
=-1
La contribution de llment 1-2-3 chacune des quantits recherches est calcule numriquement
par quadrature de Gauss ( voir contribution dun segment ).
Calcul de Y , Z :
Ces quantits sont calcules par intgration sur le contour extrieur :
Y = ZY 3 dY + Y 2 Z 2 dZ
2 2
3
Y = Y Z dY + Y Z dZ
Le contour est discrtis en lments quadratiques. La contribution dun lment chacune des
quantits recherches est calcule numriquement par quadrature de Gauss.
30
Calcul des fonctions de gauchissement :

La fonction de gauchissement de torsion (Y,Z) est la solution du problme :
2 2
+
= 0 sur le domaine de la section.
Y 2 Z 2
r r
nY +
n Z = Zn Y Y n Z sur le contour extrieur ( . n = 0 ).
Y
Z
dA = 0
n
A
o nY et nZ sont les cosinus directeurs de la normale extrieure au contour .

La fonction de gauchissement deffort tranchant suivant Y : g (Y,Z) est la solution du problme :
2g 2g
A
+ 2 = Y sur le domaine de la section.
2
IZ
Y
Z
r r
g
g
nY +
n Z = 0 sur le contour extrieur ( . n = 0 ).
Y
Z
g dA = 0
De mme, la fonction de gauchissement deffort tranchant suivant Z : h (Y,Z) est la solution du

problme :
2h 2h
A
+ 2 =
Z sur le domaine de la section.
2
IY
Y
Z
r r
h
h
nY +
n Z = 0 sur le contour extrieur ( . n = 0 ).
Y
Z
h dA = 0
Ces problmes sont rsolus par la mthode des lments finis aprs triangulation de la section par la
mthode de Delaunay. Les lments utiliss reposent sur une formulation dplacement.
Calcul de WplY , WplZ :

Ces quantits sont values en utilisant le maillage de la section.
Bibliographie :
[G4] [B5] [B6] [H2] [L6] [L4] [P1] [SEC 1] ... [SEC...]
31
B - Contraintes
Soient Y et Z les axes centraux principaux de la section. L'axe X est l'axe neutre de la poutre. Le point
O de coordonnes ( 0, 0 ) est le centre de gravit de la section.
Soit [ N TY TZ Mt MfY MfZ ] la force intrieure qui sexerce sur la section.
Z
Au point M(Y,Z), le tenseur des contraintes a pour expression :

XZ
XX XY XZ
[(M)] = XY 0
0
XZ 0
0
M
XY
o la contrainte normale est gale :

XX =
XX
Mf Y
Mf Z
N
+Z
Y
A
IY
IZ
et les contraintes de cisaillement sont gales :

Mt
T
Z + Y
J Y
A
Mt
T
=
+ Y + Y
J Z
A
XY =
XZ
g TZ h
+
Y A Y
g TZ h
+
Z A Z
o , g et h sont respectivement les fonctions de gauchissement associes Mt, TY et TZ (voir

caractristiques dune section droite).
On en dduit :
- les contraintes principales :

1
2
XX 1
2XX + 4 2 avec 2 = 2XY + 2XZ

2
2
3 = 0
- la contrainte quivalente de Von Mises :
VM = 2XX + 3 2
- la contrainte quivalente de Tresca :

T = 2XX + 4 2
Bibliographie :
[G4] [B5] [B6] [H2] [L6] [L4] [P1] [SEC 1] ... [SEC...]
Cercles de Mohr en M
32
Calculs
A - Prsentation
Les principales tapes du calcul sont :
Vrification des donnes
Construction des tables
maillage de lossature ( analyse dynamique, flambement )

construction de la table des noeuds et de la table des lments
partition des degrs de libert
calcul de la place mmoire ncessaire au calcul
segmentation des matrices globales
Calcul des matrices lmentaires et assemblage des matrices globales
matrice de rigidit
matrice masse pour une analyse dynamique
matrice de rigidit gomtrique pour la recherche des charges critiques
vecteur(s) force(s) pour une analyse statique
Factorisation de la matrice de rigidit globale : [K ] = [L] [D] [ L]T .

Analyse statique : pour chaque problme statique,
rsolution du systme dquations linaires [K] [UI] = [FI] o I est le numro du problme :
[1..NCAS, NCAS+1...NCAS+NCOMB]. NCAS et NCOMB sont respectivement le nombre de
cas simples et le nombre de combinaisons. Sil y a des dplacements imposs non nuls, le
logiciel gnre un cas de charges 0 qui se rduit aux dplacements imposs.
dition dun fichier neutre
Analyse dynamique
rsolution du problme aux valeurs propres : [K ] [ U ] = 2 [ M] [ U ] .
Flambement linaire
rsolution du problme aux valeurs propres : [K ] [ U ] = [M ] [ U ] .
A la fin du calcul, le logiciel gnre un fichier dont lextension est .RES. Ce fichier contient des
informations sur le droulement du calcul.
Calculs
33
B - Partition des degrs de libert

Les degrs de libert de la structure (d.d.l.) sont numrots de faon obtenir la partition de ceux-ci
en trois sous-ensembles :
- dplacements inconnus : (L)
- dplacements connus non nuls ( dplacements imposs ou prescrits) : (P)
- dplacements nuls : (S)
La partition des d.d.l. induit une partition des matrices globales : vecteur force [F] , vecteur
dplacement [U] , matrice de rigidit [K] ... :
FL
U L
K LL K LP K LS
[F] = FP [ U ] = U P [K ] = K LP K PP K PS
FS
U S
K SL K SP K SS
...
Les lignes et les colonnes associes aux dplacements nuls ne sont jamais assembles.
La partition des d.d.l . est effectue avant la segmentation en blocs des matrices globales ( voir
mmorisation des matrices globales ).
Analyse statique :
Les quations dquilibre [F] = [K][U] scrivent :
FL K LL K LP U L
F = K

P PL K PP U P
Les dplacements inconnus sont solution du systme dquations :

[K LL ] [ U L ] = [FL ] [K LP ] [ U P ] = [R L ]
Analyse dynamique :
Le problme aux valeurs propres [K ] [ U ] = 2 [M] [ U ] se rduit :
[K LL ] [ U L ] = 2 [M LL ] [ U L ]
34
C - Mmorisation des matrices globales

Un calcul de structure par la mthode des lments finis ncessite souvent d'importantes ressources
informatiques (temps de calcul et place mmoire).
Ces ressources peuvent tre diminues par :
- le choix de mthodes de rsolution qui utilisent les proprits des matrices globales (symtrie,
matrice bande...).
- la renumrotation des noeuds.
Dans RDM, les matrices globales (matrice de rigidit, matrice masse...) sont stockes selon la
mthode profil ou ligne de ciel . Elles sont de plus dcomposes en blocs qui sont mmoriss sur
le disque.
Considrons la structure reprsente sur la figure (10 noeuds, 10 lments) et supposons, pour
simplifier, que chaque noeud possde un seul degr de libert :
1
9
10
L'assemblage conduit la matrice de rigidit globale :

X
.
.
.
[K ] = .
X
.
X
.
. X .
X X X .
X X . .
.
.
.
.
X . X X X .
. . X X . .
.
.
. X . X X
. . . X X
.
.
. . . X .
. X . . .
X .
. X .
. . X
. . .
. X .
. . .
.
X . .
. . .
X . .
. X .
. . X
o chaque X reprsente un terme non nul. Dans cette matrice, la contribution de l'lment (4-6) est
localise en [4,4], [4,6], [6,4] et [6,6].
Le stockage de cette matrice ncessite la mmorisation de 10 x 10 =100 rels. Si l'on tient compte de
la symtrie de la matrice de rigidit, il suffit de mmoriser les termes situs sur la diagonale et au
dessus de la diagonale soit ( 10 x 11 ) / 2 = 55 rels.
Calculs
35
Lors de la factorisation de la matrice de rigidit, les termes nuls situs sur une colonne au dessus du
premier terme non nul n'interviennent pas. Il n'est pas utile de les stocker.
Les termes utiles de la matrice de rigidit sont mmoriss, colonne par colonne, dans une matrice
ligne [A] :
A(1)
A (2) A (3) A (5)
A ( 4) A ( 6)
A(7) A (8)
A (9)
[K ] =
A(10)
A (21)
A (11)
A(12)
A (22)
A (23)
A (13)
A(14)
A (24)
A (25)
A (15) A (16) A (18) A (26)

A(17) A(19) A(27)
A (20) A (28)
A (29)
A (30)
A (31)
A (32)
A (33)
A (34)
A(35)
A (36)
A (37)
A(38)
Cette mthode de stockage est appele mthode profil ou mthode ligne de ciel. Elle ncessite la
mmorisation de 38 rels.
Pour des problmes de grande taille, la mmoire centrale ne peut pas contenir toute la matrice [A].
Cette dernire est divise en blocs qui sont mmoriss sur le disque. Avec notre exemple, si on limite
la taille des blocs 15 rels, la matrice de rigidit [K] est mmorise sur le disque dans 3 blocs [A],
[B] et [C] :
A(1)
A (2) A (3) A (5)
A ( 4) A ( 6)
A(7) A (8)
A (9)
[K ] =
A(10)
B(6)
A (11)
A(12)
B(7)
B(8)
A (13)
A(14)
B(9)
B(10)
A (15) B(1) B(3) B(11)

B(2) B(4) B(12)
B(5) B(13)
B(14)
C(1)
C ( 2)
C(3)
C ( 4)
C(5)
C ( 6)
C (7 )
C(8)
C(9)
36
La place mmoire ncessaire au stockage de la matrice de rigidit [ K ] dpend de la numrotation des

noeuds.
Considrons l'exemple prcdent et renumrotons les noeuds :
8
5
1
10
Le stockage de la matrice de rigidit ncessite, avec cette numrotation, la mmorisation de 24 rels

dans deux blocs [ A ] et [ B ] :
A (3)
A (1)
A ( 2) A ( 4)
A (5) A(6)
A (7) A(8) A (10)
A (9) A (11)
[K ] =
A (12)
A (13)
B(1)
B(2)
B(3)
B(4)
B(5) B(8)
B(6) B(9)
B(7) B(10)
B(11)
La renumrotation des noeuds est automatique.
Bibliographie :
[B4] [D2] [I1] [Z3]
Remarque : entre la renumrotation des noeuds et la

segmentation en blocs, les d.d.l. sont renumrots (voir
partition des degrs de libert). Le premier d.d.l . impos
correspond toujours au dbut dun bloc. Chaque matrice
globale ( rigidit, masse... ) est donc segmente en NBL blocs
correspondant aux d.d.l. inconnus et NBI blocs correspondant
aux d.d.l. imposs :
Pour chaque matrice globale, seule la partie reprsente en

gris sur la figure ci-contre est mmorise.
1..NBL
1..NBI
Calculs
37
D - Rsolution dun systme dquations linaires

Soit rsoudre le systme dquations linaires :
[K][U] = [F]
o [K] est la matrice de rigidit, [F] le vecteur force et [U] le vecteur dplacement.
La mthode de rsolution utilise dans RDM est la mthode de Gauss.

Le calcul se fait en quatre tapes :
dcomposition de Crout de la matrice de rigidit globale : [K ] = [L] [D] [ L]T
o [L] est une matrice triangulaire infrieure dont les termes diagonaux sont gaux l'unit et
[D] est une matrice diagonale. Si le nombre de liaisons est suffisant tous les pivots sont
strictement positifs.
rsolution du systme triangulaire infrieur : [L] [Y] = [F]
rsolution du systme diagonal : [X ] = [D] 1 [Y]
rsolution du systme triangulaire suprieur : [L] T [ U ] = [X]
Bibliographie :
[B4] [D2] [I1] [Z3]
38
E - Recherche des valeurs et vecteurs propres

Prsentation
Rsoudre le problme aux valeurs propres consiste trouver les couples ( , [u]) qui satisfont
lquation :
[K][u] = [M][u] avec [u] [0] .
Les problmes traits sont les suivants :

Recherche des modes propres de vibration. Le problme rsoudre scrit :
[K ] [ U] = 2 [M] [ U]
[K] est la matrice de rigidit de la structure,

[M] est la matrice masse,
est une pulsation propre,
[u] est le vecteur des dplacements associ la pulsation .
Flambement linaire. Le problme rsoudre scrit :

[ K ] [ U ] = [ K ] [ U ]
[K] est la matrice de rigidit de la structure,

[K ] est la matrice de rigidit gomtrique,
est un coefficient de charge critique,
[u] est le vecteur des dplacements associ .
Remarque 1 : en fait, on rsout le problme aux valeurs propres :

[K LL ] [ U L ] = [M LL ] [ U L ] ( voir Partition des degrs de libert ).
Remarque 2 : les valeurs propres sont les solutions de lquation :

dterminant ( [K LL ] [M LL ] ) = 0.
Calculs
39
Remarque 3 : le nombre de modes propres demands ne doit pas excder le nombre de degrs de
libert excits. Plus prcisment, en labsence de modes rigides (la structure est isostatique ou
hyperstatique intrieurement et extrieurement), ce nombre est gal au rang de la matrice [M LL ] .
Considrons, par exemple, la structure suivante :

k
k
m
k
m
x
4
Les matrices [K LL ] et [M LL ] sont respectivement gales :

m 0 0
2k k 0
[K LL ] = k 2k k [ M LL ] = 0 m 0
0 0 m
0 k 2k
Cette structure possde autant de modes propres de vibration que de degrs de libert : 3. Les
pulsations propres sont gales :
1 = 0.7654
k
m
k
m
2 = 1.4142
3 = 1.8478
k
(rad/s)
m
Supprimons la masse porte par le noeud 3 :

k
k
m
k
m
x
4
Le nombre de degrs de libert et la matrice [K LL ] ne changent pas. Par contre la matrice [M LL ]

devient :
m 0 0
[M LL ] = 0 0 0
0 0 m
Le nombre de modes propres de cette structure est gal 2 et les pulsations propres sont gales :
1 =
k
m
2 = 1.4142
k
(rad/s)
m
40
Deux mthodes de calcul sont proposes : la mthode d'itration inverse et la mthode d'itration
sur sous-espace.
Mthode ditration inverse

Cette mthode est bien adapte au calcul de la plus petite valeur propre.
Lalgorithme est le suivant :
Choix dun vecteur de dpart : [u O ]
A la Ime itration, on effectue les opration suivantes :
Rsolution du systme linaire :
[K ] [u I ] = [M ] [u I1 ]
Calcul de I par le quotient de Rayleigh :

I =
[u I ] T [K ] [u I ]
avec a = [u I ] T [M] [u I ]
a
Mise jour du vecteur propre :

[u I ] =
[u I ]
a
I et [u I ] convergent respectivement vers la plus petite valeur propre et le vecteur propre [u]
associ.
Bibliographie :
[B4] [B5] [D2] [I1]
Calculs
41
Mthode ditration sur sous-espace

Cette mthode permet le calcul des m plus petites valeurs propres dun systme de grande
dimension.
Notre problme aux valeurs propres scrit alors :
[K][U] = [][M][U]
(1)
( 2)
o [ ] =
(m)

est la matrice diagonale des m plus petites valeurs propres
et [ U ] = [u (1) u ( 2 ) ... u ( m ) ] est la matrice des m vecteurs propres associs.
Lalgorithme est le suivant :

Choix de m vecteurs initiaux orthonoms :
[ U 0 ] = [u (01) u (02) ... u (0m ) ]
Les vecteurs de dpart sont gnrs de manire alatoire, orthogonaliss par la mthode de GramSchmidt puis norms.
A la Ime itration, on effectue les oprations suivantes :
Itration inverse simultane :
[K ] [X I ] = [M ] [ U I 1 ]
Analyse dinteraction :
- Calcul des matrices dinteraction :
[K I ] = [X I ] T [K ] [X I ] , [M I ] = [ X I ]T [M ] [X I ]
- Rsolution du problme aux valeurs propres rduit :

[ K I ] [ Q I ] = [ M I ] [Q I ] [ I ]
Ce problme est transform en sa forme classique. La matrice de rigidit rduite est

factorise par la mthode de Choleski :
[K I ] = [S]T [S]
42

Aprs multiplication des deux membres de lquation par [S] T , on obtient :
[D] [ Z] = [ Z] [ ] 1
en posant [ Z] = [S] [Q I ] et [D] = [S] T [M I ] [S]

Ce problme aux valeurs propres est rsolu par la mthode de Jacobi.
- Amlioration de lapproximation des vecteurs propres

[ U I ] = [ X I ] [Q I ]
Les matrices [ I ] et [ U I ] convergent respectivement vers la matrice des valeurs propres [] et la

matrice des vecteurs propres [U]. Les valeurs propres tant classes par ordre de module croissant, le
k-ime vecteur de la matrice [ U I ] converge vers kme vecteur propre avec un taux de convergence
gal k / m+1 .Il en rsulte que les premiers modes convergent plus vite que les modes dordre
lev. Pour amliorer la rapidit de convergence, il est donc intressant de choisir un sous-espace de
dimension q suprieur m. RDM utilise la valeur q = Min ( m+8 , 2 m ) recommande par [B4].
Dcalage du spectre de valeurs propres ( calcul des modes propres )

Remplaons le problme aux valeurs propres :
[K ] [ U] = 2 [M] [ U]
par :
( [K ] + 2 [M ] ) [ U] = ( 2 + 2 ) [M ] [ U] avec = 2 f
Les vecteurs propres du nouveau systme sont identiques ceux du systme initial. Les valeurs
propres 2 ont subi un dcalage gal : 2 .
Bibliographie :
[A6] [B4] [D2] [I1] [G6] [H4] [FLA 7] [FLA 9]
Maillage dun domaine plan 43
Maillage dun domaine plan
A - Triangulation de Delaunay
La mthode de Delaunay permet la triangulation dun domaine plan partir du maillage de sa
frontire en segments. La mise en oeuvre de cette mthode comprend les phases suivantes :
1 - Discrtisation de la frontire.
La frontire et ventuellement les lignes intrieures sont discrtises en segments.
2 - Cration de 4 points supplmentaires qui dfinissent un quadrilatre englobant le domaine

mailler. La discrtisation de cette bote en deux triangles constitue le maillage de dpart.
44
3 - Triangulation de la bote par insertion, un un, des points gnrs lors de phase 1 et
ventuellement de points supplmentaires imposs : points mailler, lignes intrieures...
Pn+1
On procde de la manire suivante :

Soit Tn la triangulation construite avec les n
premiers points et Pn+1 le point insrer dans le
maillage. Daprs la phase 2, ce point est interne
Tn. Soit E lensemble des triangles de Tn dont le
cercle circonscrit contient Pn+1. La construction
de la triangulation Tn+1 consiste :
- dtruire les triangles contenus dans E.
- remailler le sous-domaine dfini par E en
joignant Pn+1 aux artes externes de E.
4 - Elimination des triangles extrieurs au domaine.

5 - Cration des noeuds intrieurs au domaine.
Le maillage prcdent est en gnral impropre au calcul : il doit tre raffin. Une premire srie de
points est cre dans les triangles dont laire est plus grande quune certaine valeur. Ce processus est
rpt plusieurs fois. La mthode utilise est celle dcrite ci-dessus. La taille des lments gnrs
dpend :
- de la discrtisation de la frontire.
- du coefficient de dilution ( voir paramtres du maillage ).
- de la densit dlments autour des points mailler.
6 - Amlioration du maillage. Le maillage obtenu est souvent de qualit mdiocre. Il est amlior
par :
- lissage barycentrique : chaque noeud non situ sur la frontire est dplac au centre de gravit des
noeuds qui lui sont connects. Cette opration peut tre rpte plusieurs fois (voir paramtres du
maillage).
- suppression dartes.
- retournement dartes.
...
Bibliographie : [G5] [MAI 1] ... [MAI ...]
46
B - Qualit gomtrique du maillage

Plusieurs commandes de visualisation permettent de vrifier la qualit gomtrique du maillage. Les
critres utiliss sont :
Critre de distorsion :
Lmax
La qualit dun triangle est gale au rapport entre le rayon du cercle

inscrit multipli par 12 et la longueur du plus grand ct.
qualit = R 12 L max
La qualit dun quadrangle est gale au rapport entre la longueur du

plus petit ct multipli par 2 et la longueur de la plus grande
diagonale..
Dmax
qualit = C min 2 D max

Cmin
La qualit maximale est gale 1 ( triangle quilatral / carr ).
Critre du noeud milieu :

Ce critre est utilis pour les triangles 6 noeuds et les quadrangles 8 ou 9 noeuds.
L
Soit d la distance entre le noeud milieu dune arte et le segment de longueur L joignant
les noeuds extrmits de cette arte. La qualit dun lment est gale la valeur
maximale du rapport entre d et L sur llment :
qualit = max ( d / L )
La qualit maximale est gale 0 ( lment bords rectilignes ).
Critre du jacobien :
Ce critre est utilis pour les triangles 6 noeuds et les quadrangles 8 ou 9 noeuds.
A chaque lment rel du plan { x , y } est associ un lment de rfrence dans le plan { , }. La
transformation qui associe un point de llment de rfrence un point llment rel :
x = x (, )
y = y(, )
doit tre bijective.
+d
(0,1)
dA = det [J] d d
+d
+d
+d
(1,0)
Elment de rfrence
Elment rel
x

La matrice jacobienne [J] de la transformation : [J ] =
y

ne doit pas tre singulire et son dterminant doit conserver le mme signe sur tout llment. Compte
tenu des choix effectus par le logiciel ( sens trigonomtrique pour la numrotation des noeuds ) ce
dterminant doit tre positif.
La qualit dun lment est gale la valeur minimale du dterminant de la matrice jacobienne [J] sur
llment multiplie par laire de llment de rfrence et divise par laire de llment rel :
qualit = aire de llment de rfrence * ( min det [J] ) / aire de llment rel
La qualit maximale est gale 1.
48
Bibliographie
Site Internet :
J. MACKERLE, Information Retrieval on Finite Element Books : http://ohio.ikp.liu.se/fe/index.html
Ouvrages :
[A1]
[A2]
[A3]
[A4]
[A5]
[A6]
[A7]
[A8]
[A9]
[A10]
[A11]
[A12]
[B1]
[B2]
[B3]
[B4]
[B5]
[B6]
[B7]
[B8]
[B9]
[B10]
[B11]
[B12]
[B13]
[B14]
[B15]
P. AGATI, N. MATTERA, Mcanique applique, Dunod, 1990.

P. AGATI, N. MATTERA, Mcanique 2, Dunod, 1994.
J.E. AKIN, Finite elements for analysis and design, Academic Press, 1994.
M. ALBIGES, A. COIN, Rsistance des matriaux applique, 2 tomes, Eyrolles, 1976.
J. ARGYRIS, H.P. MLEJNEK, Die methode der finiten elemente, 3 tomes, Vieweg, 1986.
J. ARGYRIS, H.P. MLEJNEK, Dynamics of structures, North Holland, 1991.
J.F. AUBOUIN, Calcul des structures et informatique, Eyrolles, 1983.
B. ASANCHEYEV, Calcul des structures en BASIC, Eyrolles, 1984.
J. AVRIL, Encyclopdie danalyse des contraintes, Vishay-Micromesures.
P. AGATI, M. ROSSETTO, F. LEROUGE, Rsistance des matriaux, Dunod, 1999.
F. AXISA, Modlisation des systmes mcaniques, 4 tomes, Herms, 2001.
AISC, Load and Resistance Factor Design Specification for structural steel buildings,
AISC, 1999.
Y. BAMBERGER, Mcanique de lingnieur, 4 tomes, Hermann, 1997.
J.J. BARRAU, S. LAROZE, Calculs de structures par lments finis, Cours ENSAE, 1984.
K.J. BATHE, E.L. WILSON, Numerical methods in finite element analysis, Prentice Hall,
1976.
K.J. BATHE, Finite element procedures in engineering analysis, Prentice Hall, 1996.
J.L. BATOZ, G. DHATT, Modlisation des structures par lments finis - solides
lastiques, Herms, 1990.
J.L. BATOZ, G. DHATT, Modlisation des structures par lments finis - poutres et
plaques, Herms, 1990.
J.L. BATOZ, G. DHATT, Modlisation des structures par lments finis - coques, Herms,
1992.
A. BAZERGUI, T. BUI-QUOC, A. BIRON, G. McINTYRE, C. LABERGE, Rsistance des
matriaux, Editions de lEcole Polytechnique de Montral, 1993.
Z. BAZANT, L. CEDOLIN, Stability of structures, Oxford University Press, 1991.
D. BELLET, J.J. BARRAU, Cours dlasticit, Cpadus, 1990.
B. BOUMARD, F. LAVASTE, Rsistance des matriaux, Delagrave, 1984.
A.P. BORESI, R.J. SCHMIDT, O.M. SIDEBOTTOM, Advanced mechanics of materials,
Wiley, 1993.
J. BROZZETTI, M.A. HIRT, R. BEZ, Trait du gnie civil (Complment au volume 10) Construction mtallique - Exemples numriques adapts aux Eurocodes, Presses
Polytechniques et Universitaires Romandes, 1995.
R.D. BLEVINS, Formulas for natural frequency and mode shape, Krieger, 1993.
G. BEER, J.O. WATSON, Introduction to finite and boundary element methods for
engineers, Wiley, 1994.
Bibliographie
[B16]
[B17]
[B18]
[B19]
[B20]
[B21]
[B22]
[B23]
[B24]
[C1]
[C2]
[C3]
[C4]
[C5]
[C6]
[C7]
[C8]
[C9]
[C10]
[C11]
[C12]
[C13]
[C14]
[C15]
[C16]
[C17]
[C18]
[C19]
[C20]
[C21]
[C22]
[C23]
[C24]
[C25]
[C26]
49
A. BAZERGUI, T. BUI-QUOC, A. BIRON, G. McINTYRE, C. LABERGE, Rsistance des

matriaux- recueil de problmes - 2 tomes, Editions de lEcole Polytechnique de Montral,
1991.
K.J. BATHE, Finite element procedures, Prentice Hall, 1996.
J.J. BARRAU, S. LAROZE, Mcanique des structures Tome 1 bis : Exercices Solides
lastiques , plaques et coques, ENSAE, 1989.
R. BOUDET, P. STEPHAN, Rsistance des matriaux, Cpadus, 1998.
J. BESSON, G. CAILLETAUD, J.L. CHABOCHE, S. FOREST, Mcanique non linaire
des matriaux, Cours IPSI, 1997.
Ph. BISCH, D. BROCHARD, P. BULAND, F. GANTEBEIN, P. LABB, Comportement
sismique des btiments et quipements, Cours IPSI, 1998.
T. BELYTSCHKO, W. K. LIU, B. MORAN, Non linear finite elements for continua and
structures, Wiley, 2000.
H.D. BUI, Introduction aux problmes inverses en mcanique des matriaux, Eyrolles,
1993.
J. BESSON, G. CAILLETAUD, J.-L. CHABOCHE, S. FOREST, Mcanique non linaire
des matriaux, Herms, 2001.
G. CAIGNART, J.P. HENRY, Exercices dlasticit, Dunod, 1976.
J.A. CALGARO, M. VIRLOGEUX, Projet et construction des ponts - Analyse structurale
des tabliers de ponts, Presses de lEcole nationale des Ponts et Chausses, 1994.
Centre Scientifique et Technique du Btiment, Rgles de conception et de calcul des
charpentes en alliages daluminium - Rgles AL, Juillet 1976.
W.F. CHEN, E.M. LUI, Stability design of steel frames, CRC Press, 1991
L. CHEVALIER, Mcanique des systmes et des milieux dformables, Ellipses, 1996.
P. CHILLON, M. KERGUINAS, Rsistance des matriaux, Dunod, 1969.
W. CLOUGH, J. PENZIEN, Dynamique des structures, Pluralis, 1980.
Construire en acier, Le Moniteur, 1993.
R.D. COOK, D.S. MALKUS, M.E. PLESHA, Concepts and applications of finite element
analysis, Wiley, 1989.
J. COURBON, Calcul des structures, Techniques de lingnieur (A300).
J. COURBON, Rsistance des matriaux, Techniques de lingnieur (C2000).
J. COURBON, Formulaire de rsistance des matriaux, Techniques de lingnieur (C2060).
J. COURBON, Elments de rsistance des matriaux, Dunod, 1970.
J.C. CRAVEUR, Modlisation des structures, Masson, 1996.
M.A. CRISFIED, Finite elements and solution procedures for structural analysis, Pineridge
Press, 1986.
M.A. CRISFIED, Non-linear finite element analysis of solids and structures - 2 volumes,
Wiley, 1991,1997.
A. CURNIER, Mthodes numriques en mcanique des solides, Presses Polytechniques et
Universitaires Romandes, 1993.
R. CRAVERO, Bases pour la rsistance des matriaux, Ellipses, 1997.
L. CHEVALIER, Exercices et problmes corrigs de Mcanique des systmes et des milieux
dformables, Ellipses, 1997.
J.A. CALGARO, Poutres parois minces, Hrmes, 1998.
R. CRAVERO, Exercices pour la rsistance des matriaux, Ellipses, 1998.
R.D. COOK, Finite element modeling for stress analysis, Wiley, 1995.
A. CARPINTERI, Structural mechanics A unified approach, Chapman & Hall, 1997.
Construction mtallique et mixte acier-bton, 2 volumes, Eyrolles, 1996.
G. CAILLETAUD, J.Y. COGNARD, C. CORNUAULT, G. DHATT, Mcanique non
linaire des solides et des structures, Cours IPSI, 1995.
Physique et mcanique de la mise en forme des matriaux Ecole dt dOlron, Presses
du CNRS / IRSID, 1990.
50
[C27]
[C28]
[C29]
[C30]
[D1]
[D2]
[D3]
[D4]
[D5]
[D6]
[D7]
[D8]
[D9]
[D10]
[D11]
[D12]
[D13]
[E1]
[E2]
[E3]
[E4]
[F1]
[F2]
[F3]
[F4]
[F5]
[F6]
[F7]
[F8]
[F9]
[F10]
[F11]
[G1]
[G2]
[G3]
[G4]
[G5]
A. CAPRA, V. DAVIDOVICI, Calcul dynamique des structures en zone sismique, Eyrolles,

1982.
G. CIARLET, Introduction lanalyse numrique matricielle et loptimisation, Dunod,
1998.
J. COIRIER, Mcanique des milieux continus, Dunod, 2001.
J. COIRIER, Mcanique des milieux continus Aide mmoire, Dunod, 2001.
D. DARTUS, Elasticit linaire, Cpadus, 1995.
G. DHATT, G. TOUZOT, Une prsentation de la mthode des lments finis, Maloine
1984.
R. DOUMERG, Rsistance des matriaux, Edition de lauteur, 1967.
B. DROUIN, J.M. SENICOURT, F. LAVASTE, G. FEZANS, De la mcanique vibratoire
classique la mthode des lments finis, 2 volumes, AFNOR, 1993.
J. DUC, D. BELLET, Problmes dlasticit, Cpadus, 1984.
G. DUVAUT, Mcanique des milieux continus, Masson, 1990.
S. DUBIGEON, Mcanique des milieux continus, Technique et documentation, 1998.
H. DUMONTET, G. DUVAUT, F. LN, P. MULLER, N. TURB, Exercices corrigs de
mcanique des milieux continus, Dunod, 1998.
J. DOUIN, Mcanique des milieux continus, Diderot, 1997.
D.DESJARDINS, M. TOUZET-CORTINA, Introduction la mcanique des milieux
continus, Dunod, 1999.
D. DIDIER, M. Le BRAZIDEC, P. NATAF, R. PRALAT, G. SIMON, J. THIESSET,
Structures de gnie civil, AFNOR/Nathan, 1998.
V. DAVIDOVICI, Gnie Parasismique, Presses de lEcole Nationale des Ponts et
Chausses, 1985.
M. DUBESSET, Le manuel du Systme International dunits Lexique et conversions,
Technip, 2000.
B. EYGLUMENT, Thermique thorique et pratique, Herms, 1994.
Eurocode 3 Calcul des structures en acier , AFNOR, 1992.
D. EUVRARD, Rsolution numrique des quations aux drives partielles, Masson, 1994.
D.J. EWINS, Modal testing : theory and practice, Wiley, 1986.
M.J. FAGAN, Finite element analysis - Theory and practice, Longman, 1992.
V. FEODOSSIEV, Rsistance des matriaux, Editions Mir, 1971.
F. FREY, Trait du gnie civil (volume 1) - Analyse des structures et milieux continusStatique applique, Presses Polytechniques et Universitaires Romandes, 1990.
F. FREY, Trait du gnie civil (volume 2) - Analyse des structures et milieux continusMcanique des structures, Presses Polytechniques et Universitaires Romandes, 2000.
F. FREY, Trait du gnie civil (volume 10) - Constructions mtalliques, Presses
F. FREY, Trait du gnie civil (volume 3) - Analyse des structures et milieux continusMcanique des solides, Presses Polytechniques et Universitaires Romandes, 1998.
A. FILIATRAULT, Elments de gnie parasismique et de calcul dynamique des structures,
Editions de lEcole Polytechnique de Montral, 1996.
P. J. FREY, P.L. GEORGE, Maillages, Herms, 1999.
F. FREY, J. JIROUSEK, Trait du gnie civil (volume 6) Mthode des lments finis,
Presses Polytechniques et Universitaires Romandes, 20010.
D. FRANOIS, Essais mcaniques et lois de comportement, Herms, 2001.
J.-L. FANCHON, Guide de mcanique, Nathan, 2001.
R.H. GALLAGHER, Introduction aux lments finis, Pluralis, 1976.
R.H. GALLAGHER, W. McGUIRE, Matrix structural analysis, Wiley, 1979.
D. GAY, J. GAMBELIN, Une approche simple du calcul des structures par la mthode des
lments finis, Herms, 1989.
D. GAY, Matriaux composites, Herms, 1997.
P.L. GEORGE, Gnration automatique de maillages, Masson, 1991.
Bibliographie
[G6]
[G7]
[G8]
[G9]
[G10]
[G11]
[G12]
[G13]
[G14]
[G15]
[G16]
[G17]
[G18]
[G19]
[G20]
[G21]
[H1]
[H2]
[H3]
[H4]
[H5]
[H6]
[H7]
[H8]
[H9]
[I1]
[I2]
[I3]
[J1]
[J2 ]
[K1]
[K2]
[L1]
[L2]
[L3]
[L4]
[L5]
[L6]
51
M GERADIN, D. RIXEN, Thorie des vibrations, Masson, 1996.

P. GERMAIN, Mcanique, 2 tomes, Ellipses, 1986.
P. GERMAIN, P. MULLER, Introduction la mcanique des milieux continus, Masson,
1994.
A. GIET, L. GEMINARD, Rsistance des matriaux, 2 tomes, Dunod, 1972.
J.E. GORDON, Structures et matriaux, Belin, 1994.
P.L. GEORGE, H. BOROUCHAKI, Triangulation de Delaunay et maillage, Herms, 1997.
J.M. GERE, S.P. TIMOSHENKO, Mechanics of materials, Chapman & Hall, 1995.
T. GMR, Dynamique des structures, Presses Polytechniques et Universitaires Romandes,
1997.
G. GENTA, Vibration of structures and machines, Springer Verlag, 1995.
D. GAY, J. GAMBELIN, Dimensionnement des structures, Herms, 1999.
D. GUITARD, Mcanique du matriau bois et composites, Cpadus, 1987.
R. J. GIBERT, Vibrations des structures, Eyrolles (Collection DER EDF n 66), 1988.
Jauges extensomtriques, Groupement pour lAvancement des Mthodes dAnalyse des
Contraintes, Technique & Documentation, 2000.
M. GERADIN, A. CARDONA, Flexible Multibody Dynamics A Finite Element Approach,
Wiley, 2001.
Y. GOURINAT, Introduction la dynamique des structures, Cpadus, 2001.
P.-G. GEORGE, Maillage et adaptation, Herms, Lavoisier, 2001.
J.P. HENRY, Cours dlasticit, Dunod, 1982.
E. HINTON, D.R.J. OWEN, An introduction to the finite element computations, Pineridge
Press, 1985.
M.A. HIRT, R. BEZ, Trait du gnie civil (volume 10) - Construction mtallique, Presses
T. HUGHES, The finite element method - Linear static and dynamic finite element analysis,
Prentice Hall, 1987.
J.P. DEN HARTOG, Advanced strength of materials, Dover, 1987.
D. HENWOOD, J. BONET, Finite elements A gentle introduction, Macmillan, 1996.
K.H. HUEBNER, E.A. THORNTON, T.G. BYRON, The finite element method for
engineers, Wiley, 1995.
C.M. HARRIS, Shock and vibration handbook, McGraw-Hill, 1995.
M.A. HIRT, M. CRISINEL, Trait du gnie civil (volume 11) - Charpentes mtalliques,
Presses Polytechniques et Universitaires Romandes, 2001.
J.F. IMBERT, Analyse des structures par lments finis, Cpadus, 1995.
MODULEF : Une bibliothque modulaire dlments finis, INRIA, 1988.
F.P. INCROPERA, D.P. DEWITT, Fundamentals of Heat and Mass Transfer, Wiley, 1996.
W.A. JALIL, Calcul pratique des structures, Eyrolles, 1985.
P. JOLY, Mise en uvre de la mthode des lments finis, Ellipes, 1990.
N. KIKUCHI, Finite element methods in mechanics, Cambridge University Press, 1986.
I.A. KARNOVSKY, O.I. LEBED, Formulas for structural dynamics, McGraw-Hill, 2000.
M. LALANNE, P. BERTHIER, J. D. HAGOPIAN, Mcanique des vibrations linaires,
Masson, 1986.
L. LANDAU, E. LIFCHITZ, Thorie de llasticit, Mir, 1967.
S. LAROZE, J.J. BARRAU, Mcanique des structures - Dynamique des structures Contraintes dorigine thermique, Eyrolles/Masson, 1992.
S. LAROZE, J.J. BARRAU, Mcanique des structures - Calcul des structures en matriaux
composites, Eyrolles/Masson, 1987.
S. LAROZE, J.J. BARRAU, Mcanique des structures - Solides lastiques - Plaques et
coques, Eyrolles/Masson, 1988.
S. LAROZE, Mcanique des structures - Thorie des poutres, Eyrolles/Masson, 1988.
52
[L7]
[L8]
[L9]
[L10]
[L11]
[L12]
[L13]
[L14]
[L15]
[L16]
[L17]
[L18]
[M1]
[M2]
[M3]
[M4]
[M5]
[M6]
[M7]
[M8]
[M9]
[M10]
[M11]
[M12]
[N1]
[N2]
[N3]
[O1]
[O2]
[O3]
[O4]
[O5]
[P1]
[P2]
[P3]
[P4]
[P5]
[P6]
[P7]
[P8]
P. LASCAUX, R. THEODOR, Analyse numrique matricielle applique lart de

lingnieur, 2 tomes, Masson, 1993.
R.W. LEWIS, K. MORGAN, H.R. THOMAS, K.N. SEETHARAMU, The finite element
method in heat transfer analysis, Wiley - 1996.
B. LUCQUIN, O. PIRONNEAU, Introduction au calcul scientifique, Masson, 1996.
S. LAROZE, M. LORRAIN, Mcanique des structures Tome 2 bis : Exercices Poutres,
ENSAE, 1989.
R.K. LIVESLEY, Finite elements : an introduction for engineers, Cambridge University
Press, 1983.
A.E. LOVE, A treatise on the mathematical theory of elasticity, 1944.
M. LAREDO, Rsistance des matriaux, Dunod, 1970.
A. LEISSA, Vibration of plates, Acoustical Society of America, 1993.
A. LEISSA, Vibration of shells, Acoustical Society of America, 1993.
A. LE POURHIET, Rsolution numrique des quations aux drives partielles, Cpadus,
1988.
J. LEMAITRE, J.L. CHABOCHE, Mcanique des matriaux solides, Dunod, 1996.
P. LADEVZE, J.-P. PELLE, La matrise du calcul en mcanique linaire et non linaire,
Herms, 2001.
C. MASSONET, Rsistance des matriaux, 2 tomes, Dunod, 1968.
C. MASSONNET, G. DEPREZ, R. MAQUOI, R. MULLER, G. FONDER, Calcul des
structures sur ordinateur, Tomes 1 et 2, Eyrolles Masson, 1972.
C. MASSONNET, S. CESCOTTO, Mcanique des matriaux, De Bck Universit, 1994.
G.A. MOHR, Finite elements for solids, fluids and optimisation, Oxford Science
Publications, 1992.
J. MOREL, Calcul des structures mtalliques selon lEurocode 3, Eyrolles, 1994.
J. MOREL, Structures mtalliques - CM 66 additif 80 - Eurocode 3, Eyrolles, 1995.
P. MAITRE, Formulaire de la construction mtallique, Le Moniteur, 1997.
R.H. MACNEAL, Finite elements : their design and performance, Dekker, 1994.
P. MULLER, C. OSSADZOW, Introduction aux coques minces lastiques, Herms, 1999.
J.P. MERCIER, G. ZAMBELLI, W. KURZ, Introduction la science des matriaux,
Presses Polytechniques et Universitaires romandes, 1999.
L. MEIROVITCH, Elements of vibration analysis, Mc Graw Hill, 1986.
L. MEIROVITCH, Principles and techniques of vibrations, Prentice Hall, 1997.
W.A. NASH, Rsistance des matriaux, McGraw Hill, 1995.
J.P. NOUGIER, Mthodes de calcul numrique, Masson, 1991.
J-C. NEDELEC, Notions sur les techniques dlments finis, Ellipses, 1991.
M. OJALVO, Thin-walled bars with open profiles, Olive Press, 1990.
N.OTTOSEN, H. PETERSSON, Introduction to the finite element method, Prentice Hall,
1992.
Lexique de construction mtallique et de rsistance des matriaux, OTUA.
J. OUIN, Mcanique des structures, Casteilla, 1997.
J. OUIN, Transferts thermiques, Casteilla, 1998.
W.D. PILKEY, Fomulas for stress,strain and structural matrices, Wiley, 1994.
G. PISSARENKO, A. YAKOVLEV, V MATVEEV, Aide-mmoire de rsistance des
matriaux, Editions Mir, 1985.
J. PLUSQUELLEC, Vibrations, Techniques de lingnieur (A410).
G. PLUVINAGE, Mcanique lastoplastique de la rupture, Cpadus, 1989.
G. PLUVINAGE, 120 exercices de mcanique lastoplastique de la rupture, Cpadus,
1989.
M. PRAT, La modlisation des ouvrages, Herms, 1995.
J.S. PRZEMIENIECKI, Theory of matrix structural analysis, Dover, 1986.
M. PETYT, Introduction to finite element vibration analysis, Cambridge University Press,
1990.
Bibliographie
[P9]
[Q1]
[R1]
[R2]
[R3]
[R4]
[R5]
[R6]
[R7]
[R8]
[R9]
[R10]
[R11]
[R12]
[R13]
[R14]
[S1]
[S2]
[S3]
[S4]
[S5]
[S6]
[S7]
[S8]
[S9]
[S10]
[T1]
[T2]
[T3]
[T4]
[T5]
[T6]
[T7]
[T8]
[T9]
[T10]
[T11]
[V1]
[W1]
53
J.P. PELLE, P. BECKERS, L. GALLIMARD, Estimations des erreurs de discrtisation et

analyses adaptatives Application lautomatisation des calculs lments finis, Cours IPSI,
1996.
Q.S. NGUYEN, Stabilit et mcanique non linaire, Herms, 2000.
S.S. RAO, The finite element method in engineering, Pergamon, 1989.
N. RECHO, Rupture par fissuration des structures, Herms - 1995.
J.N. REDDY, An introduction to the finite element method, McGraw-Hill, 1993.
Rgles de construction parasismique - Rgles PS applicables aux btiments - PS92,
Eyrolles.
V. REKATCH, Problmes de la thorie de llasticit, Mir, 1980.
J.P. REZETTE, F. LELEUX, Calcul dynamique des structures par la mthode des lments
finis, CETIM, 1974.
J.J. RISLER, Mthodes mathmatiques pour la CAO, Masson, 1991.
K.C. ROCKEY, H.R. EVANS, D.W. GRIFFITHS, Elments finis, Eyrolles, 1979.
J. ROUX, Rsistance des matriaux, 3 tomes, Eyrolles, 1995.
S.S. RAO, Mechanical vibrations, Addison Wesley, 1995.
M. RAPPAZ, M. BELLET, M. DEVILLE, Modlisation numrique en science et gnie des
matriaux, Presses Polytechniques et Universitaires Romandes, 1998.
M. RAOUS, S. CESCOTTO, A. CURNIER, A. MILLARD, Modlisation mcanique et
numrique du contact et du frottement, Cours IPSI, 1996.
M. ROSEAU, Vibration des systmes mcaniques, Masson, 1984.
M. RAPPAZ, M. BELLET, M. DEVILLE, Modlisation numrique en science et gnie des
matriaux, Presses Polytechniques et Universitaires romandes, 1998.
J. SALENON, Mcanique du continu - Concepts gnraux, Ellipses, 1996.
J. SALENON, Mcanique du continu - Thermolasticit , Ellipses, 1996.
J. SALENON, Mcanique du continu - Milieux curvilignes, Ellipses, 1996.
I.M. SMITH, D.V. GRIFFITHS, Programming the finite element method, Wiley, 1988.
Socit Franaise des Mcaniciens, Guide de validation des progiciels de calcul de
structures, AFNOR, 1990.
N.Q. SON, Stabilit des structures lastiques, Springer-Verlag, 1995.
B. SZABO, I. BABUSKA, Finite element analysis, Wiley, 1991.
M.A. STUDER, F. FREY, Introduction l'analyse des structures, Presses Polytechniques et
Universitaires Romandes, 1997.
L. SAINSAULIEU, Calcul scientifique, Dunod, 2000.
C. SOIZE, Dynamique des structures, Ellipses, 2001.
S.P. TIMOSHENKO, Rsistance des matriaux, Tome 1, Dunod, 1968.
S.P. TIMOSHENKO, Rsistance des matriaux, Tome 2, Dunod, 1968.
S.P. TIMOSHENKO, J.M. GERE, Theory of elastic stability, McGraw Hill, 1961.
S.P. TIMOSHENKO, J.M. GERE, Thorie de la stabilit lastique, Dunod, 1966.
S.P. TIMOSHENKO, J.N. GOODIER, Theory of elasticity, McGraw Hill, 1970.
S.P. TIMOSHENKO, J.N. GOODIER, Thorie de llasticit, Branger, 1968.
S.P. TIMOSHENKO, History of strength of materials, Dover.
N.S. TRAHAIR, M.A. BRADFORD, The behavior and design of steel structures, Chapman
& Hall, 1991.
N.S. TRAHAIR, Flexural-torsional buckling of structures, Chapman & Hall, 1993.
P. TROMPETTE, Mcanique des structures par la mthode des lments finis, Masson,
1992.
P. TROMPETTE, C. FLEURY, C. KNOPF-LENOIR , Optimisation des structures, Cours
IPSI, 1987.
R. VALID, La mcanique des milieux continus et le calcul des structures, Eyrolles, 1977.
R. WALTHER, B. HOURIET, W. ISLER, P. MOA, Ponts haubanns, Presses
Polytechniques et universitaires romandes, 1995.
54
[W2]
[W3]
[W4]
[Y1]
[Z1]
[Z2]
[Z3]
[Z4]
W. WEAVER, S.P. TIMOSHENKO, D.H. YOUNG, Vibration problems in engineering,

Wiley, 1990.
A.L. WINDOW, G.S. HOLISTER, Strain gauge technology, Elsevier Applied Science,
1982.
C. WIELGOSZ, Rsistance des matriaux, Ellipses, 1999.
W.C. YOUNG, Roarks formulas for stress and strain, McGraw Hill, 1989.
Z.Z. ZHONG, Finite element procedures for contact-impact problems, Oxford Science
Publications, 1993.
H. ZIEGLER, Principles of structural stability, 2nd dition, Birkauser Verlag, 1977.
O.C. ZIENKIEWICZ, R.L. TAYLOR, La mthode des lments finis - Formulation de base
et problmes linaires, AFNOR, 1989.
O.C. ZIENKIEWICZ, R.L. TAYLOR, The finite element method - Solid and fluid
mechanics, Dynamics and non-linearity, McGraw Hill, 1991.
Compte rendus de confrences :

[CETIM 1994] Conception-calcul : nouveaux outils pour les bureaux dtudes, CETIM, 1994.
[CETIM 1995] Dimensionnement au flambement, CETIM, 1995.
[CTICM 1983] Stabilit des structures, CTICM, 1983.
[FEMCAD 1987] Optimisation des structures et CAO, Herms, 1987.
[GIENS 1993-2001] Colloque national en calcul des structures - Giens (Var), Herms, 1993, 1995,
1997, 1999, 2001.
[STRUCOME 1988-1995] STRUCOME - Actes du congrs, DATAID, 1988-95.
Bibliographie
55
Publications :
Section droite : caractristiques et contraintes :
[SEC 1]
[SEC 2]
[SEC 3]
[SEC 4]
[SEC 5]
[SEC 6]
[SEC 7]
[SEC 8]
[SEC 9]
[SEC 10]
[SEC 11]
[SEC 12]
[SEC 13]
[SEC 14]
[SEC 15]
[SEC 16]
G.R. COWPER, The shear coefficient in Timoshenkos beam theory, Journal of Applied
Mechanics, 335-340 (June 1966).
C. CZEKAJSKI, A. POTIRON, D. GAY, TORFLEX : Calcul des constantes
caractristiques de flexion-torsion dune section de poutre par quations intgrales,
[FEMCAD 1987] p. 258.
Z. FRIEDMAN, J.B. KOSMATKA, Torsion and flexure of a prismatic isotropic beam
using the boundary element method, Computers & Structures, vol. 74, 479-494 (2000).
F. GRUTTMANN, R. SAUER, W. WAGNER, Shear stresses in prismatic beams with
arbitrary cross-sections, IJNME, vol. 45, 865-889 (1999).
L.R. HERMANN, Elastic torsional analysis of irregular shapes, Journal of the
Engineering Mechanics Division, ASCE, vol. 91, 11-19 (1965).
M.A. JASWON, A.R. PONTER, An integral equation solution of the torsion problem,
Proceedings, Royal Society of London A, vol. 273, 237-246 (1963).
W.E. MASON, L.R. HERMANN, Elastic shear analysis of general prismatic beams,
Journal of the Engineering Mechanics Division, ASCE, vol. 91, 965-982 (1968).
S.H. NGUYEN, An accurate finite element formulation for linear elastic torsion
calculations, Computers & Structures, vol. 42, 707-711 (1992).
T. NOURI, D. GAY, J.M. CIEAUX, Calcul par lments finis des caractristiques
homognises de flexion et de torsion dans une poutre composite, [STRUCOME 1991]
p. 870.
T. NOURI, D. GAY, J.M. CIEAUX, Homognisation et contraintes de cisaillement
dans une poutre composite phases orthotropes, Revue des composites et des matriaux
avancs, vol. 2, n2, 165-181 (1992).
J. SAPOUNTZAKIS, Solution of non-uniform torsion of bars by an integral equation
method, Computers & Structures, vol. 77, 659-667 (2000).
U. SCHRAMM, V. RUBENCHIK, W. PILKEY, Beam stiffness matrix based on the
elasticity equations, IJNME, vol. 40, 211-232 (1997).
K.S. SURANA, Isoparametric elements for cross-sectional properties and stress
analysis of beams, IJNME, vol. 14, 475-497 (1979).
W. WAGNER, F. GRUTTMANN, Finite element analysis of Saint-Venant torsion
problem with exact integration of the elastic-plastic constitutive equations, CMAME 190
(2001) 3831-3848.
W. YEIH, J.R. CHANG, J.T. CHEN, A boundary formulation for calculating moments of
an arbitrary closed planar region, Engineering analysis with boundary elements, vol. 23,
611-617 (1999).
C.H. YOO, S.W. ACRA, Cross-sectional properties of thin-walled multi-cellular
section, Computers & Structures, vol. 22, 53-61 (1986).
56

Flambement :
[FLA 1]
[FLA 2]
[FLA 3]
[FLA 4]
[FLA 5]
[FLA 6]
[FLA 7]
[FLA 8]
[FLA 9]
[FLA 10]
[FLA 11]
[FLA 12]
[FLA 13]
[FLA 14]
[FLA 15]
[FLA 16]
[FLA 17]
[FLA 18]
J.H. ARGYRIS, P.C. DUNNE, D.W. SCHARFF, On large displacement-small strain

analysis of structures with rotational degrees of freedom, CMAME, vol 14, 401-451
(1978).
J.H. ARGYRIS, P.C. DUNNE, G.A. MALEJANNAKIS, D.W. SCHARFF, On large
displacement-small strain analysis of structures with rotational degrees of freedom,
CMAME, vol 15, 99-135 (1978).
J.H. ARGYRIS, H. BALMER, J.St. DOLTSINIS, P.C. DUNNE, M. HAASE, M.
KLEIBER, G.A. MALEJANNAKIS, H.P. MLEJNEK, M. MLLER, D.W. SCHARFF,
Finite element method - the natural approach, CMAME, vol 17/18, 1-106 (1979).
J.H. ARGYRIS, O. HILPERT, G.A. MALEJANNAKIS, D.W. SCHARFF, On the
geometrical stiffness of a beam in space - a consistant v.w. approach, CMAME, vol 20,
105-131 (1979).
J.H. ARGYRIS, An excursion into large rotations, CMAME, vol 32, 85-155 (1982).
J. ARGYRIS, V.F. POTERASU, Large rotations revisited, application of Lie algebra,
CMAME, vol 103, 11-42 (1993).
M. DUBOIS, Les mthodes numriques applicables au dimensionnement vis--vis du
flambement, [CETIM 1995], pages 1-53.
Y.Q. GUO, Analyse non linaire statique et dynamique des poutres tridimensionnelles
lasto-plastiques, Thse de doctorat, Universit de Compigne, 1987.
M. KNIG, D. NAGY, P. STREINER, Buckling analysis with the ASKA program
system, CMAME, vol 16, 185-212 (1978).
R. KOUHIA, On kinematical relations of spatial framed structures, Computers and
Structures, vol 40, N 5, 1185-1191 (1991).
K. KOUHIA, M. TUOMALA, Static and dynamic analysis of space frames using simple
Timoshenko type elements, IJNME, vol 36, N 5, 1189-1221 (1993).
A.F. SALEEB, T.Y.P. CHANG, A.S. GENDY, Effective modelling of spatial buckling of
beam assemblages, accounting for warping constraints and rotation-dependency of
moments, IJNME, vol 33, 469-502 (1992).
H. SHAKOURZADEH, Y.Q. GUO, J.L. BATOZ, P. PENSIRINI, Modlisation du
comportement non linaire de structures formes de poutres 3D , [GIENS 93] pp.71-81.
H. SHAKOURZADEH, Modlisation des structures poutres tridimensionnelles parois
minces et simulation du comportement non linaire gomtrique et lasto-plastique,
Thse de doctorat, Universit de Compigne, 1994.
H. SHAKOURZADEH, Y.Q. GUO, J.L. BATOZ, Sur la stabilit de poutres 3D parois
minces, [GIENS 95] pp. 335-340.
Y.B. YANG, W.McGUIRE, Stiffness matrix for nonlinear analysis, J. Struct. Eng.,
ASCE, vol 112, N 4, 853-877, 1986.
Y.B. YANG, W. McGUIRE, Joint rotation and geometric non linear analysis, J. Struct.
Eng., ASCE, vol 112, N 4, 879-905, 1986.
Y.B. YANG, S.H. KUO, Frame buckling analysis with full consideration of joint
compatibilities, J. Struct. Eng., ASCE, vol 118, N 5, 871-889, 1992.
Bibliographie
57
Maillage :
[F8] [G5] [G11]
[MAI 1]
[MAI 2]
[MAI 3]
[MAI 4]
[MAI 5]
[MAI 6]
[MAI 7]
[MAI 8]
[MAI 9]
[MAI 10]
[MAI 11]
[MAI 12]
[MAI 13]
[MAI 14]
[MAI 15]
[MAI 16]
H. BOROUCHAKI, P.L. GEORGE, Encore Delaunay, Rapport de Recherche N 2461,

INRIA, 1995.
H. BOROUCHAKI, P.L. GEORGE, Quelques conclusions sur la construction dun
mailleur 2D de type Delaunay, Rapport de Recherche N 2463, INRIA, 1995.
W.H. FREY, Selective refinement : a new strategy for automatic node placement in
graded triangular meshes, IJNME, vol 24, 2183-2200 (1987).
P.L. GEORGE, Gnration de maillage par une mthode de type Vorono Rapport de
Recherche N 1398, INRIA, 1991.
P.L. GEORGE, F. HECHT, M.G. VALLET, Creation of internal points in Voronois
type method. Control adaptation, Adv. Eng. Software, vol 13, N 5/6, 303-312, 1991.
P.L. GEORGE, F. HECHT, E. SALTEL, Automatic mesh generator with specified
boundary, CMAME, vol 92, 269-288, 1991.
P.L. GEORGE, F. HERMELINE, Delaunays mesh of a convex polyhedron in dimension
d. Application to arbitrary polyhedra, IJNME, vol 33, 975-995, 1992.
P.L. GEORGE, Gnration automatique de maillages, Mcanique Industrielle et
Matriaux, vol 46 , N 2, 1993.
P.L. GEORGE, Synthse de quelques mthodes de gnration de maillages Elments
finis, Revue europenne des lments finis, vol 2, 121-153, 1993.
P.L. GEORGE, Les problmes de ladaptation sur lexemple dun mailleur 2D soumis
une carte de taille, [STRUCOME 1994], pages 101-116.
P.L. GEORGE, Quelques mthodes de gnration de maillages, [CETIM 1994], pages
33-54.
F. HERMELINE, Triangulation automatique dun polydre de dimension N, RAIRO,
Analyse numrique, vol 16 (3) , pages 211-242, 1982.
E.B. de LISLE, P.L. GEORGE, Optimisation de maillages tridimensionnels Rapport de
Recherche N 2046, INRIA, 1993.
D.F. WATSON, Computing then n-dimensional Delaunay tessellation with application
to Voronoi polytopes, The computer journal, vol 24, n2, 167-172, 1981.
O.C. ZIENKIEWICZ, D.V. PHILLIPS, An automatic mesh generation scheme for plane
and curved surfaces by isoparametric co-ordinates, IJNME, vol 3, 519-528, 1971.
H. BOROUCHAKI, P.L. GEORGE, Aspects of 2-D Delaunay mesh generation, IJNME,
vol 40, 1957-1975, 1997.
Plaques :
[PLA 1]
[PLA 2]
[PLA 3]
J. L. BATOZ, An explicit formulation for an efficient triangular plate-bending element,

IJNME, Vol 18, 1077-1089 (1982).
J.L. BATOZ, P. LARDEUR, A discrete shear triangular nine d.o.f. element for the
analysis of thick to very thin plates, IJNME, Vol. 28, 533-560 (1989).
K.J. BATHE, E. DVORKIN, A four-node plate bending element based on MindlinReissner plate theory and a mixed interpolation, IJNME, Vol. 21, p 367-383, (1985).
Abrviations :
CMAME
IJNME
J. Struct. Eng.
IPSI
Computer Methods in Applied Mechanics and Engineering

International Journal for Numerical Methods in Engineering
Journal of Structural Engineering
Institut pour la Promotion des Sciences de lIngnieur

Rdmtheo

Transféré par

Droits d'auteur :

Formats disponibles

Rdmtheo

Transféré par

Informations du document

Description originale:

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

Rdmtheo

Transféré par

Droits d'auteur :

Formats disponibles

RDM

Institut Universitaire de Technologie

Maillage dun domaine plan : triangulation de Delaunay

A - Contraintes autour dun point

Nous appellerons " facette n en M" cet lment de surface.

RDM - Complments thoriques

est une valeur algbrique positive ( traction ) ou ngative ( compression ).

II - Etat de contrainte en un point

et le vecteur cisaillement est gal :

2 - Contrainte sur une facette quelconque : tenseur des contraintes

Sous forme matricielle :

o [(M)], appel tenseur des contraintes en M, a pour expression :

Composantes du tenseur des contraintes

3 - Rciprocit des contraintes tangentielles

4 - Etat de contrainte en un point

5 - Contrainte normale dans une direction quelconque

RDM - Complments thoriques

III - Faces et contraintes principales

soit sous forme matricielle : [(M)] [ n ] = [n ] .

Les facettes X , Y , Z sont appeles faces principales.

Les quantits X , Y et Z sont les contraintes

Dans le repre principal { M , X Y Z }, le vecteur contrainte sur la facette n a pour expression :

o a, b et c sont les composantes de n .

IV - Cercles de Mohr des contraintes

= X cos sin + Y cos sin

Les trois directions principales nous permettent de construire trois

RDM - Complments thoriques

V - Etats de contraintes particuliers

2 - Etat de cisaillement simple

Dans le repre principal { M , X Y Z }, le tenseur des contraintes est gal :

3 - Etat de contrainte isotrope

L'tat de contraintes en un point M est isotrope si, quelque soit la facette n , on a T ( M, n ) = n .

4 - Etat de contrainte plan

L'axe z est donc direction principale et la contrainte

o u,v et w sont les composantes du vecteur dplacement du point M et la masse volumique du

yz ) est l'tat de contrainte en M.

RDM - Complments thoriques

B - Dformations autour dun point

Soient { O , x y z } un repre orthonorm et M(x,y,z) un point

On appelle vecteur dplacement du point M le vecteur MM ' .

Nous noterons ses composantes :

En dcomposant la matrice [D] en sa partie symtrique [E ] =

NN ' = MM ' + [] MN + [E] MN

points situs au voisinage de M subissent une rotation dintensit autour de laxe

II Etat de dformation en un point

Considrons deux points M et N du solide voisins l'un de

Dautre part : (MN 2 ) = (MN.MN) = MN.MN + MN.MN = 2 l 2 [n ] T [E] [n ]

quantits xx , yy et zz reprsentent donc respectivement lallongement unitaire en M dans les

RDM - Complments thoriques

2 - Glissement ou dformation angulaire

Soient deux points N1 et N2 voisins de M et tels que les

eux les deux vecteurs MN '1 et MN' 2 .

Soit langle que font entre eux les vecteurs MN1 et MN 2 .

(MN 1 . MN 2 ) = (l1l 2 cos ) = (l1 l 2 ) cos l1 l 2 sin = l1l 2 ( M, n 1 , n 2 ) car = .

Dautre part : (MN 1 . MN 2 ) = MN 1 . MN 2 + MN 1 . MN 2 = 2 l1l 2 [n 1 ]T [E ][n 2 ] .

3 - Etat de dformation en un point

III - Faces principales

Les quantits X , Y et Z sont appeles dformations principales.