| 1 |
francois |
941 |
clc;
|
| 2 |
|
|
syms x Ra Rc Mc F L A E I;
|
| 3 |
|
|
Mab=Ra*x;
|
| 4 |
|
|
Mdb=-F*x;
|
| 5 |
|
|
Mcb=-Mc;
|
| 6 |
|
|
Pab=0;
|
| 7 |
|
|
Pdb=0;
|
| 8 |
|
|
Pcb=Rc;
|
| 9 |
|
|
d1Mab=x;
|
| 10 |
|
|
d1Mdb=0;
|
| 11 |
|
|
d1Mcb=-L;
|
| 12 |
|
|
d1Pab=0;
|
| 13 |
|
|
d1Pdb=0;
|
| 14 |
|
|
d1Pcb=-1;
|
| 15 |
|
|
dU=d1Pab*Pab*L/A/E+d1Pdb*Pdb*L/A/E+d1Pcb*Pcb*L/A/E+int(Mab*d1Mab/E/I,x,0,L)+int(Mdb*d1Mdb/E/I,x,0,L)+int(Mcb*d1Mcb/E/I,x,0,L)
|
| 16 |
|
|
eq1=Ra+Rc-F;
|
| 17 |
|
|
eq2=Rc*L+Mc-2*F*L;
|
| 18 |
|
|
R=solve(eq1,eq2,dU,Ra,Rc,Mc);
|
| 19 |
|
|
disp('Ra=');
|
| 20 |
|
|
simplify(R.Ra)
|
| 21 |
|
|
disp('Rc=');
|
| 22 |
|
|
simplify(R.Rc)
|
| 23 |
|
|
disp('Mc=');
|
| 24 |
|
|
simplify(R.Mc)
|
| 25 |
|
|
d2Mab=-x;
|
| 26 |
|
|
d2Mdb=-x;
|
| 27 |
|
|
d2Mcb=0;
|
| 28 |
|
|
d2Pab=0;
|
| 29 |
|
|
d2Pdb=0;
|
| 30 |
|
|
d2Pcb=2;
|
| 31 |
|
|
dU2=d2Pab*Pab*L/A/E+d2Pdb*Pdb*L/A/E+d2Pcb*Pcb*L/A/E+int(Mab*d2Mab/E/I,x,0,L)+int(Mdb*d2Mdb/E/I,x,0,L)+int(Mcb*d2Mcb/E/I,x,0,L) |
