| 1 |
francois |
941 |
clc;
|
| 2 |
|
|
disp('Solution exercice 3 du chapitre deformation');
|
| 3 |
|
|
syms epsx epsy gammaxy
|
| 4 |
|
|
eq1=(epsx+epsy)/2+(epsx-epsy)/2*cos(2*0)+gammaxy/2*sin(2*0)-1050e-6;
|
| 5 |
|
|
eq2=(epsx+epsy)/2+(epsx-epsy)/2*cos(2*pi/4)+gammaxy/2*sin(2*pi/4)-800e-6;
|
| 6 |
|
|
eq3=(epsx+epsy)/2+(epsx-epsy)/2*cos(2*pi/2)+gammaxy/2*sin(2*pi/2)+750e-6;
|
| 7 |
|
|
R=solve(eq1,eq2,eq3);
|
| 8 |
|
|
disp('epsx');
|
| 9 |
|
|
vpa(R.epsx)
|
| 10 |
|
|
disp('epsy');
|
| 11 |
|
|
vpa(R.epsy)
|
| 12 |
|
|
disp('gammaxy');
|
| 13 |
|
|
vpa(R.gammaxy)
|
| 14 |
|
|
epsilon=[R.epsx R.gammaxy/2;R.gammaxy/2 R.epsy;];
|
| 15 |
|
|
disp('tenseur des deformations');
|
| 16 |
|
|
vpa(epsilon)
|
| 17 |
|
|
[V,D]=eig(epsilon);
|
| 18 |
|
|
norm1=V(1,1)*V(1,1)+V(2,1)*V(2,1);
|
| 19 |
|
|
norm1=sqrt(norm1);
|
| 20 |
|
|
norm2=V(1,2)*V(1,2)+V(2,2)*V(2,2);
|
| 21 |
|
|
norm2=sqrt(norm2);
|
| 22 |
|
|
direction=[V(1,1)/norm1 V(1,2)/norm2;V(2,1)/norm1 V(2,2)/norm2;];
|
| 23 |
|
|
disp('Directions principales du tenseur des deformations');
|
| 24 |
|
|
vpa(direction)
|
| 25 |
|
|
disp('Deformations principales');
|
| 26 |
|
|
vpa(D) |
