| 1 |
francois |
941 |
clc;
|
| 2 |
|
|
disp('Solution exercice 2 du chapitre deformation');
|
| 3 |
|
|
disp('tenseur des deformations en um/m');
|
| 4 |
|
|
epsilon=[400 240/2 480/2;240/2 -200 -240/2;480/2 -240/2 160]
|
| 5 |
|
|
disp('Vecteurs et valeurs propes en um/m');
|
| 6 |
|
|
[V,D]=eig(epsilon)
|
| 7 |
|
|
disp('Cisaillement maximal en um/m');
|
| 8 |
|
|
gammamax=D(3,3)-D(1,1)
|
| 9 |
|
|
disp('Direction du cisaillement maximal');
|
| 10 |
|
|
direction=[0.5*(V(1,1)+V(1,3)) 0.5*(V(2,1)+V(2,3)) 0.5*(V(3,1)+V(3,3) )]
|
| 11 |
|
|
disp('Normal plan du cisaillement max');
|
| 12 |
|
|
n=[V(1,2) V(2,2) V(2,3)]
|
