blob: 0dcb9281c750598be06c0bd25a652d2e077b1d70 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
|
pathname=get_absolute_file_path('5_1.sce')
filename=pathname+filesep()+'5_1data.sci'
exec(filename)
L=[(h^2 +AB^2)^0.5;BC;(h^2 +CD^2)^0.5;CD;BC;AB;h;(h^2 +BC^2)^0.5;h];
A=atan(h/AB),B=atan(h/BC),D=atan(h/CD);
lFC=Load2;
P=[
0 0 sin(D) 0 0 0 0 sin(B);
cos(A) 0 0 0 0 1 0 0;
cos(A) -1 0 0 0 0 0 0;
sin(A) 0 0 0 0 0 1 0;
0 0 0 -1 1 0 0 0;
0 0 0 0 1 -1 0 cos(B);
0 0 0 0 0 0 1 -sin(B);
0 1 -cos(D) 0 0 0 0 cos(B);];
X=[-Load2;0;0;-Load1;0;0;Load1;0];
D1=[0;0;0;1;1;1;0;0;0];
FC=0;
P1=[cos(A) 0 0 0 0 1 0 0;
0 0 cos(D) 1 0 0 0 0;
cos(A) -1 0 0 0 0 0 0;
sin(A) 0 0 0 0 0 1 0;
0 0 0 -1 1 0 0 0;
0 0 0 0 1 -1 0 cos(B);
0 0 0 0 0 0 1 sin(B);
0 0 1 0 0 0 0 1]
X1=[0;0;0;0;0;0;1;0]
FA=[inv(P)*X;lFC];
B1=[inv(P1)*X1;FC];
for i=1:9
K1(i)=FA(i)*B1(i)*L(i);
K2(i)=FA(i)*D1(i)*L(i);
end
delB=(sum(K1)*10^3)/(E*CS);
delD=(sum(K2)*10^3)/(E*CS);
printf("\nδB,v: %f mm",delB);
printf("\nδD,h: %f mm",delD);
|
