Analyse Tensorielle en Génie Civil
Analyse Tensorielle en Génie Civil
Introduction à l’
Analyse tensorielle
D. Razafindralandy
Licence 3 Génie Civil
Université de La Rochelle
Razafindralandy (Cours en construction) 2 Analyse tensorielle (L3 Génie Civil)
Introduction
∂T
On sait que la chaleur se propage selon l’équation = κ∆T .
∂t
28°
15°
Plan du cours
Systèmes de coordonnées et métrique
Tenseurs d’ordre 1 : vecteurs et formes linéaires
Tenseurs d’ordre supérieur et d’ordre 0
Opérations tensorielles
Analyse tensorielle
Pré-requis
Algèbre linéaire : changement de base, dual, matrice, . . .
Calcul différentiel élémentaire : gradient, divergence, . . .
Doc
A la B.U. : livres d’algèbre linéaire et d’analyse tensorielle
Razafindralandy (Cours en construction) 5 Analyse tensorielle (L3 Génie Civil)
Chapitre 1
Coordonnées et
métrique
Razafindralandy (Cours en construction) 6 Analyse tensorielle (L3 Génie Civil)
Coordonnées cartésiennes
z
Base canonique de R3
M Base C = (~ex ,~ey ,~ez ) fixés, base
ez orthonormée
y
Coordonnées : x , y , z
x ey
ex
−−→
OM = x~ex + y~ey + z~ez
Plus généralement : base cartésienne
~e1 , ~e2 , ~e3 fixés
mais pas forcément orthogonaux ni normés x3 M
Coordonnées cartésiennes ou rectilignes
e3 x2
x 1, x 2, x 3
e2
On
peut écrire (matrice de passage) e1
x = p11 x 1 + p21 x 2 + p31 x 3 , x1
y = p12 x 1 + p22 x 2 + p33 x 3 ,
p13 x 1 + p23 x 2 + p33 x 3
z =
Razafindralandy (Cours en construction) 7 Analyse tensorielle (L3 Génie Civil)
Coordonnées usuelles
z
M Coordonnées
cylindriques
ez x = ρ cos θ
ey y
y = ρ sin θ
z = z
x
ex
z
Coordonnées sphériques
M
x = r cos θ sin ϕ r
ez
y = r sin θ sin ϕ
ey y
z = r cos ϕ
ϕ=colatitude. On peut aussi utiliser la latitude x
ex
Razafindralandy (Cours en construction) 8 Analyse tensorielle (L3 Génie Civil)
Base naturelle
Déplacement infinitésimale
(ex : calcul de déformation dans un milieu continu)
M'
−−→ −−→ −−→ dM
−→ −−−→ ∂ OM 1 ∂ OM 2 ∂ OM 3 M
dM = dOM = dx + dx dx
∂x 1 ∂x 2 ∂x 3
O
Base naturelle
−−→
∂ OM
~ei = , i = 1, · · · , n=2 ou 3
∂x i
−→
de telle sorte que dM = dx 1~e1 + dx 2~e2 + dx 3~e3
−→ −→ −→
−→ ∂ OM ∂ OM ∂ OM
Ex : dM = dρ~
eρ + dθ~
eθ + dϕ~
eϕ , ~
eρ = , ~
eθ = , ~
ez =
∂ρ ∂θ ∂z
Razafindralandy (Cours en construction) 9 Analyse tensorielle (L3 Génie Civil)
ey
dy dM
M ex
ey dx
ex −→
O dM = dx~ex + dy~ey + dz~ez
Razafindralandy (Cours en construction) 10 Analyse tensorielle (L3 Génie Civil)
Tenseur métrique
g11 g12 g13
Tenseur métrique (gij ) = g21 g22 g23 où gij = ~ei ~ej
g31 g32 g13
Tenseur ? On verra plus tard.
Métrique ? Permet de calculer des distances ei ||2 )
(gii = ||~ et des angles v = ||~
u ~
(~ u|| ||~
v || cos α)
Pour
la base naturelle cylindrique
~eρ ~eρ ~eρ ~eθ ~eρ ~ez 1 0 0
G = ~eθ ~eρ ~eθ ~eθ ~eθ ~ez = 0 ρ2 0
~ez ~eρ ~ez ~eθ ~ez ~ez 0 0 1
(
1 si i = j
Base orthonormée si gij =
0 6 j
si i =
Razafindralandy (Cours en construction) 11 Analyse tensorielle (L3 Génie Civil)
Notations
(
1 si i = j
Symbole de Kronecker δij = δji = δ ij =
0 6 j
si i =
Ex : δ8,8 = 1, δ23 = 0
Une base est orthonormée si gij = δij ∀i, j
Ex : δij
= symbole de
Kronecker
1 0 0
(δij ) = 0 1 0
0 0 1
Razafindralandy (Cours en construction) 12 Analyse tensorielle (L3 Génie Civil)
Remarques
Vecteur
déplacement
ux , uy , uz
poids de la croix+structure
6= u r , u θ , u ϕ
cartésien
ez
M ey M 6= u 1 , u 2 , u 3
ex
u 6= · · ·
u
u=uxex+uyey+uzez Mais c’est le même
vecteur ~u ⇒ Il y a
forme finale?
e1
cylindrique
sphérique
u e u e3 u
e
r z
u=u er+u e +u e u=u e +u e +u ez u=u e1+u2e2+u3e3
1
Razafindralandy (Cours en construction) 14 Analyse tensorielle (L3 Génie Civil)
Cas général
Deux bases B = (~e1 ,~e2 ,~e3 ) et B 0 = (~e10 ,~e20 ,~e30 ) avec
Matrice de passage de B à B 0
~e 0 ~e20 ~e30
11
p21 p31 → sur →
−
p1 e1
[pji ] = p12 p22 p32
→ sur −
→
e2
p13 p23 p33 → sur →
−
e3
u 0j = qij u i N • = QA• Q = P −1
P
i ou
Les composantes d’un vecteur sont dites contravariantes N • = P −1 A•
Razafindralandy (Cours en construction) 16 Analyse tensorielle (L3 Génie Civil)
Exercice
Chapitre 2
Tenseurs
d’ordre 1
Razafindralandy (Cours en construction) 18 Analyse tensorielle (L3 Génie Civil)
u1 v1 X3
~u ~v = u 2 v 2 = u 1 v 1 + u 2 v 2 + u 3 v 3 = ui v i
u3 v3 i=1
Rappel :
g11 g12 g13 Convention :
(gij ) = g21 g22 g23 .
1e indice=num. ligne,
g31 g32 g33 2e indice=num. colonne
IIIExiste-t’il
X des quantités u1 , u2 , u3 telles que
~u ~v = uj v j comme dans le cas orthonormé ?
j
3 X
3 3 3
!
X X X
~u ~v = gij u i v j = gij u i vj
i=1 j=1 j i=1
| {z }
uj
n
X n
X
~u ~v = uj v j où uj = gij u i
j=1 i=1
n n
! rapport au vecteur ~u ?
IIIQue représentent les nombres uj par
X X X
uj = u i gij = u i (~ei ~ej ) = u i ~ei ~ej , ⇒
i=1 i=1 i
e2 e2
uj = ~u ~ej
u u
u2
u2
u i = Coordonnées dans la direction ~ei
u1 e1 u1 e1
ui = Mesure de ~u vue par ~ei
||~ei || = 1
Razafindralandy (Cours en construction) 21 Analyse tensorielle (L3 Génie Civil)
X
Mais : Tpp =?
p
Razafindralandy (Cours en construction) 22 Analyse tensorielle (L3 Génie Civil)
Espace dual
E ∗ = { formes linéaires E → R}
Matrice d’une application linéaire
x
x x x x
M : R3 → R 2 M~u = x =
x x x x
x
x
f : R3 → R f (~u ) = x x x x = x
x
Matrice d’une forme
1 linéaire f : (a b c)
u
f (~u ) = (a b c) u 2 = au 1 + bu 2 + cu 3
u3
Ex : un déplacemet se mesure directement
Une force peut
se “mesurer” indirectement à partir de son travail
ux
F : ~v = u y 7→ F (~u ) = W = Fx u x + Fy u y + Fz u z
uz
Force : forme linéaire de matrice (Fx Fy Fz )
Razafindralandy (Cours en construction) 23 Analyse tensorielle (L3 Génie Civil)
Base duale
f (~u ) = au 1 + bu 2 + cu 3
Les formes linéaires ~εi forment une base de E ∗ , appelée base duale
La base duale vérifie :
~εi ~ej = δji
Tenseur d’ordre 1
Changement de base
Vecteur
déplacement
ux , uy , uz
poids de la croix+structure
6= u r , u θ , u ϕ
cartésien
ez
M ey M 6= u 1 , u 2 , u 3
ex
u 6= · · ·
u
u=uxex+uyey+uzez Mais c’est le même
vecteur ~u ⇒ Il y a
forme finale?
e1
cylindrique
sphérique
u e u e3 u
e
r z
u=u er+u e +u e u=u e +u e +u ez u=u e1+u2e2+u3e3
1
Razafindralandy (Cours en construction) 27 Analyse tensorielle (L3 Génie Civil)
Matrice de passage
~e10 ~e20 ~e30
→ sur →
−
1
p21 p31
p1 e1
[pji ] = p12 p22 p32 → sur −
→
e2
p13 p23 p33 → sur →
−
e3
Razafindralandy (Cours en construction) 28 Analyse tensorielle (L3 Génie Civil)
De manière matricielle
Règle de changement de base pour un vecteur
A• = PN • ou N • = QA•
Razafindralandy (Cours en construction) 29 Analyse tensorielle (L3 Génie Civil)
A• = QN• ou N• = PA•
En résumé,
Règles de changement de base
u i = pji u 0j ou u 0j = qij u i
ui = qij u 0j ou uj0 = pji ui
Définition physique d’un tenseur d’ordre 1
Un tenseur d’ordre 1 est une suite de nombres, à un indice, vérifiant
une de ces règles de changement de base
u• =? u• , u• = g u• , ui = gi u• , ui = gij uj
Chapitre 3
Tenseurs
d’ordre supérieur
Razafindralandy (Cours en construction) 32 Analyse tensorielle (L3 Génie Civil)
Tenseurs d’ordre 2
1 1
u v
~u = ~
,v=
u2 v2
= u 1 v 1 + 2u 1 v 2 + 0 + 3u 2 v 2
Razafindralandy (Cours en construction) 33 Analyse tensorielle (L3 Génie Civil)
Changement de variance
F t •• =???t ••
F t •• =???t••
Razafindralandy (Cours en construction) 35 Analyse tensorielle (L3 Génie Civil)
Changement de base
Définition physique
Un tenseur d’ordre 2 est une suite/un tableau de nombres à 2 indices qui
obéit à une de ces règles de changement de base
Le tableau [gij ]
Voyons comment il change de base
gkl0 = ~ek0 ~el0 = (pki ~ei ) (plj~ej ) = pki plj gij
Les g kl , définis par (g kl ) = (gkl )−1 , n’en seraient pas les composantes
contravariantes ?
Notons hij les composantes 2 fois covariantes de G.
h•• =??g•• hk• = g ki gi• , hkl = g ki g lj gij
Or g est symétrique, donc h = g gij g = δjk g jl = g kl .
lj kl ki jl
g kl = δlk et gl k = δlk .
Un contre-exemple de tenseur.
Une matrice de passage ne forme pas un tenseur
Ce n’est pas la peine de décaler les indices : pji
Razafindralandy (Cours en construction) 40 Analyse tensorielle (L3 Génie Civil)
Tenseur d’ordre 0.
Changement de variance : t aucun
Changement de base : r = t
Exemple
T''=18oC
d d''=2.2m
M
cylindrique
e
T=18oC ez
cartésien
d=2.2m M e
T = T 0 = T 00
ez e d = d 0 = d 00
M ey
er
ex M T'=18oC
d'=2.2m
e e sphérique
Chapitre 4
Opérations
algébriques
Razafindralandy (Cours en construction) 43 Analyse tensorielle (L3 Génie Civil)
Dans une autre base, t 0ij = u 0i v 0j = (qki u k )(qlj v l ) =⇒ t 0ij = qki qlj t ij
ij
Donc, les t sont les composantes 2 fois contravariantes d’un tenseur d’ordre 2,
noté ~u ⊗ ~v .
Ex : Dans Navier-Stokes incompressible : (~u ∇)~u = div(u ⊗ u)
1 −1
Exemple. ~u = , ~v = . Composantes de t = ~u ⊗ ~v et r = ~v ⊗ ~u ?
1 0
t 11 = u 1 v 1 = − 1, t 12 = u 1 v 2 = 0, t 21 = u 2 v 1 = − 1, t 22 = u 2 v 2 = 0
r 11 =
Contraction
Contraction = Egalisation de 2 indices de variances différentes (et
sommation)
Ex : T = (t ij k l ). Contraction(T ) selon le 1er et le 3è indices = t ij i l
Ex. Contraction(t ij ) = t ii = t 11 + t 22 + t 33 = tr t
Produit contracté
T ⊗S =Contraction(T ⊗ S) T et S d’ordres > 0
~u ⊗~v = contraction(ui v j ) = ui v i != ~u ~v
v1
=Produit matriciel (u1 u2 ) 2
v
Matrice⊗vecteur : (mi j )⊗(v k ) = contraction(mi j v k ) = mi j v j
=Produit matrice-vecteur
Matrice⊗Matrice : (mi j )⊗(s kl ) = contraction(mi j s kl ) = mi j s jl
=Produit matrice-matrice
Trace et déterminant
Trace. tr T = contraction(T ) = T ii = Ti i = g ij Tij = gij T ij .
Transposition
> >
T est la forme bi-linéaire (donc tenseur) tq T (~u , ~v ) = T (~v , ~u )
(>
T )ij = T ji , >
( T )ij = Tji , (>
T )i j = Tj i, (> T )i j = T j i
>
Symétrique si T =T Antisymétrique si > T = −T
Décomposition en parties symétrique et antisymétrique
T = sym(T ) + antisym(T )
T +> T T −> T
sym(T ) = , antisym(T ) =
2 2
Déviateur si tr T = 0
Décomposition en parties sphérique et déviatorique
1
T = ( tr T )G + dev (T ) n = dim E
n
Razafindralandy (Cours en construction) 49 Analyse tensorielle (L3 Génie Civil)
Chapitre 5
Analyse tensorielle
Razafindralandy (Cours en construction) 51 Analyse tensorielle (L3 Génie Civil)
∂v i
Soit ~v = v i ~ei un champs de vecteurs. Est-ce que les forment un tenseur ?
∂x j
Changement de base, de passage P. x i = pji x 0j
∂v i ∂v i ∂x 0k ∂pli v 0l k i k
∂v 0l 0l
∂pli k
= = q j = p qj + v q .
∂x j ∂x 0k ∂x j ∂x 0k l
∂x 0k ∂x 0k j
A cause du dernier terme, qui ne s’annule pas lorsque la base n’est pas rectiligne,
∂v i
les ne forment pas un tenseur.
∂x j
(
~er = cos θ~ex + sin θ~ey
Ex.
~eθ = −r sin θ~ex + r cos θ~ey
∂~e1
= Γ1j 1~ej = Γ111~e1 + Γ121~e2
∂x 1
∂~e1 ∂~er →
−
1
= = 0 = 0~e1 + 0~e2 ⇒ Γ111 = Γ121 = 0
∂x ∂r
Razafindralandy (Cours en construction) 54 Analyse tensorielle (L3 Génie Civil)
Et Γj i k = g il Γjlk
Symétrie : Γi j k = Γkj i , Γijk = Γkji
−→ −→
∂ 2 OM ∂ 2 OM
(se déduit du fait que = .)
∂x i ∂x j ∂x j ∂x i
Ce ne sont pas des tenseurs ne vérifient pas les formules de changement de base
∂~εj
= −Γi j k ~εk .
∂x i
∂~ej →
−
Dans une base rectiligne, k
= 0 , ∀j, k. Donc Γi j k = 0 et Γjik = 0
∂x
Razafindralandy (Cours en construction) 55 Analyse tensorielle (L3 Génie Civil)
d~v δvi δv i dv i dx l
= ~ei où = + v k Γik l .
dt δt δt dt dt
δv i
!
est un tenseur d’ordre 1, appelé dérivée covariante de ~v par
δt
rapport à t. Les composantes covariantes vérifient
δvi dvi dx l
= − vk Γki l .
δt dt dt
Razafindralandy (Cours en construction) 56 Analyse tensorielle (L3 Génie Civil)
δv i dv i dx l δvi dvi dx l
= + v k Γik l = − vk Γki l
δt dt dt δt dt dt
δT ijk dT ijk dx m dx m dx m
Ex : = + Γl i m T ljk − Γj l m T ilk − Γml k T ijl
δt dt dt dt dt
δf ∂f
Pour un scalaire : = .
δt ∂t
Razafindralandy (Cours en construction) 57 Analyse tensorielle (L3 Génie Civil)
δv i ∂v i dx l ∂v i
Lorsque t = x j , = + v k Γki l = + v k Γki l δjl
δx j ∂x j dx j ∂x j
δv i
!
Gradient d’un vecteur : ( grad ~v ij ) =
δx j
δv i ∂v i δvi ∂vi
= + v k Γki j et = − vk Γi kj
δx j ∂x j δx j ∂x j
Divergence d’un vecteur
√
δv i 1 ∂ vi g
div ~v = tr( grad ~v ) = i et div ~v = √
δx g ∂x i
Gradient et divergence d’un tenseur d’ordre 2
δT ij
!
−−→ δT ij ∂T ij
( gradT )ij k = où = + T mj Γki m + T im Γkj m
δx k δx k ∂x k
−→ δT ij −−→
div T = ~ei (contraction des 2 derniers indices de gradT )
δx j
Se généralise à un tenseur d’ordre quelconque
Razafindralandy (Cours en construction) 58 Analyse tensorielle (L3 Génie Civil)
Chapitre 6
Quelques
applications
Razafindralandy (Cours en construction) 60 Analyse tensorielle (L3 Génie Civil)
Mécanique du point
S1
S2
Robotique
Coordonnées : 3 angles + 1 distance S3
Razafindralandy (Cours en construction) 62 Analyse tensorielle (L3 Génie Civil)
Etude de surface/coque
1
x = cos θ
x = cos θ cos ϕ
x = f (r , s)
y = sin θ y = sin θ cos ϕ y = f 2 (r , s)
z =z z = sin ϕ z = f 3 (r , s)
La surface n’est pas un espace vectoriel, mais son plan tangent en est un. On peut y appliquer
l’analyse tensorielle.
∂f 1 ∂f 2 ∂f 3
~er = ~ex + ~ey + ~ez
∂r ∂r ∂r ~er · ~er ~er · ~es
∂f 1
∂f 2
∂f 3 g=
~es = ~ex + ~ey + ~ez ~es · ~er ~es · ~es
∂s ∂s ∂s
(base de l’espace tangent)
E = λG⊗G + µI où Iijkl = g ik g jl + g il g jk
σ = λ( tr ε)G + 2µε
Fonction invariante, . . .
Razafindralandy (Cours en construction) 65 Analyse tensorielle (L3 Génie Civil)
Formules (compliquées)