| 1 |
francois |
941 |
clc;
|
| 2 |
|
|
disp('Solution exercice 3 du chapitre introduction');
|
| 3 |
|
|
disp('1 Resolution du probleme avec charge fictive P en B');
|
| 4 |
|
|
syms Ra Rc w L t P x E I v Ma teta;
|
| 5 |
|
|
eq1=Ra+Rc+w*L+P;
|
| 6 |
|
|
eq2=Rc*L+int(w*t,t,0,L)+P*L/2;
|
| 7 |
|
|
R=solve(eq1,eq2,Ra,Rc);
|
| 8 |
|
|
disp('Ra=');
|
| 9 |
|
|
R.Ra
|
| 10 |
|
|
disp('Rc=');
|
| 11 |
|
|
R.Rc
|
| 12 |
|
|
Mab=R.Ra*x+int(w*(x-t),t,0,x);
|
| 13 |
|
|
Mab
|
| 14 |
|
|
Mbc=Mab+P*(x-L/2);
|
| 15 |
|
|
Mbc
|
| 16 |
|
|
v=int(Mab*diff(Mab,P)/E/I,x,0,L/2)+int(Mbc*diff(Mbc,P)/E/I,x,L/2,L);
|
| 17 |
|
|
disp('vb=');
|
| 18 |
|
|
limit(v,P,0)
|
| 19 |
|
|
disp('2 Resolution du probleme statique avec moment fictif en A');
|
| 20 |
|
|
eq3=Ra+Rc+w*L;
|
| 21 |
|
|
eq4=Rc*L+int(w*t,t,0,L)+Ma;
|
| 22 |
|
|
R=solve(eq3,eq4,Ra,Rc);
|
| 23 |
|
|
disp('Ra=');
|
| 24 |
|
|
R.Ra
|
| 25 |
|
|
disp('Rc=');
|
| 26 |
|
|
R.Rc
|
| 27 |
|
|
M=-Ma+R.Ra*x+int(w*(x-t),t,0,x);
|
| 28 |
|
|
M
|
| 29 |
|
|
teta=int(M*diff(M,Ma)/E/I,x,0,L);
|
| 30 |
|
|
disp('tetea=');
|
| 31 |
|
|
limit(teta,Ma,0) |
