1 / 95
Rappels
Dfinitions
MEC6418 - NOTES DE COURS
Oprations
Symtries
Symtries mat.
Projecteurs.
Notions lmentaires sur les tenseurs
Const. Ingnieur
Par: Martin Lvesque
professeur du dpartement de gnie mcanique
Hiver 2011
Espaces vectoriels
Dfinitions
Corps
Rappels
Espace
App. n-Lin.
Ensemble de nombres (R, C, etc.) utiliss pour dfinir un
espace vectoriel.
Espace vectoriel linairel
Oprations
Symtries
2 / 95
Symtries mat.
Projecteurs.
Ensemble dentits mathmatiques (scalaires, vecteurs,
fonctions, matrices, etc.) dfinis sur un corps et ayant
certaines proprits communes.
Const. Ingnieur
Par exemple, R3 , qui pourrait tre not R R R,
pourrait tre lensemble des triplets (x1 , x2 , x3 ) o
xi R.
Les rgles pour un espace vectoriel linaire sont:
1. Si x et y E, alors x + y E.
2. Si R et x E, alors x E
Il est intressant de noter que les objets dun espace
vectoriel (x, y, etc.) sont appels vecteurs mme sil sagit
de scalaires, vecteurs, matrices, fonctions, etc.
Espaces vectoriels
Rappels
Espace
App. n-Lin.
Dfinitions
Oprations
Symtries
3 / 95
Produit scalaire
Pour x, y E, le produit scalaire dnot < x, y > est une
fonction qui associe le doublet (x, y), compris dans lespace
E E, un lment de R. En langage mathmatique,
cette phrase scrit: < , >: E E R .
Les proprits du produit scalaire sont les suivantes:
Symtries mat.
1. < x, y >= < x, y >
2. < x + y, z >=< x, z > + < y, z >
3. < x, x >> 0 si x 6= 0, < x, x >= 0 x = 0
Projecteurs.
Const. Ingnieur
Des exemples de produit scalaires sont:
Pn
1. < x, y >= k=1 xk yk o x, y Rn
R1
2. < f, g >= 0 f (x)g(x)dx o f, g sont des fonctions
continues pour x [0, 1]
Application n-linaire
Rappels
Espace
App. n-Lin.
4 / 95
Une application n-linaire est une application qui fait
correspondre un vecteur w, n vecteurs. La Figure 1 illustre
schmatiquement une application bi-linaire.
Dfinitions
Oprations
Symtries
Symtries mat.
Projecteurs.
Const. Ingnieur
Application
n-linaire
Figure 1: Reprsentation dune application bi-linaire qui prend le
couple de vecteurs (u, v) de lespace E F et y fait correspondre
un vecteur w dans un espace S. Si on nomme l cette application
linaire on crira: l : E F S.
Application n-linaire
Rappels
Espace
App. n-Lin.
Dfinitions
Oprations
Symtries
Symtries mat.
Projecteurs.
5 / 95
Les proprits de n-linarit sont les suivantes. Soient
(u1 , u2 , . . . , un ) n vecteurs appartenant diffrents espaces
vectoriels, 1 , 2 , . . . n R et T cette application n-linaire.
On aura:
T (1 u1 , 2 u2 , . . . , n un ) =
(1 2 . . . n ) T (u1 , u2 , . . . , un ) (1a)
Const. Ingnieur
T (u1 + un+1 , u2 + un+2 , . . . , un + un+n ) =
iX
1 =1 iX
2 =1
i1 =0 i2 =0
...
iX
n =1
in =0
T (uni1 +1 , uni2 +2 , . . . , unin +n ) (1b)
Application n-linaire et forme n-linaire
6 / 95
Par exemple, pour une application bi-linaire, on aurait
Rappels
Espace
App. n-Lin.
T (1 u1 , 2 u2 ) = 1 2 T (u1 , u2 )
(2a)
Dfinitions
T (u1 + u3 , u2 + u4 ) =
Oprations
Symtries
T (u1 , u2 ) + T (u1 , u4 ) + T (u3 , u2 ) + T (u3 , u4 ) (2b)
Symtries mat.
Projecteurs.
Const. Ingnieur
Une forme n-linaire est une application n-linaire qui fait
correspondre n vecteurs un scalaire appartenant un corps.
Particularits notre cours
On notera par E lensemble de tous les triplets de scalaires
rels permettant de dfinir tous les vecteurs de dimension 3.
En dautres termes E est R R R, o R est considr
comme un espace vectoriel.
On notera par {~
e1 , ~e2 , ~e3 } les vecteurs formant une base
orthonorme de E.
Tenseur dordre 1
Dfinitions
Tens. 1
Tens. 2
Tens. 4
Oprations
Un tenseur dordre 1 permet de dfinir une forme linaire
sur E.
Si on note par t un tenseur dordre 1, alors on aura:
t(u) = o u E et R
Symtries
Symtries mat.
Projecteurs.
Dfinition:
Rappels
7 / 95
Const. Ingnieur
(3)
Le produit scalaire classique est une forme linaire. Si u E,
dans ce cas on aurait:
t(u) =
i=3
X
t i ui
i=1
= t 1 u1 + t 2 u2 + t 3 u3
(4)
On peut montrer que les composantes de t sont donnes par:
t (~ei ) = ti
Voir dmonstration au tableau pour la linarit et (5).
(5)
Tenseur dordre 2
8 / 95
Un tenseur dordre 2 permet de dfinir une forme bi-linaire sur
E E.
Si on note par un tenseur dordre 2, alors on aura:
Rappels
Dfinitions
Tens. 1
Tens. 2
Tens. 4
(u, v) = o u, v E et R
Oprations
Symtries
La forme bi-linaire suivante est un tenseur dordre 2:
Symtries mat.
(u, v) =
Projecteurs.
(6)
Const. Ingnieur
j=3
i=3 X
X
ui ij vj
i=1 j=1
i=3
X
ui i1 v1 + ui i2 v2 + ui i3 v3
i=1
= u1 11 v1 + u2 21 v1 + u3 31 v1 +
u1 12 v2 + u2 22 v2 + u3 32 v2 +
u1 13 v3 + u2 23 v3 + u3 33 v3
=
Voir dmonstration au tableau pour la bi-linarit
(7)
Tenseur dordre 4
9 / 95
Un tenseur dordre 4 peut avoir plusieurs significations.
Forme 4-linaire de E 4 : T(u, v, w, a) = o u, v, w et a E
Rappels
Dfinitions
Tens. 1
Tens. 2
Tens. 4
et R.
Par exemple,
Oprations
T(u, v, w, a) =
Symtries
Symtries mat.
ui vj wk al Tijkl =
(8)
i=1 j=1 k=1 l=1
Projecteurs.
Const. Ingnieur
3 X
3 X
3 X
3
X
La proprit de 4-linarit se vrifie de manire analogue
la bi-linarit pour les tenseurs dordre 2.
Cette interprtation est rarement utilise en mcanique.
Tenseur dordre 4
10 / 95
Un tenseur dordre 4 peut avoir plusieurs significations.
Forme bi-linaire de F F , o F est lespace des tenseurs
dordre 2.
Rappels
Dfinitions
Tens. 1
Tens. 2
Tens. 4
Oprations
T , =
Symtries
Symtries mat.
Projecteurs.
Const. Ingnieur
Par exemple,
3
3 X
3 X
3 X
X
ij Tjikl kl =
i=1 j=1 k=1 l=1
Si et sont la dformation et si T est le tenseur de
rigidit C, alors un tenseur dordre 4 peut servir dfinir
lnergie de dformation.
(9)
Tenseur dordre 4
11 / 95
Un tenseur dordre 4 peut avoir plusieurs significations.
Application linaire de F F , o F est lespace des tenseurs
dordre 2.
Rappels
Dfinitions
Tens. 1
Tens. 2
Tens. 4
3
3 X
X
T =
Tijkl kl = ij =
Oprations
Symtries
(10)
k=1 l=1
Symtries mat.
Projecteurs.
Const. Ingnieur
Par exemple,
Si est la dformation et si T est le tenseur de rigidit C,
est le tenseur des contraintes .
Un tenseur dordre 4 peut donc faire le lien entre les
dformations et les contraintes.
Produit tensoriel
Rappels
Dfinitions
Oprations
Produit
Einstein
Contraction
Delta
Identits
Drive
Ch. base
Symtries mat.
Projecteurs.
Const. Ingnieur
Deux tenseurs dordre m et n permettent de crer un tenseur
dordre m + n par lopration suivante:
(11)
Ai1 i2 ...im Bj1 j2 ...jn = Ai1 i2 ...im Bj1 j2 ...jn
Symtries
12 / 95
Par exemple:
a1 b1 a1 b2 a1 b3
a b = = ai bj = ij = a2 b1 a2 b2 a2 b3
a3 b1 a3 b2 a3 b3
o a t reprsent comme une matrice.
(12)
Produit contract et notation dEinstein
13 / 95
Le produit
P3scalaire de deux vecteurs peut sexprimer par
~u ~v = i=1 ui vi .
Avec la notation dEinstein, on suppose quil y a somme sur
tous les indices rpts. On pourra donc crire:
Rappels
Dfinitions
Oprations
Produit
Einstein
Contraction
Delta
Identits
Drive
Ch. base
Symtries
Symtries mat.
Projecteurs.
Const. Ingnieur
3
X
ui vi = ui vi
(13)
i=1
On qualifiera les indices rpts de muets car on aurait pu
obtenir le mme rsultat en crivant uj vj , uz vz , etc.
On qualifiera de franc un indice qui ne disparat pas lors de
lopration.
Produit contract et notation dEinstein
Le produit n-fois contract consiste dans un premier temps
effectuer le produit tensoriel entre deux tenseurs et par la suite
rpter n paires dindices o un indice appartient chaque
tenseur originel.
Par exemple, considrons le produit simplement contract, not
par un point (), entre deux tenseurs dordre 1, qui scrira:
u v:
Rappels
Dfinitions
Oprations
Produit
Einstein
Contraction
Delta
Identits
Drive
Ch. base
Symtries
Symtries mat.
Projecteurs.
Const. Ingnieur
14 / 95
1. On ralise le produit tensoriel: u v = ui vj
2. On rpte une paire dindices: ui vj devient ui vi
Alors, en vertu de la convention dEinstein,
u v = u1 v1 + u2 v2 + u3 v3 .
Produit contract et notation dEinstein
Rappels
Dfinitions
Oprations
Produit
Einstein
Contraction
Delta
Identits
Drive
Ch. base
Symtries
Symtries mat.
Projecteurs.
Const. Ingnieur
15 / 95
Considrons le produit simplement contract entre un tenseur
dordre 2 et un tenseur dordre 1, qui scrira: u:
1. On ralise le produit tensoriel: u = ij uk
2. On rpte une paire dindices: ij uk devient ij uj = vi
On aura un indice franc (i) et un indice rpt (j), ce qui
donne comme rsultat un tenseur dordre 1.
Pour le produit doublement contract, not par (:), on
contracte deux paires dindices. Par exemple, si on calculait
: , on aurait:
1. On ralise le produit tensoriel: = ij kl
2. On rpte deux paires dindices: ij kl devient ij ji = a
ce qui est un scalaire.
Produit contract et notation dEinstein
Rappels
Dfinitions
Oprations
Produit
Einstein
Contraction
Delta
Identits
Drive
Ch. base
Symtries
Symtries mat.
Projecteurs.
Const. Ingnieur
16 / 95
Le produit C : conduit
1. On ralise le produit tensoriel: C = Cijkl mn
2. On rpte deux paires dindices: Cijkl mn devient
Cijkl kl = ij
On aura deux indices francs (ij) et deux indices rpts (kl), ce
qui donne comme rsultat un tenseur dordre 2.
Le produit C : M conduit
1. On ralise le produit tensoriel: C M = Cijkl Mmnop
2. On rpte deux paires dindices: Cijkl Mmnop devient
Cijkl Mklop = Aijop
On aura quatre indices francs (ijop) et deux indices rpts
(kl), ce qui donne comme rsultat un tenseur dordre 4.
Produit contract et notation dEinstein
Rappels
Dfinitions
Oprations
Produit
Einstein
Contraction
Delta
Identits
Drive
Ch. base
Symtries
Symtries mat.
Projecteurs.
Const. Ingnieur
17 / 95
Pour le produit quatre fois contract, not par (::), on
contracte quatre paires dindices. Par exemple, si on calculait
A :: B, on aurait:
1. On ralise le produit tensoriel: A B = Aijkl Bmnop
2. On rpte quatre paires dindices: Aijkl Bmnop devient
Aijkl Bklij
ce qui est un scalaire.
Il doit tre not que la position des indices contracts vient
dune convention. On aurait pu contracter des indices des
positions diffrentes. On verra pourquoi plus tard cette
convention est utile.
Il est intressant de noter que les produits contracts
permettent de dfinir des formes et des applications linaires.
Le Delta de Kronecker
Rappels
Dfinitions
Oprations
Produit
Einstein
Contraction
Delta
Identits
Drive
Ch. base
18 / 95
Par dfinition, le Delta de Kronecker, ij est:
(
1 pour i = j
ij =
0 autrement
(14)
Cest une quantit fort utile pour les manipulations des tenseurs
et il est important de la matriser.
Par exemple, le produit:
Symtries
: = ij jk = ik =
Symtries mat.
Projecteurs.
(15)
Voir dmonstration complte au tableau.
Const. Ingnieur
La contraction dun tenseur dordre 2 avec le fait en sorte que
le tenseur rsultant est celui qui a t contract o lon
remplace lindice contract par celui qui ne lest pas dans le
Delta de Kronecker.
Les tenseurs identits
19 / 95
Lidentit dordre deux, nomme i, est dfinie par i : = . On
voit avec lquation (15) que ikl = kl .
Lidentit dordre 4 gnrale, nomme Ig , fait en sorte que
Rappels
Dfinitions
Oprations
Produit
Einstein
Contraction
Delta
Identits
Drive
Ch. base
Ig : A = A
g
= ik jl et que lon calcule Ig : A, on
Si lon dfinit Iijkl
aura:
1. Produit tensoriel Ig A = ik jl Amnop
2. Contraction sur deux paire dindices:
ik jl Amnop ik jl Aklop .
3. Application de la proprit du Delta de Kronecker:
ik jl Aklop = Aijop = A.
Symtries
Symtries mat.
Projecteurs.
Const. Ingnieur
On a donc dmontr que la dfinition didentit propose
rencontre la proprit de lidentit.
Les tenseurs identits
Rappels
Dfinitions
Oprations
Produit
Einstein
Contraction
Delta
Identits
Drive
Ch. base
20 / 95
Lidentit gnrale nest pas utilise en mcanique car elle ne
conduit pas une reprsentation matricielle (on verra plus loin).
Comme les tenseurs dordre 4 utiliss en mcanique sont
symtriques, on utilisera lidentit dfinie par:
Iijkl
ik jl + il jk
=
2
(16)
Si on calcule I : A, on trouve que:
Symtries
Symtries mat.
Projecteurs.
Const. Ingnieur
Aijkl + Ajikl
I:A=
2
ce qui est gal A uniquement si Aijkl = Ajikl
(17)
Drive par rapport un tenseur
Rappels
Dfinitions
Symtries mat.
Projecteurs.
Const. Ingnieur
(18)
o, comme pour le produit tensoriel, le tenseur rsultant est un
tenseur dordre m + n.
Symtries
La drive dun tenseur A E m par rapport un tenseur
B E n est dfinie par:
Ai1 i2 ...im
= Ci1 i2 ...im j1 j2 ...jn
Bj1 j2 ...jn
Oprations
Produit
Einstein
Contraction
Delta
Identits
Drive
Ch. base
21 / 95
Voir exemple au tableau portant sur la drive de
lnergie de dformation.
Changement de base dun tenseur
22 / 95
Problme
Rappels
Dfinitions
Oprations
Produit
Einstein
Contraction
Delta
Identits
Drive
Ch. base
On souhaiterait, partir de la connaissance des composantes de
~ exprimes dans le repre xyz obtenir celles dans le repre
V
x y z
Cela peut se faire laide de la connaissance des cosinus
directeurs. En 2D on aura les angles directeurs:
Symtries
Symtries mat.
Projecteurs.
Const. Ingnieur
Figure 2: Illustration des angles directeurs
Changement de base dun tenseur
23 / 95
Rappels
Dfinitions
Oprations
Produit
Einstein
Contraction
Delta
Identits
Drive
Ch. base
Symtries
Symtries mat.
Projecteurs.
En 3D, on aura:
Const. Ingnieur
T
Vx
Vy
=
Vz
Vx Vy Vz
cos (x , x)
cos (y , x)
cos (z , x)
cos (x , y) cos (y , y) cos (z , y) (19)
cos (x , z) cos (y , z) cos (z , z)
Changement de base dun tenseur
24 / 95
Rappels
Dfinitions
Oprations
Produit
Einstein
Contraction
Delta
Identits
Drive
Ch. base
Symtries
Symtries mat.
Projecteurs.
Const. Ingnieur
8 9T
< Vx =
Vy
= Vx
: ;
Vz
Vy
Vz
cos (x , x)
4 cos (x , y)
cos (x , z)
cos (y , x)
cos (y , y)
cos (y , z)
cos (z , x)
cos (z , y) 5
cos (z , z)
On peut voir que ces oprations consistent projeter les
~ sur les vecteurs de la base x y z .
composantes de V
On peut aussi voir que les colonnes de la matrice sont les
composantes de x , y et z , respectivement, exprimes dans la
base {x, y, z}
Cette matrice, nomme matrice des cosinus directeurs, est
note par [P ] et permet de dfinir les relations suivantes:
Vi = Pji Vj
(20a)
Vi = Pij Vj
(20b)
Changement de base dun tenseur
25 / 95
Vi = Pji Vj
Rappels
Dfinitions
Vi = Pij Vj
Oprations
Produit
Einstein
Contraction
Delta
Identits
Drive
Ch. base
lquation prcdente on a utilis le fait que Pij Pjk = ik (ou
autrement dit [P ]1 = [P ]T ).
Cela est facilement dmontrable car la matrice [P ] est
orthogonale, cest dire quelle est forme de vecteurs colonne
qui sont orthogonaux et unitaires.
Symtries
Symtries mat.
Projecteurs.
Const. Ingnieur
Le produit Pij Pjk est en fait le produit scalaire entre le
vecteur i et le vecteur k. Comme les vecteurs sont unitaires,
alors le produit du vecteur i avec lui-mme donne 1.
Comme les vecteurs sont orthogonaux, le produit entre un
vecteur i et un vecteur k, pour i 6= k, donne 0. On retrouve
la dfinition du Delta de Kronecker, qui est aussi lidentit
dordre 2.
Changement de base dun tenseur
Rappels
Pour un tenseur dordre 2, on avait que:
(u, v) = ui ij vj = o u, v E et R
Dfinitions
Oprations
Produit
Einstein
Contraction
Delta
Identits
Drive
Ch. base
26 / 95
(22)
Imaginons maintenant que u et v sont exprims dans une autre
base {x , y , z }. On aura ui = Pik uk et vi = Pik vk .
Lquation (22) devient alors:
(u, v) = ui ij vj
Symtries
Symtries mat.
Projecteurs.
Const. Ingnieur
ij Pjl vl
(Pik Pjl ij ) uk vl
uk vl
kl
Pik uk
(23)
= u ,v
=
Puisque les vecteurs u, v et u , v reprsentent la mme
quantit, il est normal que (u, v) = (u , v ). On tire donc la
formule de changement de base:
= Pik Pjl ij
kl
(24)
Changement de base dun tenseur
Rappels
Dfinitions
Oprations
Produit
Einstein
Contraction
Delta
Identits
Drive
Ch. base
Pour un tenseur dordre 4, on procde dune manire
compltement analogue. On peut utiliser nimporte quelle
interprtation dj donne. On aura par exemple, pour u, v, w
et a E:
T(u, v, w, a) = Tijkl ui vj wk al
= Tijkl Pim um Pjn vn Pko wo Plp ap
= (Pim Pjn Pko Plp Tijkl ) um vn wo ap
Symtries
Symtries mat.
27 / 95
Tmnop
um vn wo ap
(25)
= T u ,v ,w ,a
Projecteurs.
Const. Ingnieur
On tire donc la formule de changement de base pour un tenseur
dordre 4:
Tmnop
= Pim Pjn Pko Plp Tijkl
(26)
Symtries majeures et mineures
28 / 95
Tenseur dordre 2
Rappels
Dfinitions
Oprations
Symtries
Symtries
Reprsentation
Conventions
Un tenseur dordre 2 est dit symtrique si ij = ji . Cela fait
passer le nombre de constantes indpendantes de 9 6.
La majorit des tenseurs dordre 2 utiliss en mcanique sont
symtriques (contraintes, dformations).
Tenseur dordre 4
Symtries mat.
Projecteurs.
Const. Ingnieur
Un tenseur dordre 4 qui a les symtries dites mineures est tel
que Aijkl = Ajikl = Aijlk = Ajilk .
Un tenseur dordre 4 qui a les symtries dites majeures est tel
que Aijkl = Aklij .
Un tenseur dordre 4 qui a les symtries majeures et mineures a
jusqu 21 constantes indpendantes (au lieu de 81 quand il ny
a pas de symtrie.
La grande majorit des tenseurs dordre 4 en mcanique du
solide possdent les symtries majeures et mineures.
Reprsentation matricielle des tenseurs symtriques
29 / 95
Lorsque les tenseurs dordre 2 et dordre 4 prsentent les
symtries numres au transparent prcdent, il est possible de
reprsenter un tenseur dordre 2 comme un vecteur 6
composantes et un tenseur dordre 4 par une matrice 6 6.
Par exemple, un tenseur dordre 2 peut tre reprsent par:
Rappels
Dfinitions
Oprations
Symtries
Symtries
Reprsentation
Conventions
Symtries mat.
Projecteurs.
Const. Ingnieur
11
22
33
223
231
212
= I
(27)
Dans toute la suite, les indices majuscules (e.g. I, J, etc) seront
utiliss pour indiquer que lont fait rfrence un tenseur
exprim sous la notation matricielle.
Reprsentation matricielle des tenseurs symtriques
Rappels
30 / 95
En notation matricielle, un tenseur dordre 4 pourra scrire
sous la forme:
Dfinitions
Oprations
Symtries
Symtries
Reprsentation
Conventions
Symtries mat.
Projecteurs.
Const. Ingnieur
A=
A1111
A2211
A3311
2A2311
2A3111
2A1211
A1122
A2222
A3322
2A2322
2A3122
2A1222
A1133
A2233
A3333
2A2333
2A3133
2A1233
2A1123
2A2223
2A3323
2A2323
2A3123
2A1223
2A1131
2A2231
2A3331
2A2331
2A3131
2A1231
2A1112
2A2212
2A3312
2A2312
2A3112
2A1212
(28)
On crira alors que A = AIJ
Reprsentation matricielle des tenseurs symtriques
31 / 95
Cette notation est appele notation de Voigt modifie.
Avec cette notation, le produit doublement contract devient
un produit scalaire entre vecteurs ou entre un vecteur et une
matrice.
Par exemple, : = I I
Pour le produit entre un tenseur dordre 4 et un tenseur dordre
2, cela quivaut au produit matriciel: A : = AIJ J .
Rappels
Dfinitions
Oprations
Symtries
Symtries
Reprsentation
Conventions
Symtries mat.
Projecteurs.
Const. Ingnieur
Les mmes rgles sappliquent
Les tenseurs identits sont:
1
1
I=
i = ;
pour la drivation.
1
0
0
0
0
0
0
1
0
0
0
0
0
0
1
0
0
0
0
0
0
1
0
0
0
0
0
0
1
0
0
0
0
0
0
1
(29)
Conventions
Dans la littrature scientifique actuelle, la convention est
dutiliser des lettres romaines majuscules pour les tenseurs
dordre 4 (e.g. A, B, C, etc), des lettres grecques minuscules
pour les tenseurs dordre 2 (e.g. , , , etc) et des lettres
romaines minuscules (e.g. a, b, c, etc.) pour les tenseurs dordre
1.
Lordre dun tenseur peut tre reprsent par le nombre de
barres (ou de tildes dans certains ouvrages) en-dessous ou
en-dessus des lettres. Cest ce qui a t adopt jusqu
maintenant.
Pour des raison de typographie, dans les journaux scientifiques,
on utilise les caractres gras pour reprsenter un tenseur. Par
exemple, a, , et A reprsenteraient des tenseurs dordre 1, 2
et 4, respectivement, tandis que a, et A reprsenteraient des
scalaires (car ils ne sont pas en gras).
Dans le cours, nous allons utiliser indpendamment la notation
avec les barres en-dessous (pour le tableau essentiellement) et
les lettres grasses pour le polycopi.
Rappels
Dfinitions
Oprations
Symtries
Symtries
Reprsentation
Conventions
Symtries mat.
Projecteurs.
Const. Ingnieur
32 / 95
Introduction
Rappels
Certains monocristaux ont une microstructure qui prsente
des plans de symtrie (i.e. les atomes sont disposs de sorte
que lon peut tablir un plan de symtrie)
Certains composites ont les fibres alignes selon certaines
directions, ce qui fait que lon peut avoir soit des plans ou
des axes de symtrie.
Oprations
Symtries
Projecteurs.
Const. Ingnieur
Certains matriaux ont des microstructures qui prsentent
certaines symtries.
Dfinitions
Symtries mat.
Introduction
Triclinique
Monoclinique
Orthotropie
Quadratique
Cubique
Rhombodrique
Iso. Transverse
Isotropie
33 / 95
Ces microstructures vont avoir une effet sur les matrices de
rigidit et de souplesse, notes C et S respectivement
Ces tenseurs dordre 4 vont prendre des formes particulires
et on verra que le nombre de paramtres indpendants pour
les caractriser diminue.
Triclinique
Rappels
Dfinitions
Oprations
Symtries
Symtries mat.
Introduction
Triclinique
Monoclinique
Orthotropie
Quadratique
Cubique
Rhombodrique
Iso. Transverse
Isotropie
Projecteurs.
Const. Ingnieur
34 / 95
Un matriau triclinique est un matriau qui ne prsente aucune
symtrie dans sa microstructure
Le matriau aura donc 21 constantes qui sont distinctes
En fait, il est possible de montrer que le nombre de
contantes non nulles peut tre rduit 18 par un
changement de base appropri.
Triclinique
35 / 95
Exemple
Rappels
Dfinitions
Oprations
Symtries
Symtries mat.
Introduction
Triclinique
Monoclinique
Orthotropie
Quadratique
Cubique
Rhombodrique
Iso. Transverse
Isotropie
Projecteurs.
Const. Ingnieur
Figure 3: Exemple de matriau triclinique. Fibres (f) longues
dorientation quelconque dans une matrice (m) isotrope. Ef = 300GPa,
Em = 70GPa, s = m = 0.3, fv = 30%. C obtenue par calcul lments
finis dune cellule de base sur laquelle on a appliqu des conditions aux
rives priodiques. Gracieuset de Amine El-Mourid.
Monoclinique
Rappels
Dfinitions
Oprations
Symtries
Symtries mat.
Introduction
Triclinique
Monoclinique
Orthotropie
Quadratique
Cubique
Rhombodrique
Iso. Transverse
Isotropie
Projecteurs.
Const. Ingnieur
36 / 95
Un matriau monoclinique prsente un plan de symtrie dans sa
microstructure
Cela veut dire que si on exprime ce tenseur dans une base
dont un vecteur est selon la normale au plan de symtrie, ce
tenseur aura la mme expression dans une base qui est la
rflexion de la premire
Supposons par exemple un tenseur A exprim dans la base
forme par les vecteurs {e1 , e2 , e3 } (qui sont en fait
orients selon les axes x, y et z). Si A est monoclinique et
que son plan de symtrie est normal e2 , alors, il aura la
mme expression dans la base {e1 , e2 , e3 }.
La matrice de passage de la premire base la seconde est:
1
0 0
(30)
[P ] = 0 1 0
0
0 1
Monoclinique
Rappels
Dfinitions
Supposons A sous la forme gnrale suivante (en se rappelant
quil possde les symtries majeures et mineures):
Oprations
A=
Symtries
Symtries mat.
Introduction
Triclinique
Monoclinique
Orthotropie
Quadratique
Cubique
Rhombodrique
Iso. Transverse
Isotropie
37 / 95
Projecteurs.
Const. Ingnieur
A11
A12
A13
A14
A15
A16
A12
A22
A23
A24
A25
A26
A13
A23
A33
A34
A35
A36
A14
A24
A34
A44
A45
A46
A15
A25
A35
A45
A55
A56
A16
A26
A36
A46
A56
A66
Si on applique la relation de changement de base
Aijkl = Pim Pjn Pko Plp Amnop , en notant que
P = (1) (pas de somme sur ) et en appliquant les
proprits du produit simplement contract de la fonction de
Dirac, on a:
Aijkl = (1)i+j+k+l Aijkl (pas de somme sur indices)
(31)
Pour que A soit monoclinique, il faut que Aijkl = Aijkl
(32)
Monoclinique
Effectuons les vrifications en supposant des termes non-nuls:
Rappels
Dfinitions
Oprations
Symtries
Symtries mat.
Introduction
Triclinique
Monoclinique
Orthotropie
Quadratique
Cubique
Rhombodrique
Iso. Transverse
Isotropie
Projecteurs.
Const. Ingnieur
38 / 95
A11 = A1111 = (1)4 A1111 = A1111 = A11
A12 = A1122 = (1)6 A1122 = A1122 = A12
A13 = A1133 = (1)8 A1133 = A1133 = A13
A14 = 2A1123 6= (1) 2A1123 = 2A1123 = A14
A15 = 2A1131 = (1) 2A1131 = 2A1131 = A15
A16 = 2A1112 6= (1) 2A1112 = 2A1112 = A16
A22 = A2222 = (1)8 A2222 = A2222 = A22
A23 = A2233 = (1)10 A2233 = A2233 = A23
A24 = 2A2223 6= (1) 2A2223 = 2A2223 = A24
A25 = 2A2231 = (1) 2A2231 = 2A2231 = A25
A26 = 2A2212 6= (1) 2A2212 = 2A2212 = A26
A34
A33 = A3333 = (1)12 A3333 = A3333 = A33
11
= 2A3323 6= (1)
2A3323 = 2A3323 = A34
Monoclinique
Rappels
Dfinitions
Oprations
Suite des vrifications
10
A35 = 2A3331 = (1)
2A3331 = 2A3331 = A35
A36 = 2A3312 6= (1) 2A3312 = 2A3312 = A36
A44 = 2A2323 = (1)10 2A2323 = A2323 = A44
Symtries
Symtries mat.
Introduction
Triclinique
Monoclinique
Orthotropie
Quadratique
Cubique
Rhombodrique
Iso. Transverse
Isotropie
A45 = 2A2331 6= (1)9 2A2331 = 2A2331 = A45
A55 = 2A3131 = (1)8 2A3131 = 2A3131 = A55
A66 = 2A1212 = (1)6 2A1212 = 2A1212 = A66
A46 = 2A2312 = (1)8 2A2312 = 2A2312 = A46
A56 = 2A3112 6= (1)7 2A3112 = 2A3112 = A56
Projecteurs.
Const. Ingnieur
39 / 95
Bas sur ces rsultats, on peut voir que pour avoir que
Aijkl = Aijkl , il faut que:
A14 = A16 = A24 = A26 = A34 = A36 = A45 = A56 = 0
(33)
Monoclinique
Rappels
Dfinitions
Oprations
Symtries
Symtries mat.
Introduction
Triclinique
Monoclinique
Orthotropie
Quadratique
Cubique
Rhombodrique
Iso. Transverse
Isotropie
Projecteurs.
Const. Ingnieur
40 / 95
Alors, si A est monoclinique et que le plan de symtrie a sa
normale selon e2 , alors la forme gnrique de A est:
A11 A12 A13 0 A15 0
A12 A22 A23 0 A25 0
A13 A23 A33 0 A35 0
A=
(34)
0
0
0 A44 0 A46
A15 A25 A35 0 A55 0
0
0
0 A46 0 A66
Pour le cas o le plan de symtrie a sa normale selon e1 ou e3 ,
il suffit dintervertir certaines lignes et colonnes pour les indices
4 6.
On peut voir que lon aura 13 constantes non-nulles. On peut
montrer que lon peut ramener cela 12 par un changement de
base appropri, ce qui fait que lon a en fait 12 constantes
indpendantes.
Monoclinique
Rappels
Dfinitions
Oprations
Symtries
Symtries mat.
Introduction
Triclinique
Monoclinique
Orthotropie
Quadratique
Cubique
Rhombodrique
Iso. Transverse
Isotropie
Projecteurs.
Const. Ingnieur
41 / 95
Si un matriau monoclinique de souplesse S est soumis un
chargement axial dans la direction e1 , alors la dformation est
donne par:
11
S11 11
S12 11 22
S13 11 33
=
(35)
S:=
0
0
S15 11 231
0
0
On peut voir que contrairement un matriau isotrope,
22 6= 33 et quen plus, malgr que lon ait un chargement
uni-axial, une dformation de cisaillement 31 apparat.
Un exemple dun matriau monoclinique serait par exemple une
range de fibres dans un composite qui seraient comprises dans
un plan mais qui ne serait pas alignes selon les autres axes.
Monoclinique
42 / 95
Exemple
Rappels
Dfinitions
Oprations
Symtries
Symtries mat.
Introduction
Triclinique
Monoclinique
Orthotropie
Quadratique
Cubique
Rhombodrique
Iso. Transverse
Isotropie
Projecteurs.
Const. Ingnieur
Figure 4: Exemple de matriau monoclinique. Fibres (f) longues alignes
et faisant un angle dans le plan xz dans une matrice (m) isotrope. Plan
se symtrie selon la direction y. Ef = 300GPa, Em = 70GPa,
s = m = 0.3, fv = 30%. C obtenue par calcul lments finis dune
cellule de base sur laquelle on a appliqu des conditions aux rives
priodiques. Gracieuset de Amine El-Mourid.
Orthotrope
43 / 95
Un matriau orthotrope a une microstructure qui prsente deux
plans de symtrie.
Si on suppose que les vecteurs e1 et e3 sont les normales ces
plans de symtrie, il faut que A soit invariant par changement
de base avec les matrices de passage suivantes:
1 0 0
1 0
0
h
i
h
i
P (1) = 0 1 0
P (2) = 0 1
0 (36)
0 0 1
0 0 1
Rappels
Dfinitions
Oprations
Symtries
Symtries mat.
Introduction
Triclinique
Monoclinique
Orthotropie
Quadratique
Cubique
Rhombodrique
Iso. Transverse
Isotropie
Projecteurs.
Const. Ingnieur
(1)
o on aura que Pij = (1)n ij , o n = 1 pour i = j = 1 et
(2)
n = 0 autrement, et Pij = (1)m ij , o m = 1 pour
i = j = 3 et m = 0 autrement.
Orthotrope
Rappels
Dfinitions
Oprations
Symtries
Symtries mat.
Introduction
Triclinique
Monoclinique
Orthotropie
Quadratique
Cubique
Rhombodrique
Iso. Transverse
Isotropie
Projecteurs.
Const. Ingnieur
44 / 95
Pour le matriau orthotrope, on doit vrifier les deux quations
suivantes:
(1)
(1)
(1)
(1)
Aijkl = Pim Pjn Pko Plp Amnop
w
= (1) Aijkl
(2)
(2)
(2)
(37a)
(2)
Aijkl = Pim Pjn Pko Plp Amnop
z
= (1) Aijkl
(37b)
o w = 1 pour un nombre impair dindices gaux 1 et w = 0
autrement, z = 1 pour un nombre impair dindices gaux 3 et
z = 0 autrement.
En effectuant les vrifications comme prcdemment ( faire en
exercice), on trouve que lon doit avoir:
A14 = A15 = A16 = A24 = A25 = A26 = A34 =
A35 = A36 = A45 = A46 = A56 = 0 (38)
Orthotrope
Rappels
Dfinitions
Oprations
Symtries
Symtries mat.
Introduction
Triclinique
Monoclinique
Orthotropie
Quadratique
Cubique
Rhombodrique
Iso. Transverse
Isotropie
Projecteurs.
Const. Ingnieur
45 / 95
Alors, si A est orthotrope la
A11 A12
A12 A22
A13 A23
A=
0
0
0
0
0
0
forme gnrique de A est:
A13 0
0
0
A23 0
0
0
A33 0
0
0
0 A44 0
0
0
0 A55 0
0
0
0 A66
(39)
On peut voir que si on avait pris deux autres plans de symtrie,
on aurait obtenu le mme rsultat. Cela revient dire que le
fait davoir deux plans de symtrie dans la microstructure
implique la prsence dun troisime.
On voit quun matriau orthotrope est dfini par 9 constantes.
Un exemple de matriau est le mme composite que lexemple
monoclinique mais o les fibres sont orientes selon un des axes
de la base.
Orthotrope
46 / 95
Exemple
Rappels
Dfinitions
Oprations
Symtries
Symtries mat.
Introduction
Triclinique
Monoclinique
Orthotropie
Quadratique
Cubique
Rhombodrique
Iso. Transverse
Isotropie
Projecteurs.
Const. Ingnieur
Figure 5: Exemple de matriau orthotrope. Fibres (f) longues isotropes
de modules diffrents disposes perpendiculairement dans une matrice
(m) isotrope. Dans la direction x, Ef = 3000GPa et f = 0.43; dans la
direction y, Ef = 300GPa et f = 0.26; dans la direction z,
Ef = 1500GPa et f = 0.43. Em = 70GPa et m = 0.3. fv = 17.45%.
C obtenue par calcul lments finis dune cellule de base sur laquelle on a
appliqu des conditions aux rives priodiques. Gracieuset de Amine
El-Mourid.
Quadratique
47 / 95
Si la microstructure dun matriau est invariante par rotations
de k 2 autour dun axe, on dit quil a la symtrie quadratique.
Par exemple, si on tudie les rotations autour de z, la matrice
de passage devient:
cos k 2 sin k 2 0
[P ] = sin k 2
(40)
cos k 2 0
0
0 1
Rappels
Dfinitions
Oprations
Symtries
Symtries mat.
Introduction
Triclinique
Monoclinique
Orthotropie
Quadratique
Cubique
Rhombodrique
Iso. Transverse
Isotropie
Projecteurs.
Const. Ingnieur
On peut donc dfinir trois matrices de passage pour
k = {1, 2, 3} et on obtiendra:
1
0 0
0 1 0
h
i
h
i
P (2) = 0 1 0
P (1) = 1 0 0
0
0 1
0 0 1
(41)
0 1 0
h
i
0 0
P (3) = 1
0
0 1
Quadratique
Rappels
(z)
Oprations
Symtries mat.
Introduction
Triclinique
Monoclinique
Orthotropie
Quadratique
Cubique
Rhombodrique
Iso. Transverse
Isotropie
Projecteurs.
Const. Ingnieur
Pour que A possde la symtrie quadratique, il faut que:
(z)
(z)
(z)
Aijkl = Pim Pjn Pko Plp Amnop
Dfinitions
Symtries
48 / 95
Si on effectue les calculs, on
de la forme:
A11 A12
A12 A11
A13 A13
A=
0
0
0
0
0
0
pour z = {1, 2, 3}
(42)
se rendra compte que A doit tre
A13 0
0
0
A13 0
0
0
A33 0
0
0
0 A44 0
0
0
0 A44 0
0
0
0 A66
(43)
On pourra interprter un tenseur quadratique comme un tenseur
orthotrope pour lequel deux directions sont quivalentes.
La forme de A peut facilement se dduire les matriaux tant
invariants par rotation de k 2 autour des autres axes.
On a maintenant 6 composantes indpendantes.
Quadratique
49 / 95
Exemple
Rappels
Dfinitions
Oprations
Symtries
Symtries mat.
Introduction
Triclinique
Monoclinique
Orthotropie
Quadratique
Cubique
Rhombodrique
Iso. Transverse
Isotropie
Projecteurs.
Const. Ingnieur
Figure 6: Exemple de matriau quadratique. Fibres (f) longues isotropes
alignes disposes selon un arrangement carr dans une matrice (m)
isotrope. Ef = 300GPa, f = 0.3, Em = 70GPa et m = 0.3. fv = 30%.
C obtenue par calcul lments finis dune cellule de base sur laquelle on a
appliqu des conditions aux rives priodiques. Gracieuset de Amine
El-Mourid.
Cubique
50 / 95
Rappels
Dfinitions
Oprations
Pour que A possde la symtrie cubique, il faut quil soit
invariant par rotation de k 2 autour des trois axes. La matrice
de passage pour cette transformation est obtenue par rotation
successive autour des axes x, y et z et est donne par:
Symtries
Symtries mat.
Introduction
Triclinique
Monoclinique
Orthotropie
Quadratique
Cubique
Rhombodrique
Iso. Transverse
Isotropie
Projecteurs.
Const. Ingnieur
1
0
0
[P ] = 0 cos 2 sin 2
0 sin 2
cos 2
cos 2 0 sin 2
0
1
0
sin 2 0 cos 2
cos 2 sin 2 0
sin
cos
0 (44)
2
2
0
0 1
o le () est le produit matriciel (i.e. simplement contract) et
o , et peuvent prendre les valeurs 1, 2 et 3.
Cubique
51 / 95
Il faut donc que A soit invariant par changement de base avec
les 27 matrices de passage qui seront gnres partir de
lquation (44).
Si on effectue les calculs, on obtiendra
Rappels
Dfinitions
Oprations
Symtries
Symtries mat.
Introduction
Triclinique
Monoclinique
Orthotropie
Quadratique
Cubique
Rhombodrique
Iso. Transverse
Isotropie
Projecteurs.
Const. Ingnieur
A=
A11 A12 A12 0
0
0
A12 A11 A12 0
0
0
A12 A12 A11 0
0
0
0
0
0 A44 0
0
0
0
0
0 A44 0
0
0
0
0
0 A44
(45)
On aura 3 constantes indpendantes.
La symtrie cubique peut-tre interprte comme lorthotropie
o les trois directions sont quivalentes.
Cubique
52 / 95
Exemple
Rappels
Dfinitions
Oprations
Symtries
Symtries mat.
Introduction
Triclinique
Monoclinique
Orthotropie
Quadratique
Cubique
Rhombodrique
Iso. Transverse
Isotropie
Projecteurs.
Const. Ingnieur
Figure 7: Exemple de matriau composite symtrie cubique. Sphres
(s) lastiques et isotropes distribues selon un arrangement rgulier et
cubique dans une matrice (m) lastique et isotrope. Es = 300GPa,
Em = 70GPa, s = m = 0.3, fv = 30%. C obtenue par calcul lments
finis dune cellule de base sur laquelle on a appliqu des conditions aux
rives priodiques. Gracieuset de Amine El-Mourid.
Rhombodrique
53 / 95
Similairement la symtrie quadratique, un tenseur est
rhombodrique lorsque quil est invariant par des changement
de base correspondant des rotations de k 2
3 autour dun seul
axe.
Si par exemple laxe de rhombodrie est laxe z, on aura les
matrices de passage suivantes:
2
2
cos k 3 sin k 3 0
2
cos
k
0
(46)
[P ] = sin k 2
3
3
0
0 1
Rappels
Dfinitions
Oprations
Symtries
Symtries mat.
Introduction
Triclinique
Monoclinique
Orthotropie
Quadratique
Cubique
Rhombodrique
Iso. Transverse
Isotropie
Projecteurs.
Const. Ingnieur
o k = {1, 2, 3, 4, 5}. Il faut donc que A ait la mme
expression dans les 5 bases obtenues laide des matrices de
passage dfinies lquation (46).
Rhombodrique
Rappels
Dfinitions
Oprations
Symtries
Symtries mat.
Introduction
Triclinique
Monoclinique
Orthotropie
Quadratique
Cubique
Rhombodrique
Iso. Transverse
Isotropie
Projecteurs.
Const. Ingnieur
Si on effectue les calculs,
A11 A12 A13
A12 A11 A13
A13 A13 A33
A=
A14 A14 0
0
0
0
0
0
0
54 / 95
on aura
A14
0
0
A14
0
0
0
0
0
A44
0
0
2A14
0
A44
0
2A14 A11 A12
On aura donc 6 composantes indpendantes.
(47)
Isotropie transverse
Rappels
Dfinitions
Oprations
Symtries
Symtries mat.
Introduction
Triclinique
Monoclinique
Orthotropie
Quadratique
Cubique
Rhombodrique
Iso. Transverse
Isotropie
55 / 95
Similairement aux symtries cubiques et rhombodriques, un
tenseur isotrope transverse A est invariant par une rotation
arbitraire autour dun axe donn. Par exemple, si laxe
disotropie transverse est laxe z, alors, la matrice de passage
est:
cos sin 0
[P ] = sin
cos 0
(48)
0
0 1
pour tout angle .
En rptant les mmes oprations que pour les autres symtries,
on trouvera quun tenseur isotrope transverse est de la forme:
Projecteurs.
Const. Ingnieur
A=
A11 A12 A13 0
0
0
A12 A11 A13 0
0
0
A13 A13 A33 0
0
0
0
0
0 A44 0
0
0
0
0
0 A44
0
0
0
0
0
0 A11 A12
o on aura 5 composantes indpendantes.
(49)
Isotropie
56 / 95
Un tenseur A est isotrope sil a la mme expression dans toutes
les bases qui sont obtenues par des rotations arbitraires autour
des tois axes
La matrice de passage est comme celle de lquation (44) mais
o les 2 ont t enlevs.
En rptant les mmes oprations que pour les autres symtries,
on trouvera quun tenseur isotrope transverse est de la forme:
A11 A12 A12
0
0
0
A12 A11 A12
0
0
0
A12 A12 A11
0
0
0
A=
0
0 A11 A12
0
0
0
0
0
0
A11 A12
0
0
0
0
0
0
A11 A12
(50)
Rappels
Dfinitions
Oprations
Symtries
Symtries mat.
Introduction
Triclinique
Monoclinique
Orthotropie
Quadratique
Cubique
Rhombodrique
Iso. Transverse
Isotropie
Projecteurs.
Const. Ingnieur
On aura donc 2 constantes indpendantes.
Projecteurs isotropes
57 / 95
Si un tenseur dordre 4 est isotrope, on peut le dcomposer en:
Rappels
Oprations
Symtries
Symtries mat.
Projecteurs.
Isotropie
Cubique
Iso. Transverse
Const. Ingnieur
(51)
A = J + K
Dfinitions
o J = 13 i i et K = I J, avec i et I les tenseurs identit
dordre 1 et 2.
(voir dmonstration au tableau o lon montre que J et
K sont isotropes)
Si on fait les calculs, on aura que:
2
6
16
6
J= 6
36
4
1
1
1
0
0
0
1
1
1
0
0
0
1
1
1
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
7
6
7
16
7
6
7 et K = 6
7
36
5
4
2
1
1
0
0
0
1
2
1
0
0
0
1
1
2
0
0
0
0
0
0
3
0
0
0
0
0
0
3
0
0
0
0
0
0
3
7
7
7
7 (52)
7
5
Projecteurs isotropes
58 / 95
En utilisant lquation (52) et si lon fait les calculs, on a le
rsultat suivant: J : J = J
Il en dcoule que:
Rappels
Dfinitions
Oprations
Symtries
Symtries mat.
Projecteurs.
Isotropie
Cubique
Iso. Transverse
Const. Ingnieur
J : K = J : (I J) = J J = 0
K : K = K : (I J) = K 0 = K
Les tenseurs J et K sont qualifis dorthogonaux.
Pour un tenseur de rigidit isotrope, on aura que:
C = 3kJ + 2K
(53)
o k est le module de compressibilit et le module de
cisaillement.
Pour un tenseur de souplesse isotrope, on aura que:
S=
1
1
J+
K
3k
2
(voir la dmonstration au tableau)
(54)
Projecteurs isotropes
Rappels
Oprations
Symtries
Symtries mat.
Projecteurs.
Isotropie
Cubique
Iso. Transverse
Const. Ingnieur
Il dcoule de ces observations les proprits suivantes. Soient A
et A deux tenseurs dordre 4. Alors:
Dfinitions
59 / 95
A + A = ( + ) J + ( + ) K
A : A = ( ) J + ( ) K
A = J + K
A1 = 1 J + 1 K
Soit A = J + K, alors:
= J :: A
5 = K :: A
(voir dmonstrations au tableau)
Projecteurs isotropes
Rappels
Dfinitions
Oprations
Symtries
Symtries mat.
Projecteurs.
Isotropie
Cubique
Iso. Transverse
Const. Ingnieur
60 / 95
Soit le tenseur des dformations, alors
1
J : = sph () = tr () i
3
1 0 0
11 + 22 + 33
0 1 0
=
3
0 0 1
(55)
( dmontrer en exercice)
sph () est appele partie sphrique de .
Soit le tenseur des dformations, alors
1
K : = dev () = tr () i
3
#
"
2
3
3
11
22
33
12
13
1
312
222 11 33
323
=
3
313
323
233 11 22
(56)
( dmontrer en exercice)
dev () est appele partie dviatorique de .
Projecteurs isotropes
61 / 95
Interprtation
Rappels
On a que
Dfinitions
Oprations
Symtries
Symtries mat.
Projecteurs.
Isotropie
Cubique
Iso. Transverse
Const. Ingnieur
tr () = 11 + 22 + 22
x y z
+
+
=
L
L
L
(57)
o x, y et z sont les allongements infinitsimaux dans les
directions x, y et z dun cube de longueur L.
V0
Si on calcule le changement de volume relatif, V V
, on aura
0
que:
V V0
(L + x) (L + z) (L + z) L3
=
V0
L3
L2 (x + y + z)
+
=
L3
L (xy + xz + yz)
+
3
L
xyz
L3
(58)
Projecteurs isotropes
Comme x, y et z sont trs petits devant L, alors:
x + y + z
V V0
V0
L
tr ()
Rappels
Dfinitions
Oprations
Symtries
Symtries mat.
Projecteurs.
Isotropie
Cubique
Iso. Transverse
62 / 95
Il est facile de montrer que:
V V0
tr (sph ()) = tr ()
V0
tr (dev ()) = 0
Const. Ingnieur
(59)
(60a)
(60b)
On peut donc voir que J permet dextraire la partie dun
tenseur dordre 2 qui correspond au changement de volume,
mais forme constante. Cest pour cette raison que dans
lexpression de C, le module de compressibilit k multiplie J.
On peut donc voir que K permet dextraire la partie dun
tenseur dordre 2 qui correspond au changement de forme
volume constant. Cest pour cette raison que dans lexpression
de C, le module de cisaillement multiplie K.
Projecteurs isotropes
On aussi vu plus haut que J et K permettent dextraire les
paramtres et (k et quand A = C).
Il est intressant de remarquer que k et multiplient deux
tenseurs orthogonaux. Cela permet dinterprter k et comme
des constantes associes deux modes de dformation
indpendants.
Ce nest pas le cas avec le module dYoung E et le coefficient
de Poisson qui font appel avec des dformations qui changent
la fois la forme et le volume (sauf pour = 0.5).
Pour ces raisons, on va prfrer travailler dans ce cours avec
(k, ) quavec (E, ) pour caractriser des matriaux isotropes.
On verra que cela simplifie de beaucoup certains calculs.
Rappels
Dfinitions
Oprations
Symtries
Symtries mat.
Projecteurs.
Isotropie
Cubique
Iso. Transverse
Const. Ingnieur
63 / 95
Projecteurs pour la symtrie cubique
64 / 95
Dfinissons le tenseur
Rappels
Dfinitions
Oprations
Symtries
Symtries mat.
Projecteurs.
Isotropie
Cubique
Iso. Transverse
Const. Ingnieur
Z = e1 e1 e1 e1 + e2 e2 e2 e2
+ e3 e3 e3 e3 (61)
o on peut voir que ce tenseur a la symtrie cubique. En effet,
toutes les rotations de 2 autour dun axe vont faire en sorte
que les deux autres vecteurs vont tre interchangs.
Par exemple, si on effectue une rotation autour de e3 de 2 , on
aura que e1 = e2 et que e2 = e1 . En appliquant la dfinition
de Z, on se rend compte que Z est invariant par les rotations
de 2 , donc quil est de symtrie cubique.
Comme J et K sont isotropes, ils possdent aussi la symtrie
cubique.
Pour la symtrie cubique, on introduit les tenseurs suivants:
Ka = Z J et Kb = K Ka = I Z
(62)
Projecteurs pour la symtrie cubique
65 / 95
On peut montrer que J, Ka et Kb sont orthogonaux (voir
preuve au tableau).
Il en dcoule quun tenseur A symtrie cubique peut
sexprimer par:
A = J + Ka + Kb
(63)
Rappels
Dfinitions
Oprations
Symtries
Symtries mat.
Projecteurs.
Isotropie
Cubique
Iso. Transverse
Const. Ingnieur
On aura les mmes proprits (addition, inverse, double
contraction, multiplication par un scalaire) que pour le tenseur
isotrope, ce qui peut tre trs utile pour les calculs.
Pour interprter les tenseurs J, Ka et Kb , nous allons procder
de la mme manire que pour lisotropie. Soit le tenseur des
dformations.
Comme pour lisotropie, J : = sph () i
1
Ka : =
3
"
211 22 33
0
0
0
222 11 33
0
0
0
233 11 22
(64)
Projecteurs pour la symtrie cubique
Rappels
Dfinitions
Oprations
Symtries
Symtries mat.
Projecteurs.
Isotropie
Cubique
Iso. Transverse
Const. Ingnieur
66 / 95
0 12 13
1
12 0 23
Kb : =
3
13 23 0
(65)
On peut donc voir que Ka extrait la partie axiale de la partie
dviatorique des dformations tandis que Kb extrait les de
cisaillement de la partie dviatorique du tenseur des
dformations.
Alors, si C est un tenseur de rigidit symtrie cubique, on
aura que:
C = 3kJ + 2a Ka + 2b Kb
(66)
o a et b pourraient tre interprts comme deux modules de
cisaillement: a pour le changement de forme dun cube un
prisme et b pour le changement des angles droits du prisme.
De manire similaire lisotropie, on aura les relations
suivantes:
J :: A = ;
Ka :: A = 2;
Kb :: A = 3
(67)
Projecteurs pour lisotropie transverse
67 / 95
Des relations similaires existent pour lisotropie transverse.
Soit n, le vecteur disotropie transverse. On dfinit iT le
tenseur didentit dans le plan transverse par:
Rappels
Dfinitions
Oprations
iT = i n n
Symtries
Symtries mat.
Projecteurs.
Isotropie
Cubique
Iso. Transverse
Const. Ingnieur
Par exemple, si n = e3 , on aura:
0 0 0
1 0 0
iT = 0 1 0 0 0 0
0 0 1
0 0 1
1 0 0
= 0 1 0
0 0 0
(68)
(69)
On voit que iT est i laquelle on a retranch la composante
correspondant la direction de laxe disotropie transverse.
Avec leur expressions matricielles, on peut facilement voir que
iT iT = iT et que (n n) iT = 0.
Projecteurs pour lisotropie transverse
Rappels
Dfinitions
Oprations
Symtries
Symtries mat.
Projecteurs.
Isotropie
Cubique
Iso. Transverse
Const. Ingnieur
68 / 95
Introduisons le tenseur EL = n n n n
Ce tenseur est isotrope transverse. Pour le montrer,
considrons n = e3 . On peut voir que ELijkl = 0 sauf pour
EL3333 = 1. Considrons le tenseur EL exprim dans une
base qui est la base dorigine laquelle on a fait subir une
rotation arbitraire autour de n. Alors, on aura:
EL ijkl = Pim Pjn Pko Plp ELmnop
= Pi3 Pj3 Pk3 Pl3 EL3333
(70)
Avec la dfinition de Pij de lquation (48), on voit que
Pi3 = 0, sauf pour P33 = 1. Alors, EL ijkl = 0, sauf pour
i = j = k = l = 3 o il vaut EL3333 . On vient donc de
prouver que EL est isotrope transverse.
Projecteurs pour lisotropie transverse
Rappels
69 / 95
Introduisons le tenseur JT = 21 iT iT
Dfinitions
Oprations
Ce tenseur est isotrope transverse. Pour le montrer,
considrons n = e3 . En utilisant la dfinition de iT , on a
que:
Symtries
Symtries mat.
Projecteurs.
Isotropie
Cubique
Iso. Transverse
Const. Ingnieur
1
JT = i i i (n n) (n n) i
2
+nnnn
1
= 3J + EL i (n n) (n n) i
2
(71)
On a dj montr que J et EL taient au moins isotropes
transverse.
Montrez que i (n n) + (n n) i est isotrope
transverse. faire en exercice
Projecteurs pour lisotropie transverse
Rappels
Dfinitions
Oprations
Symtries
Symtries mat.
Projecteurs.
Isotropie
Cubique
Iso. Transverse
Const. Ingnieur
70 / 95
Introduisons le tenseur IT . Ce tenseur permet dextraire les
composantes dans le plan disotropie transverse dun tenseur
dordre deux. Par exemple, si e3 est laxe disotropie transverse,
alors:
1 0 0 0 0 0
0 1 0 0 0 0
0 0 0 0 0 0
IT =
(72)
0 0 0 0 0 0
0 0 0 0 0 0
0 0 0 0 0 1
Si on effectue les calculs, on aura que:
11 12 0
IT : = 12 22 0
0
0 0
ce qui est bien le rsultat voulu.
(73)
Projecteurs pour lisotropie transverse
Rappels
Dfinitions
Oprations
Symtries
Symtries mat.
Projecteurs.
Isotropie
Cubique
Iso. Transverse
Const. Ingnieur
71 / 95
Introduisons les tenseurs suivants:
1
KE = (2n n iT ) (2n n iT )
6
KT = IT JT
(74a)
(74b)
KL = K KT KE
(74c)
2
(iT n n)
(74d)
F=
2
o F nest pas symtrique.
On peut montrer que tout tenseur A isotrope transverse peut
scrire sous la forme
A = EL + JT + F + FT + KT + KL
o = si A est symtrique. Cette relation gnrale est
prsente car il est possible que pour certains produits
doublement contracts, le tenseur rsultant ne soit plus
symtrique.
(75)
Projecteurs pour lisotropie transverse
Les tenseurs EL , JT , KE , KT et KL sont des projecteurs
orthogonaux, comme J et K, par exemple.
Soit le tenseur des dformations. Nous allons interprter ces
diffrents tenseurs en considrant leur effet sur . On considre
toujours que laxe disotropie transverse est e3 .
On peut voir que:
0 0 0
(76)
EL : = 0 0 0
0 0 33
Rappels
Dfinitions
Oprations
Symtries
Symtries mat.
Projecteurs.
Isotropie
Cubique
Iso. Transverse
Const. Ingnieur
72 / 95
On peut voir que EL a leffet dextraire la dformation
uniaxiale selon laxe disotropie transverse.
Projecteurs pour lisotropie transverse
Rappels
Symtries
Const. Ingnieur
1 0 0
11 + 22
0 1 0
JT : =
2
0 0 0
Oprations
Projecteurs.
Isotropie
Cubique
Iso. Transverse
Si on fait les calculs, on aura que:
Dfinitions
Symtries mat.
73 / 95
(77)
On peut voir que JT extrait le changement de surface dans le
plan transverse.
Si on fait les calculs, on aura que:
1 0 0
233 11 22
KE : =
(78)
0 1 0
6
0
0 2
On peut voir que KE extrait lallongement volume constant
selon laxe disotropie transverse.
Projecteurs pour lisotropie transverse
Rappels
Dfinitions
Oprations
Symtries
Symtries mat.
Projecteurs.
Isotropie
Cubique
Iso. Transverse
Const. Ingnieur
Si on fait les calculs, on aura que:
22
212
0
1 11
212
22 11 0
KT : =
2
0
0
0
74 / 95
(79)
On peut voir que KT extrait les composantes de cisaillement
dans le plan transverse ( surface constante).
Si on fait les calculs, on aura que:
0
0 13
0 23
(80)
KL : = 0
13 23 0
On peut voir que KL extrait les composantes de cisaillement
faisant intervenir laxe disotropie transverse.
Projecteurs pour lisotropie transverse
Rappels
Dfinitions
Oprations
Symtries
Symtries mat.
Projecteurs.
Isotropie
Cubique
Iso. Transverse
Const. Ingnieur
Si on fait les calculs, on aura que:
33 0
0
1
33
0
F + FT : =
2
0
0 11 + 22
75 / 95
(81)
En effectuant les calculs on voit F fait apparatre 33 aux
positions 11 et 22 du tenseur rsultant tandis que FT fait
apparatre 11 et 22 la position 33 du tenseur rsultant.
Le tenseur F fait donc apparatre le couplage quil y a entre
lallongement dans la direction disotropie transverse et celle
dans le plan transverse, et vice-versa.
F nest pas un projecteur orthogonal par rapport aux autres
tenseurs introduits pour dfinir lisotropie transverse. Les rgles
dadditivit et de contraction ne sont pas aussi simples que
celles dcrites pour lisotropie et la symtrie cubique.
Projecteurs pour lisotropie transverse
76 / 95
Si on a deux tenseurs isotropes transverses A(1) et A(2) , alors:
Rappels
A(1) + A(2) = (1 + 2 ) EL
Dfinitions
FT
+ (1 + 2 ) KT + 1 + 2 KL (82)
Oprations
+ (1 + 2 ) JT + (1 + 2 ) F +
Symtries
Symtries mat.
Projecteurs.
Isotropie
Cubique
Iso. Transverse
Const. Ingnieur
Le produit contract devient:
(1)
(2)
:A
= 1 2 + 1 2 EL
1 2 + 1 2 JT + (1 2 + 1 2 ) F
T
+ 1 2 + 1 2 F + (1 2 ) KT
+ 1 2 KL
(83)
Projecteurs pour lisotropie transverse
Rappels
Dfinitions
Oprations
Symtries
Symtries mat.
Projecteurs.
Isotropie
Cubique
Iso. Transverse
Const. Ingnieur
Il est donc commode de reprsenter un tenseur isotrope
transverse par le triplet:
A = , , , =
Les rgles opratoires deviennent donc:
n
o n
o
A(1) + A(2) = (1) , 1 , 1 + (2) , 2 , 2
n
o
= (1) + (2) , 1 + 2 , 1 + 2
A(1) : A(2)
n
o n
o
= (1) , 1 , 1 : (2) , 2 , 2
n
o
= (1) (2) , 1 2 , 1 2
1
1 1 1
A = , ,
77 / 95
(84)
(85a)
(85b)
(85c)
Projecteurs pour lisotropie transverse
Rappels
Dfinitions
Oprations
Symtries
Symtries mat.
Projecteurs.
Isotropie
Cubique
Iso. Transverse
Const. Ingnieur
78 / 95
En effectuant les calculs, on trouvera que:
EL :: A = ;
F :: A = ;
JT :: A = ;
KT :: A = 2;
FT :: A =
KL :: A = 2
(86)
ce qui peut tre trs utile pour trouver la dcomposition dun
tenseur isotrope transverse lorsque lon connat son expression
tensorielle.
Introduction
79 / 95
On a souvent recours aux constantes dingnieur pour dcrire
des tenseurs de rigidit ou souplesse.
Ces constantes sont physiquement plus faciles comprendre
que celles prsentes jusqu maintenant, mais ont des
proprits mathmatiques moins intressantes (on verra
pourquoi plus loin).
On connat trois types de constantes dingnieur
Rappels
Dfinitions
Oprations
Symtries
Symtries mat.
Projecteurs.
Const. Ingnieur
Introduction
Orthotropie
Isotropie
Cubique
Iso. Transverse
Df. Positif
1. Module dYoung Ei dfini par: Ei = iiii lorsque le matriau
est soumis une contrainte ii , toutes les autres tant
nulles.
2. Coefficient de Poisson ij dfini par ij = jj
lorsque le
ii
matriau est soumis une contrainte ii , toutes les autres
tant nulles.
3. Les modules de cisaillement Gij dans le plan ij.
Orthotropie
Rappels
Dfinitions
Oprations
Symtries
Symtries mat.
Projecteurs.
Const. Ingnieur
Introduction
Orthotropie
Isotropie
Cubique
Iso. Transverse
Df. Positif
80 / 95
Avec les constantes dingnieur, la souplesse S dun tenseur
orthotrope est donne par:
21
31
1
E3
0
0
0
E1 E2
32
1
12
0
0
0
E1
E2
E3
13 23
1
0
0
0
E1
E2
E3
S=
(87)
1
0
0
0
0
0 2G23
0
0
0
0
0
2G13
1
0
0
0
0
0
2G12
La rigidit C sera donne par C = S1 (le rsultat nest pas
donn ici car il est relativement complexe).
On verra plus tard que les rigidits et les souplesses sont
toujours symtriques. Cela permet de dduire les rsultats
suivants:
21
12
=
;
E2
E1
31
13
=
;
E3
E1
32
23
=
E3
E2
(88)
Isotropie
81 / 95
Rappels
Dfinitions
Oprations
Symtries
Symtries mat.
Projecteurs.
Const. Ingnieur
Introduction
Orthotropie
Isotropie
Cubique
Iso. Transverse
Df. Positif
Avec les constantes dingnieur, la
isotrope est donne par:
1
E
E
E
E
E
E
E
E
E
S=
0
0
0
0
0
0
0
0
0
souplesse S dun tenseur
0
0
0
1+
E
0
0
0
0
0
0
1+
E
0
0
0
0
0
1+
E
(89)
En se rappelant que S = J + K, que = J :: S et que
5 = K :: S, on aura que (aprs avoir effectu les calculs):
1 2
1
=
=
E
3k
(90a)
1+
1
=
E
2
(90b)
o k et sont les modules de compressibilit et de cisaillement,
respectivement.
Isotropie - note sur le module de compressibilit
Rappels
Dfinitions
Symtries
Symtries mat.
Projecteurs.
Const. Ingnieur
Introduction
Orthotropie
Isotropie
Cubique
Iso. Transverse
Df. Positif
Par dfinition, le module de compressibilit relie la pression
hydrostatique au changement de volume par la relation:
= 3k
Oprations
82 / 95
(91)
o est la pression hydrostatique et est la dformation axiale
dans une des trois directions lorsque le solide est soumis une
pression hydrostatique.
En utilisant la souplesse dun tenseur isotrope soumis une
pression hydrostatique telle ii = , on trouve que
E
k=
3(1 2)
(92)
On peut donc voir que lon peut relier les constantes devant les
projecteurs aux constantes dingnieur.
Symtrie cubique
Rappels
Dfinitions
Oprations
Symtries
Symtries mat.
Projecteurs.
Const. Ingnieur
Introduction
Orthotropie
Isotropie
Cubique
Iso. Transverse
Df. Positif
83 / 95
Avec les constantes dingnieur, la souplesse S
symtrie cubique est donne par:
1
0
0
E
E
E
1
0
0
E
E
E
0
0
E
E
E
S=
1
0
0
0 2G
0
0
0
0 0 2G
0
0
0 0
0
dun tenseur
0
0
0
0
0
1
2G
(93)
En se rappelant que S = J + Ka + Kb , que = J :: S,
que 2 = Ka :: S et que 3 = Kb :: S, on aura que:
1 2
1
=
E
3k
1+
1
=
= a
E
2
1
1
= b
=
2G
2
o k, a et b sont les modules de compressibilit et de
cisaillement.
(94a)
(94b)
(94c)
Isotropie transverse
Rappels
Dfinitions
Oprations
Symtries
Symtries mat.
Projecteurs.
Const. Ingnieur
Introduction
Orthotropie
Isotropie
Cubique
Iso. Transverse
Df. Positif
84 / 95
Avec les constantes dingnieur, la souplesse S dun tenseur
isotrope transverse est donne par:
l
t
1
0
0
0
Et Et El
l
1
t
0
0
0
Et
Et
El
l l
1
0
0
0
El
El
El
S=
(95)
1
0
0
0
0
0 2G
0
0
0
0 0
2Gl
1+t
0
0
0 0
0
Et
o laxe disotropie transverse est e3 , Et = E1 = E2 , El = E3 ,
t = 12 = 21 , l = 31 = 32 , Gl = G23 = G13 .
Il est aussi intressant de dfinir les quantits suivantes:
Et
Le module de cisaillement transverse Gt = 2(1+
t)
Le module de compressibilit transverse tat plan de
l
dformation Kt = 2 (1 E)Et E2
(voir dmonstration
2E
[
]
t
l
l t
au tableau).
Isotropie transverse
Rappels
Dfinitions
85 / 95
En se rappelant que
T
S = EL + JT + F + F + KT + KL et en effectuant
les calculs, on aura que:
Oprations
1
= EL :: S =
El
Symtries
Symtries mat.
Projecteurs.
1 t
= JT :: S =
Et
2l
= = F :: S =
El
Const. Ingnieur
Introduction
Orthotropie
Isotropie
Cubique
Iso. Transverse
Df. Positif
(96b)
(96c)
1 + t
1
= KT :: S =
2
Et
(96d)
1
1
= KL :: S =
2
2Gl
(96e)
(96a)
Il serait intressant dinterprter le lien entre les constantes
dingnieur et les modules utiliss pour dfinir une rigidit
isotrope transverse.
Isotropie transverse
Rappels
Dfinitions
Oprations
Symtries
Symtries mat.
Projecteurs.
Const. Ingnieur
Introduction
Orthotropie
Isotropie
Cubique
Iso. Transverse
Df. Positif
86 / 95
c ,
c,
c,
Si on nomme par
et
les modules dun tenseur
de rigidit isotrope transverse C = S1 et en utilisant la
relation (85c), on obtient:
c
c
1
=
(97a)
2
c c
1
=
(97b)
1
=
En utilisant les rsultats de lquation (96) on a que:
c
c
=
=
2
= El + 4l2 Kt
(97c)
1t
Et
1 1t
El Et
2l2
El2
(98)
Isotropie transverse
Rappels
Dfinitions
Oprations
Symtries
Symtries mat.
Projecteurs.
Const. Ingnieur
Introduction
Orthotropie
Isotropie
Cubique
Iso. Transverse
Df. Positif
87 / 95
On aura de plus:
1
=
=
2
El
1 1t
El Et
2l2
El2
El Et
=
El (1 t ) 2Et l2
(99a)
= 2Kt
c =
=
2
= 2 2l Kt
2l
1 1t
El Et
2l2
El2
(99b)
1
= = 2Gt
(99c)
1
= 2Gl
(99d)
c =
Isotropie transverse
Dfinitions
Oprations
Symtries
Symtries mat.
Projecteurs.
Const. Ingnieur
Introduction
Orthotropie
Isotropie
Cubique
Iso. Transverse
Df. Positif
Il est intressant de remarquer que
et
sont relis aux
modules de cisaillement dun matriau isotrope transverse.
c , qui multiplie JT , fait aussi apparatre un module de
compressibilit, comme dans le cas de lisotropie et de la
symtrie cubique.
Les termes c et c sont plus difficiles interprter car ils font
intervenir un couplage entre El , l et Kt .
Rappels
88 / 95
Tenseurs dfinis positifs
Rappels
Symtries
Symtries mat.
Projecteurs.
Const. Ingnieur
Introduction
Orthotropie
Isotropie
Cubique
Iso. Transverse
Df. Positif
Un tenseur A est dit dfini positif si
: A : > 0
Dfinitions
Oprations
89 / 95
(100)
Il est commun de noter A > 0 pour indiquer quil est dfini
positif, bien que cette relation nait pas de sens.
Lorsque A possde les symtries majeures et mineures, on a vu
que la relation prcdente scrit sous la forme
{}T [A] {} > 0 {}
(101)
o {} est un vecteur 6 composantes et [A] est une matrice
6 6 symtrique.
Lorsque lon a une matrice carre et symtrique, il est toujours
possible de trouver une base dans laquelle cette matrice est
diagonale. Si on nomme par ii les termes sur la diagonale, qui
sont en fait les valeurs propres, la relation (101) devient:
i ii i > 0 ii > 0
(102)
Tenseurs dfinis positifs
Alors, pour dterminer si A > 0, il suffit de calculer ses valeurs
propres et de sassurer quelles soient toutes strictement
positives.
En mcanique, on montrera plus tard que S > 0 et C > 0,
sinon, le matriau nest pas admissible thermodynamiquement.
Si on entre une matrice de rigidit qui nest pas dfinie positive
dans un code de calculs par lments finis, on peut obtenir des
erreurs. Cest ce qui arrive typiquement lorsque lon ne connat
pas toutes les proprits mcaniques (constantes dingnieur) et
que lon fait des suppositions.
Dterminer si la rigidit ou la souplesse de notre matriau est
dfinie positive est donc trs important et a des implications
pratiques.
Rappels
Dfinitions
Oprations
Symtries
Symtries mat.
Projecteurs.
Const. Ingnieur
Introduction
Orthotropie
Isotropie
Cubique
Iso. Transverse
Df. Positif
90 / 95
Tenseurs dfinis positifs - isotropie
Rappels
Dfinitions
Oprations
Symtries
Symtries mat.
Projecteurs.
Const. Ingnieur
Introduction
Orthotropie
Isotropie
Cubique
Iso. Transverse
Df. Positif
91 / 95
Pour un matriau isotrope, la matrice de rigidit C est donne
par:
+2
0 0 0
3
3
3
+2
0
0
0
3
3
3
2
+2
0 0 0
3
3
(103)
C= 3
0
0
0 0
0
0
0
0 0
0
0
0
0 0
Les valeurs propres de C sont obtenues en rsolvant le systme
det[C I] = 0. Compte tenu de la forme de C, on peut voir
que lon aura trois valeurs propres rptes qui valent (trois
derniers termes sur la diagonale) et les autres valeurs propres
seront donnes par:
1
det
3
+ 2
+ 2
+ 2
= ( )( )2 = 0 (104)
Tenseurs dfinis positifs - isotropie
Il faut donc que = 3k > 0 et = 2 > 0 pour que C > 0.
E
Avec les constantes dingnieur, on a que k = 3(12)
. On voit
que E > 0 et que 0.5.
E
Avec les constantes dingnieur, on a que = 2(1+)
. On voit
que E > 0 et que 1.
Donc, un tenseur isotrope sera dfini positif uniquement si
E > 0 et 1 0.5.
Rappels
Dfinitions
Oprations
Symtries
Symtries mat.
Projecteurs.
Const. Ingnieur
Introduction
Orthotropie
Isotropie
Cubique
Iso. Transverse
Df. Positif
92 / 95
Tenseurs dfinis positifs - symtrie cubique
Rappels
Dfinitions
Oprations
Symtries
Symtries mat.
Projecteurs.
Const. Ingnieur
Introduction
Orthotropie
Isotropie
Cubique
Iso. Transverse
Df. Positif
93 / 95
Pour un matriau symtrie cubique, la matrice de rigidit C
est donne par:
+2
0 0 0
3
3
3
+2
0
0
0
3
3
3
2
+2
3
0 0 0
3
3
C=
(105)
0
0
0
0 0
0
0
0
0 0
0
0
0
0 0
Les valeurs propres de cette matrice sont {, , , , , }.
Donc, pour un tenseur symtrique cubique dfini positif, on
aura:
1. E > 0
2. G > 0
3. 1 0.5
Tenseurs dfinis positifs - isotropie transverse
Rappels
Dfinitions
C=
Symtries
Symtries mat.
Projecteurs.
Const. Ingnieur
Introduction
Orthotropie
Isotropie
Cubique
Iso. Transverse
Df. Positif
Pour un matriau isotrope transverse, la matrice de rigidit C
est donne par:
Oprations
94 / 95
+
2
2
+
2
0
0
0
0
0
0
0
0
0
0 0
0 0
0
0
0 0
0
0
0
0
(106)
Comme pour lisotropie et la symtrie cubique, on peut voir que
{ , , } seront des valeurs propres. Pour les autres valeurs
propres, on doit solutionner pour
2
2
2
(107)
det
=0
2
2
2
2
2
Tenseurs dfinis positifs - isotropie transverse
Rappels
Dfinitions
Oprations
Symtries
Symtries mat.
Projecteurs.
Const. Ingnieur
Introduction
Orthotropie
Isotropie
Cubique
Iso. Transverse
Df. Positif
95 / 95
Si on fait les calculs, on aura lquation caractristique suivante:
2
2
(108)
( ) ( + ) + = 0
On voit que est une valeur propre et les deux autres seront
donnes par la solution de lquation de la parenthse de droite,
soit
q
+ ( + )2 4 ( 2 )
5 , 6 =
(109)
2
On voit donc que les deux dernires valeurs propres seront
positives si > 0, > 0 et 2 .
Alors, si on rcapitule, C > 0 si et seulement si:
> 0;
> 0;
> 0;
> 0;
(110)
En connaissant le lien entre les constantes dingnieur et les
modules donns plus haut, on peut tirer les ingalits
respecter entre les diverses constantes.