| 1 |
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) |
