| 1 |
francois |
941 |
clc;
|
| 2 |
|
|
disp('Solution exercice 2 du chapitre contraintes');
|
| 3 |
|
|
disp('tenseur des contraintes en MPa');
|
| 4 |
|
|
sigma=[9 -30 0;-30 0 90;0 90 -81]
|
| 5 |
|
|
disp('Vecteurs et valeurs propes en MPa');
|
| 6 |
|
|
[V,D]=eig(sigma)
|
| 7 |
|
|
disp('Vecteurs contraintes de la facette ABC en MPa');
|
| 8 |
|
|
A=[1; 0; 0];
|
| 9 |
|
|
B=[0; 2; 0];
|
| 10 |
|
|
C=[0; 0; -3];
|
| 11 |
|
|
AB=B-A;
|
| 12 |
|
|
AC=C-A;
|
| 13 |
|
|
n=cross(AB,AC);
|
| 14 |
|
|
n=n/norm(n,2);
|
| 15 |
|
|
vecteur_sigma=sigma*n
|
| 16 |
|
|
disp('Angle entre le vecteur contrainte et la normale a la face en degre');
|
| 17 |
|
|
beta=acos(dot(n,vecteur_sigma)/norm(n,2)/norm(vecteur_sigma,2));
|
| 18 |
|
|
beta=beta/pi*180
|
| 19 |
|
|
disp('Contrainte normale en MPa');
|
| 20 |
|
|
sigman=dot(vecteur_sigma,n)
|
| 21 |
|
|
disp('Contrainte de cisaillement en MPa');
|
| 22 |
|
|
tau=sqrt(norm(vecteur_sigma,2)^2-sigman^2)
|
