Scilab Code
Scilab Code
PROGRAM: RESULT:
disp(a) 2. 2. 9. 7.
0. 1. 5. - 2.
disp('R2=R2-2*R1')
R2=R2-2*R1
a(2,:)=a(2,:)-2*a(1,:);
1. 0. - 3. 8.
disp(a)
0. 2. 15. - 9.
disp('interchange R2 and R3')
0. 1. 5. - 2.
a([2,3],:)=a([3,2],:);
interchange R2 and R3
disp(a)
1. 0. - 3. 8.
disp('R3=R3-2*R2') 0. 1. 5. - 2.
a(3,:)=a(3,:)-2*a(2,:); 0. 2. 15. - 9.
disp(a) R3=R3-2*R2
disp('R3=(1/5)*R3') 1. 0. - 3. 8.
a(3,:)=(1/5)*a(3,:); 0. 1. 5. - 2.
disp(a) 0. 0. 5. - 5.
R3=(1/5)*R3
disp('R2=R2-5*R3 and R1=R1+3*R3')
1. 0. - 3. 8.
a(2,:)=a(2,:)-5*a(3,:);
0. 1. 5. - 2.
a(1,:)=a(1,:)+3*a(3,:);
0. 0. 1. - 1.
disp(a)
R2=R2-5*R3 and R1=R1+3*R3
s=[a(1,4);a(2,4);a(3,4)];
1. 0. 0. 5.
disp('solution is')
0. 1. 0. 3.
disp(s) 0. 0. 1. - 1.
solution is
5.
3.
- 1.
PROGRAM: RESULT:
// Calculation
// (a)
delta_st=delta_st*12;// in
// (b)
// (c)
M_C=(-P*(L_AB+L_BC))+(C(1)*L_BC)-
DOWNLOAD BOOK
(UDL*L_BC*(L_BC/2)) // kNm
- Chapter 25 : Shear Force and
Bending Moment M_D=0 // kNm
Example 25.7 // LOCATION OF MAXIMUM BM
// Initilization of variables // Max BM occurs at E at a distance of 2.5 m
from B i.e x=L_AE=4.5 m from free end A.
L_AD=8 // m // length of the beam
Thus max BM is given by taking moment at B
L_AB=2 // m
L_AE=4.5 // m // given
L_BC=4 // m
M_E=(-2*L_AE)+(4.5*(L_AE-2))-((1/2)*(L_AE-
L_CD=2 // m 2)^2) // kNm
// Caalculations y=[-2;-2;2.5;0;-1.5;2;0]
// REACTIONS a=[0;2;4.5;6;8]
F_D=0
// BENDING MOMENT
Fdb=((Fad*sind(theta))-(W1))/sind(theta) //N
// (T) // Using eq'n 6
Chapter 9 : Analysis of Plane
Trusses and Frames Fde=(Fdb*cosd(theta))+(Fad*cosd(theta)) //N
// (C) // Using eq'n 5
Example 9.1
// (4) JOINT E
Feb=((Fce*cosd(theta))-(Fde))/cosd(theta) //N
// Initilization of variables // (C) // Using eq'n 7
W1=2000 //N // load at joint D of the truss // Results
W2=4000 //N // load at joint E of the truss clc
Lac=6 //m // length of the tie printf('The Axial Force in member AD (Fad) is
%f N \n',Fad)
Lab=3 //m
printf('The Axial Force in member AB (Fab) is
Lbc=3 //m %f N \n',Fab)
theta=60 //degree // interior angles of the printf('The Axial Force in member CE (Fce) is
truss %f N \n',Fce)
// Calculations printf('The Axial Force in member CB (Fcb) is
// Here A is simply supported & B is roller %f N \n',Fcb)
support. Now the SUPPORT REACTIONS are printf('The Axial Force in member DB (Fdb) is
given as, %f N \n',Fdb)
Rc=((W1*(Lab/2))+(W2*(Lab+(Lbc/2))))/Lac printf('The Axial Force in member DE (Fde)is
//N // Taking moment at A %f N \n',Fde)
Ra=W1+W2-Rc //N // Take sum Fy=0 printf('The Axial Force in member EB (Feb) is
// ANALYSIS OF TRUSS BY METHOD OF JOINT %f N \n',Feb)
// (1) JOINT A
// (2) JOINT C
// (3) JOINT D