II Tenseur des contraintes

II Tenseur des contraintes

Nous allons prsenter la notion de contrainte par extension de la mcanique des solides indformables. Cette
prsentation est directement lie la signification physique des contraintes, et offre l'avantage d'utiliser un
bagage que vous maitrisez.
A savoir : Le mouvement dun solide indformable est compltement dtermin ds que lon connat chaque
instant deux vecteurs : la force rsultante et le moment rsultant appliqus sur le solide (tat actuel).
En effet, la position du solide est dfinie par deux vecteurs galement (la position de son centre de
gravit et la rotation autour de ce centre de gravit) et le principe fondamental de la dynamique
permet de relier la force rsultante au dplacement du centre de gravit et le moment rsultant la
rotation autour du centre de gravit.
Nous allons isoler par la pense un petit lment de matire quelconque et crire le Principe Fondamental de
la Dynamique "PFD".

II-1 Tenseur des contraintes

Imaginons un corps soumis des efforts extrieurs, les contraintes sont des efforts intrieurs (ou efforts de
cohsion) qui se dveloppent au sein du milieu, ils sont dus aux interactions entre les particules du milieu.
Soit un corps en quilibre que nous coupons en deux parties. La 2


n
partie 1 reste l'quilibre sous les actions mcaniques extrieures

qui lui sont appliques et sous les actions de la partie 2 sur la dF


partie 1. = T ( P, n ) ds
ds
Nous admettrons que sur chaque lment de la surface de coupe
s'exerce une force lmentaire





1
dF = T ( P, n ) ds avec T ( P, n ) : pression ( N / m2 )




Le vecteur T ( P, n ) est appel vecteur contrainte en P par rapport la facette de normale n . Ce vecteur a la
dimension d'une pression.





Le principe de l'action raction
T ( P, n ) = T ( P, n )




Montrons que le vecteur contrainte T ( P, n ) dpend linairement de n

Ide : consiste isoler par la pense un petit lment de matire pour faire apparatre les efforts de cohsion
dans la matire.
Un lment de volume sera un infiniment petit d'ordre 3 alors que les lments de surface sont des infiniment
petits d'ordre 2, ce que nous notons :
dv
0
0( 3 ) et ds
0
0( 2 )
crivons les 3 quations, de rsultante, dduites du PFD pour un petit tridre de matire. La quantit
d'acclration est un infiniment petit d'ordre 3, ce qui revient crire les quations d'quilibre suivantes:.


e3

S : surface incline du tridre de normale n




Si : faces des plans (e j , ek ) du tridre, surfaces de normale ei





3
PFD
T ( P, n ) ds + T (Qi , ei ) dsi = 0( )
S1

i
S




T (Qi , ei ) = T (Qi , ei )
S2


e2








T ( P, n ) ds = T (Qi , ei ) dsi or dsi = ds [Link]

S3
e1 i

17
Mcanique des Milieux Continus & calcul des structures

T ( P, n ) = T (Qi , ei )ni Il existe donc un oprateur linaire n
T ( P, n ) = ( P ) n
i

Les 3 quations de moment conduisent la symtrie de l'oprateur : ij = ji

Dfinition :
Le tenseur des contraintes est dfini dans ltat dform, not C : tenseur des contraintes de Cauchy



Il permet d'exprimer la tension sur un lment de surface de normale n n




T ( P, n ) = C ( P ) n / m2 )
cest une pression. ( N





T ( P, n )
La force lmentaire exerce sur cette surface est : dF = T ( P, n ) ds P
T
ds
Ce tenseur est symtrique =

Dans un systme de coordonnes cartsienne C est reprsent

par :
11 12 13
= 22 23
C 12
13 23 33
ii contraintes normales
ij contraintes de cisaillement


Le matriau est dit en traction suivant la direction ei si la
contrainte normale ii > 0.
on parlera de compression dans le cas contraire.

La signification physique des coefficients de cette matrice est donne par e3
la figure ci-contre


T ( P, e1 )
11 11 31
21




T ( P, e1 ) = 21 = 12 11 e1

31 13
Les contraintes correspondent le plus souvent des pressions trs importantes (fonction de la cohsion du
matriau), c'est pourquoi le Mga-pascal est couramment utilis comme unit.
Pa = 1N / m 2 et MPa = 106 Pa = 1N / mm 2
Historiquement :
Vous rencontrerez, dans certains ouvrages anciens, l'unit traditionnelle "l'hectobar" qui reprsente 10 N par
mm 2 soit approximativement 1Kgf par mm 2 .
1hbar = 10 MPa .

Le tenseur des contraintes est une grandeur spatiale (dfinie dans l'tat actuel) pour pouvoir effectuer les
calculs dans une reprsentation Lagrangienne nous devons exprimer ce tenseur sur l'tat de rfrence.

T

Utilisons la transformation de surface prsente dans le chapitre prcdent : nds = J F Ndso
Nous obtenons L : Tenseur des contraintes de Piola-Lagrange ou 1er Tenseur de Piola-Kirchhoff


T
dF ( P ) = L Ndso avec L = J C F

18
 II Tenseur des contraintes

1

Puis exprimons dF ( P ) sur l'tat initial : dFo ( P ) = F dF ( P )

Nous obtenons K : Tenseur des contraintes de Piola-Kirchhoff

1 T

dFo ( P ) = L Ndso avec K = J F C F dFo pseudo tension
Ce tenseur est symtrique

En petites perturbations ces 3 tenseurs sont confondus, et on note : = C = K = L

ExerExercice1I-1 : quilibre de la matire

On considre une plaque d'paisseur uniforme h petite devant les autres dimensions. La plaque est soumise un


champ de force volumique f (pesanteur) et des actions mcaniques extrieures. On isole un petit paralllpipde
de matire dx1 dx2 h

1- Expliquer pourquoi ce paralllpipde est l'quilibre, faire une figure pour reprsenter tous les efforts
sur cet lment de matire.
2- Exprimer les efforts de cohsion en fonction du tenseur des contraintes.
3- En dduire les deux quations d'quilibre de rsultante et l'quation d'quilibre de moment

II-2 Proprits du tenseur des contraintes

La tension sur un lment de surface de normale n s'exprime par :
T ( P , n
)




n
T ( P , n
) = ( P ) n T = n
Ce vecteur peut se dcomposer en contrainte normale composante suivant la
normale de la tension sur la facette et en un vecteur tangentiel, dit contrainte ds
tangentielle ou contrainte de cisaillement.





= n
. n

Soit : T = n + avec




= T n = t

II-2.1 Thorme de Cauchy

Ce thorme souvent dit thorme de rciprocit des contraintes tangentielles, vient de la symtrie du
tenseur des contraintes (absence de couple de masse dans le milieu).

19
 Mcanique des Milieux Continus & calcul des structures

Thorme de Cauchy : n T ( P , n
)



pour toutes directions n , n '





Nous avons T ( P , n
) . n ' = T ( P , n
') . n

T ( P , n
')
La figure suivante donne une reprsentation du thorme dans le cas



particulier o n n '









Le thorme s'crit : n n n ' n . n ' = n ' . n n'

II-2.2 contraintes principales

Le tenseur des contraintes tant un tenseur symtrique d'ordre 2, il peut s'exprimer sur une base privilgie
dite base principale dans laquelle la matrice sera diagonale. Cette base propre orthonorme directe sera note




(eI , eII , eIII )
I 0 0
[ ]( e
I ,e
II ,e
III ) = 0 II 0
0 0 III
Les facettes, dont un des vecteurs de la base principale est la normale, ne subissent aucun cisaillement. Ce sont
des facettes en traction compression pure.

Si pour toute facette les contraintes de cisaillement sont nulles, c'est que l'tat de contrainte est
hydrostatique
0 0
[ ] = 0 0 avec = P si P est la pression hydrostatique
0 0
Un tat de contrainte est dit uni-axial (c'est une hypothse de modlisation) si deux valeurs propres sont
supposes ngligeables :
0 0
[ ] = 0 0 0 cas de l'hypothse des barres en traction = F /S
0 0 0
Un tat de contrainte est dit plan (c'est une hypothse de modlisation) si une valeur propre est suppose
ngligeable :
I 0 0 11 12 0
[ ] = 0 II 0
, ou dans une autre base du plan [ ] = 12 22 0



0 0 0 0 0 0


C'est l'tat de contrainte sur une surface, de normale eIII , non charge. "dans le vide"

Application :
Plaques et coques minces : Le milieu tant mince l'tat de contrainte, sur l'paisseur, s'cartera peu de
l'tat de contraintes sur les faces. Or si la pression atmosphrique est nglige devant les chargements
mcaniques extrieurs, l'tat de contrainte sur les faces d'une plaque ou d'une coque est plan. L'hypothse
consiste donc considrer que l'tat de contrainte est plan en tout point du milieu.

ExerExercice1I-2 : tat de contraintes particuliers

1- A quels tats de contraintes correspondent les situations suivantes ?
La contrainte normale toute facette est nulle.
La contrainte de cisaillement sur toute facette est nulle.

20
 II Tenseur des contraintes

10 5
= en MPa
5
2- Soit un tat de contrainte de surface dfini par :
5
Dterminer les contraintes principales et directions associes.
Dans quelle position une fissure dbouchant sur cette surface sera-t-elle la moins pnalisante ?

II-2.3 Tri-cercle de Mohr

Ide : on cherche une reprsentation de ltat contrainte ( , ) quand n varie.
Pour simplifier considrons dans un premier temps un tat de contraintes
planes, et plaons-nous dans la base principale.

t
I
I cos


cos n
On note n= et =
T =

sin II II sin eII
D'o la contrainte normale = I cos 2 + II sin 2


T ( P , n
)
et la contrainte de cisaillement = ( II I ) cos sin eI

I + II I II
= 2
+
2
cos 2
Soit :
= I II sin 2
2


le vecteur n dcrivant l'ensemble des directions du plan II 2 I
le point ( , ) dcrit un cercle.


+ II II T
de centre I de rayon I tat de contrainte
2 2 ( , )
Les contraintes principales sont ordonnes I > II
C'est le cercle de Mohr des contraintes planes

Pour un tat de contrainte tridimensionnel

Ce qui vient d'tre fait dans un plan principal de contrainte peut l'tre dans les deux autres plans ce qui nous
conduira au tri-cercle de Mohr.

1 n1 n


1 1
= 2 T = 2 n2
Notons :
et n = n2
sur la base principale.
3 n n
3 3 3





et la dcomposition T = n + on pose =


Nous avons : n =1
n1 + n2 + n3 = 1
2 2 2

T. n =
1n1 + 2 n2 + 3 n3 =
2 2 2

2
Et T = 2 + 2
12 n12 + 22 n22 + 32 n32 = 2 + 2

21
 Mcanique des Milieux Continus & calcul des structures

1 1 1 n12 1
n 2 =
Soit sous forme matricielle :
1 2 3 2
2 2 2
1 2 3 n3 +
2 2 2

( 2 + 3 2
) ( )
2
+ 2 2 2 3
n =2 2

1 ( 2 1 )( 3 1 )

( 3 + 1 2
) ( )
2
+ 2 3 21
n2 =
2 2
D'o l'on tire :
( 3 2 )( 1 2 )

( 1 + 2 2
) ( )
2
2 2 + 2 1 2 2
n1 = ( 1 3 )( 2 3 )

Pour une configuration 1 > 2 > 3 il faut vrifier ni

2
0 ce qui nous donne 3 conditions :

(
2 + 3
) + ( )
2 2 3 2
2
2 2

( ) + ( )
3 + 1 2 3 1 2
2
2
2

( ) + ( )
2 1 2 2
1 +2 2 2

Dans le plan ( , ) , la zone est donc dlimite par 3 cercles correspondants aux 3 cercles de Mohr des plans


principaux ni = 0 , nous limitons la reprsentation du tri-cercle de Mohr au demi-plan car 0 ( = )


n
MAX

n2 = 0

n3 = 0
n1 = 0
3 2 1

(n1 , n2 ) = 0 (n1 , n3 ) = 0 (n2 , n3 ) = 0

Les trois points sur l'axe horizontal correspondent aux 3 contraintes principales. La figure montre que le
cisaillement maximal est gal au rayon du plus grand cercle, et s'exerce sur une facette situe 45 dans le plan
dfini par les directions principales correspondantes aux contraintes principales extrmales.




Pour notre figure 1 > 2 > 3
MAX = ( 1 3 ) / 2 Sur une facette de normale n = 2
2
(n1 + n3 )
Sur cette facette la contrainte normale est = ( 1 + 3 ) / 2

Dans le cas gnral MAX = ( max( i ) min( i ) ) 2

ExerExercice1I-3 : tat de contrainte plane

On considre l'tat homogne de contrainte plane dfini par la figure ci dessous en ( MPa )

22
 II Tenseur des contraintes

30

20



n /3 e2 20
30


e1

1- Dterminer les contraintes principales et les directions associes

Par le calcul
En utilisant le cercle de Mohr


2- Calculer la contrainte normale et celle de cisaillement sur la facette de normale n (cf. figure).
Par le calcul
En utilisant le cercle de Mohr

II-3 Critres de limite lastique

Les critres servent dfinir le domaine des sollicitations pour lesquelles le comportement du matriau sera
lastique c'est--dire instantan et rversible. Ces critres s'expriment en fonction des invariants du tenseur
des contraintes c'est pourquoi nous en parlons dans ce chapitre.

Notion de critre
L'exprience montre que la notion de critre doit tenir compte de l'historique H 1du matriau. Le
critre de limite lastique est une fonction relle f ( , H ) telle que un tat de contrainte sera dit
admissible si De .
De = ;{ }
f ( , H ) 0

La surface f ( , H ) = 0 est dite surface Limite

Un tat de contrainte extrieur au domaine admissible est physiquement impossible atteindre
Pour un matriau isotrope f doit tre invariant dans tout changement de base, la reprsentation de surface
limite dans l'espace des contraintes principales est une surface pour laquelle, les plans (1,1,0), (0,1,1), (1,0,1)
sont des plans de symtrie, et laxe (1,1,1) est un axe de symtrie ternaire.

II-3.1 Invariants scalaires du tenseur des contraintes

Les coefficients du polynome caratristique P = det ( 1) sont les invariants du tenseur des contraintes.
1
Le premier invariant est la contrainte moyenne dfinie par m = tr
3
En contraintes principales : m = ( 1 + 2 + 3 ) / 3

1
L'historique ou paramtres de mmoire du matriau sont des notions relatives au cours sur les matriaux.
Sachez que pour les mtaux (matriaux ductiles, ils peuvent se dformer plastiquement) nous
parlerons d'crouissage. Par opposition un matriau est dit fragile s'il ne se dforme pas plastiquement (verre).

23
Mcanique des Milieux Continus & calcul des structures

Le second et troisime invariants s'expriment en fonction du dviateur des contraintes : s = m 1 . Le

terme m1 est appel partie sphrique du tenseur des contraintes, et s la partie dviatorique du tenseur.
En contraintes principales : si = i m
2 3
J 2 = tr ( s ) et J 3 = tr ( s ) sont les second et troisime invariants du tenseur.
En contraintes principales : J 2 = s12 + s22 + s32
J 3 = s13 + s23 + s33
On peut noter que : tr s = 0
J 3 = 3det s
Un peu de physique :
L'exprience montre qu'un lment de matire homogne isotrope (mtaux) conserve un comportement
lastique quelque soit l'tat de contrainte hydrostatique qui lui est appliqu. or un tat de contrainte
hydrostatique, correspond un dviateur des contraintes nul, c'est pourquoi les critres d'lasticit des
mtaux sont indpendant du premier invariant, ils ne dpendent que du dviateur des contraintes.
Notez qu'un sable ou un sol ne peut pas tre assimil un solide homogne isotrope (grains), c'est pourquoi
pour ces matriaux les critres seront fonction du premier invariant "pression moyenne".

II-3.2 critres pour les mtaux

Ces critres sont anciens, ils restent cependant largement utiliss car les mthodes didentification fiables des
paramtres du modle caractrisant le comportement des matriaux ont fait leurs preuves pour ces critres.
Critre de Tresca (1870)
Ce critre, le plus ancien, est d'origine physique il consiste observer que l'tat limite est atteint pour une
valeur de cisaillement maximale qui est une caractristique du matriau.

La surface limite est dfinie par : Max o 0

1
En contraintes principales : Max = max i j
2 ij

Ce critre est bien indpendant de la partie sphrique du tenseur des contraintes. Le fait ajouter un tat de
contrainte hydrostatique ne modifie pas la valeur du seuil.

Pour un essai de traction, le seuil est dfini par o = o / 2 .

Si l'on utilise cet essai comme rfrence, la contrainte de Tresca est : T = max i j
ij

La surface limite est dfinie par T = o sa reprsentation graphique dans l'espace des contraintes principales
est une surface cylindrique de base hexagonale non borne daxe (1,1,1).

3
1,1,1 3
Von Mises

fo() Tresca

2
1 2
1
Critre de Von Mises (1910)
Ayant observ que le critre de Tresca tait trop svre, Von Mises proposa un critre plus mathmatique
bas sur le second invariant du tenseur des contraintes.

La surface limite est dfinie par : J2 s 0

24
 II Tenseur des contraintes

s a la dimension d'une contrainte et reprsente la valeur critique (seuil)

En pratique on souhaite identifier le seuil la valeur de o de lessai de traction.

Exprimons J 2 en fonction des contraintes principales

3
J 2 = ( i m ) 2 avec m = ( 1 + 2 + 3 ) / 3
i =1

J2 =
1
3
(
( 1 2 ) + ( 2 3 ) + ( 3 1 )
2 2 2
)
2 2
Pour l'essai de traction 2 = 3 = 0
J 2 = 1
3
3
La grandeur VM = J 2 est appele contrainte quivalente de Von Mises.
2
1
(1 2 )2 +( 2 3 ) +( 3 1 )
2 2
En contraintes principales : VM =
2
La surface limite est dfinie par VM = o sa reprsentation graphique dans l'espace des contraintes
principales est une surface cylindrique non borne daxe (1,1,1) de rayon 2 o
1,1,1 3
3 R = 2 o

fo()
1 2
2 Vue (1,1,1) dans le plan dquation :
1 + 2 + 3 = 0
1
Lcart maximal entre le critre de Von Mises et celui de Tresca est 1 3 2 soit un cart de lordre de 13%.
De ce fait, le critre de Tresca plus restrictif est moins utilis pour le dimensionnement des pices mcaniques,
et la plus part des codes de calcul (lments finis) en mcanique proposent, par dfaut, la visualisation de la
contrainte quivalente de Von Mises.

II-3.3 critres pour les matriaux granulaires (sols)

Ces critres doivent traduire le fait que la rsistance au cisaillement est dautant plus forte que la compression
(confinement) du milieu est leve.
Max ( m ) est fonction du premier invariant, et Max augmente avec m il faut penser m < 0
Critre de Mohr Coulomb : C m tan
C : Constante de cohsion, cest la rsistance au cisaillement en absence de compression
C = 0 le matriau est dit pulvrulent (sable)
: Angle de frottement, il traduit laugmentation de la rsistance au cisaillement avec le confinement
(frottement de contact inter-granulaire)
= 0 le matriau est parfaitement cohrant
on retrouve le critre de Tresca
La surface limite est dfinie par
(1 + sin ) Max( i ) (1 sin ) Min( i ) = 2C cos
sa reprsentation graphique dans l'espace des contraintes principales est une surface conique non borne daxe
(1,1,1).

25
 Mcanique des Milieux Continus & calcul des structures

1,1,1
3 3
C cotg

2
1 2
1

Critre de Drcker Prager : J 2 CD 3 m tan D

La signification physique de ce critre est identique celle du critre de Mohr Coulomb avec comme support
mathmatique le critre de Von-Mises

ExerExercice1I-4 : Essai de traction

On tudie ici lessai classique de traction. Lprouvette est cylindrique de diamtre d=10 mm et de hauteur h=10
cm. Elle est soumise un effort F=200 N

1- Dans la zone centrale de l'prouvette l'tat de contrainte est suppos uniforme et uniaxial, donner son
expression.
2- Exprimer, en fonction de F et S, la contrainte normale et celle de cisaillement sur une facette incline
d'un angle quelconque.
Caractriser les tats de contrainte pour = 0 et 90.
Quelle est la valeur du cisaillement maximal sur quels plans agit-il ?
Retrouver ces rsultats en utilisant le cercle de Mohr
3- Reprsenter le critre de Tresca dans le plan ( , ) .
En dduire la valeur du cisaillement limite en fonction de la limite lastique du matriau o .
4- Reprsenter le critre de Mohr-Coulomb " C " dans le plan ( , ) .
Dans quel plan se produira la rupture de l'prouvette ?
Qu'elle est la charge limite de cette prouvette ?
Analyser les rsultats.

II-4 quations du mouvement

II-4.1 lments de calcul indispensables
Nous rappelons ici sans en donner la dmonstration les thormes mathmatiques qui vont tre utiliss dans
ce chapitre.
Formules de Green Ostrogradsky
La formule de Green Ostrogradsky ou thorme de la divergence permet de passer d'une intgrale de
volume une intgrale de surface, c'est une formule d'intgration par parties.



Pour a scalaire grada dV =
D
an dS
D





Pour a vecteur diva dV =
D
a. n dS
D


Pour a tenseur diva dV = a n dS
D D

Drivation d'une intgrale : Thorme de Reynolds

26
 II Tenseur des contraintes

d dA


A dV = ( + A divv ) dV Formule valable avec a scalaire ou a vecteur
dt D D dt
Application : crivons le TH de Reynolds dans le cas d'une grandeur scalaire.
d

g dV = ( g + g divv ) dV g 


dt D Utilisons La drive Eulrienne g =+ gradg .v
D
t





Nous obtenons : Et gradg .v + gdivv = div( gv )
d g


dt D
g dV = + div( gv ) dV
D
t
d g


Appliquons le Th d'Ostrogradsky:
dt D
g dV = dV +
D
t g v . n dS
D

II-4.2 conservation de la masse

C'est une application directe des formules prcdentes
dM

La masse du systme est dfinie par M = dV
= ( + divv ) dV
D
dt D
dM

crire = 0 pour tout domaine
+ divv = 0
dt



g

Avec la drive Eulrienne : + grad .v + divv = 0 soit + div( v ) = 0
t t
C'est l'quation d'Euler
En utilisant la transformation gomtrique la conservation de la masse peut s'crire :
1
dm = o dVo = dV or dV = JdVo
=
o J

II-4.3 quations du mouvement

PFD
crivons le Principe Fondamental de la Dynamique pour un domaine matriel quelconque nous avons :




Les 3 quations de rsultantes sont : a dV = f dV + T dS
D D D
Ces quations sont crites dans l'tat dform






Or T = n
a dV = f dV + n dS
D D D




Le Th d'Ostrogradsky
a dV = f dV + div ( ) dV
D D D

Soit : (
D




)
a f div( ) dV = 0 et ceci D



Nous obtenons les 3 quations locales : a = f + div( )
Ces quations sont dfinies l'intrieur du domaine, dans son tat actuel.

Les 3 quations de moment donnent : ij = ji Symtrie du tenseur des contraintes

27
Mcanique des Milieux Continus & calcul des structures

Sur la frontire du domaine nous devons crire des conditions aux limites (changes avec l'extrieur), ces
conditions aux limites sont de deux types :



Dplacements imposs u = ud sur la frontire Du condition dite de Dirichlet

Forces imposes n = Td sur la frontire D condition dite de Neumann
Bien entendu D = Du D et Du D =
Il est impossible d'imposer la fois le dplacement et les efforts
Du point de vue pratique l'ingnieur parlera de conditions aux limites cinmatiques (dplacements ou vitesses)
et naturelles (forces ou pressions). Cette classification sera essentielle pour lorsque nous aborderons la
rsolution de ces problmes.

La mise en quations par le PFD conduit un systme d'quations aux Drives Partielles "EDP"
Pour le PFD : Nous obtenons le problme aux limites reprsent par la figure ci-dessous




u = f + div
M Du



u = ud M D
M D



n = Td



n Td

PTV


crivons maintenant le Principe des Travaux Virtuels : u A = Te + Ti tout instant



A = a. u dV Reprsente le travail virtuel des quantits d'acclration
D





Te = f . u dV + T . u dS Reprsente le travail virtuel des efforts extrieurs
D D

Ti = : dV Reprsente le travail virtuel des efforts intrieurs, les contraintes

D

ii travaillent dans des dformations ii , et les cisaillements ij dans des "glissements" ij .

Le produit doublement contract des deux tenseurs signifie que l'on multiplie terme terme les
coefficients des matrices correspondantes. Compte tenu de la symtrie des tenseurs, l'nergie de
dformation lmentaire est :
2 Ed = : = 1111 + 22 22 + 33 33 + 2 1212 + 2 1313 + 2 23 23
On peut utiliser les glissements ij = 2 ij dans cette expression

= grad ( u ) + grad ( u )
1 T


Le tenseur taux de dformation est dfini par
2
Pour le PTV : Nous obtenons la formulation variationnelle du problme








u u. u dV = : dV + f . u dV + T . u dS
D D D D
tout instant, pour tout systme
L'intrt de cette formulation est qu'elle se prte bien aux calculs numriques. Nous prsentons dans un autre
cours "MEF : Mthode des lments Finis " comment passer de cette forme intgrale une quation
matricielle calcule partir d'une discrtisation du domaine en lments de formes gomtriques simples

28
 II Tenseur des contraintes

Les deux principes de la mcanique noncs prcdemment sont bien entendu quivalents





Partons de lquation locale a f div( ) = 0 sur le domaine






"EL" u u.( u div( ) - f ) dV = 0
D





Or div( u ) = : grad ( u ) + div . u
En cartsiennes ( pi ij ), j = ij pi , j ij , j pi
De plus est un tenseur symtrique donc

T



div( u ) = : grad ( u ) + div . u

grad T u
+ grad u




Nous avons donc div . u = div ( u ) + :
2

En reportant dans la premire relation nous obtenons :








"EL" u ( u . u + : grad s ( u ) div ( u ) - u . f ) dV = 0
D Utilise la symtrie du tenseur
des contraintes
Appliquons le TH d'Ostrogradsky






div( u ) dV = u. n dS = u. n dS
D D D









"EL" u ( u . u + : grad s ( u ) - u . f ) dV ( u. n ) dS = 0
D D

En tenant compte des conditions aux limites sur la frontire D M D n = Td
"EL."

[Link] sur











u ( u . u + : grad s ( u ) - u . f ) dV u . n dS u .Td . dS = 0
D Du D

Soit le PTV : n = TI Champ des efforts








u u. u dV + : dV = f . u dV + T . u dS inconnus sur cette frontire

D D D D

Si nous utilisons un champ de dplacement virtuel valeur nulle sur la frontire Du

Champ de dplacement virtuel
M Du u( M ) =0
u. n dS = 0 cinmatiquement admissible
Du

Nous obtenons la formulation variationnelle du problme

M Du u = ud soit uCA








uCA u. u dV + : dV = f . u dV + T D . u dS
D D D D

29
Mcanique des Milieux Continus & calcul des structures

Notes personnelles

30