Résolution de l'équation de Navier-Lamé

Nous arrivons enfin dans le vif du sujet. Pour arriver à calculer notre déformation élastique,
nous allons résoudre l'équation de Navier-Lamé,ci-dessous, en 3D.

2 u + 1 1 2 ν . u = 0 nabla ^2 vec u + {1} over {1-2%nu} nabla nabla . vec u =0

Nous allons prendre une déformation u de la forme suivante :

u = e iKx X e iKyY A ( z ) vec u = func e^{iKx X}func e^{iKyY} A(z)

Nous voulons donc déterminer le vecteur A(z), qui nous donnera la forme de la déformation.

A ( z ) = A1 ( z ) A2 ( z ) A3 ( z ) vec A(z)= left ( stack{A1(z) # A2(z) # A3(z)} right )

Passons les longs calculs (mais qui sont détaillés ICI pour les intéréssés). Voici le résultat obtenu :

Kx 2 Ky 2 Kx 2 1 2 ν KxKy 1 2 ν iKx 1 2 ν λ Kx 2 KxKy 1 2 ν Ky 2 iKy 1 2 ν λ iKx 1 2 ν λ iKy 1 2 ν λ Kx 2 Ky 2 + λ 2 A1 A2 A3 = 0 0 0 font sans bold size *1 left[matrix{-Kx^2-Ky^2-{Kx^2} over {1-2%nu}#{-KxKy} over {1-2%nu}#{iKx} over {1-2%nu} %lambda##-Kx^2-{KxKy} over {1-2%nu}#-Ky^2# {iKy} over {1-2%nu}%lambda ###{iKx} over {1-2%nu}%lambda#{iKy} over {1-2%nu}%lambda#-Kx^2-Ky^2+%lambda^2}right]*font sans bold size*1 left[stack{A1 # A2 # A3}right]= font sans bold size *1 left[stack{0 # 0 # 0}right]