summaryrefslogtreecommitdiff
path: root/1205/CH4/EX4.8/S_4_8.sce
blob: 7aed83b7152312dce1468818f3bd6363e2f4f684 (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
clc;
W=-1200;//N,j  Weight 
BD=[-2.4,1.2,-2.4];//m, Vector BD
EC=[-1.8,0.9,0.6];//m, Vector EC
//T_BD=norm(T_BD)*BD/norm(BD);// m, vector of tension in BD
//T_EC=norm(T_EC)*EC/norm(EC);// m, vector of tension in EC
// Applying equillibrium conditions we get
// Sum_F=0, and Sum(M_A)=0 and setting co-efficient equal to zero
A=[0.8,0.771;1.6,-0.514];//MAtrix of co-efficient
b=[-1440;0];//matrix b
x=linsolve(A,b);// solution matrix
T_BD=x(1);// N,Tension in BD
T_EC=x(2);//N, Tension in EC
printf("T_BD= (%.0f N) and T_EC= (%.0f N) \n",x(1),x(2));

Ax=2/3*T_BD+6/7*T_EC;//N, x component of reaction at A
Ay=-(1/3*T_BD+3/7*T_EC+W);//N, Y component of rection at A
Az=2/3*T_BD-2/7*T_EC;//N, z component of reaction at A

printf("Reaction at A is A=(%.0f N)i +(%.0f N)j +(%.1f N)k \n",Ax,Ay,Az);