initial commit / add all books

4 files changed, 96 insertions, 0 deletions

diff --git a/3705/CH1/EX1.1/Ex1_1.sce b/3705/CH1/EX1.1/Ex1_1.sce
new file mode 100644
index 000000000..d86a5575b
--- /dev/null
+++ b/3705/CH1/EX1.1/Ex1_1.sce
@@ -0,0 +1,22 @@
+
+clear//

+

+//NOTE:The notation has been changed to simplify the coding process

+

+//Variable Declaration

+P_AB=4000 //Axial Force at section 1 in lb

+P_BC=5000 //Axial Force at section 2 in lb

+P_CD=7000 //Axial Force at section 3 in lb

+A_1=1.2 //Area at section 1 in in^2

+A_2=1.8 //Area at section 2 in in^2

+A_3=1.6 //Area at section 3 in in^2

+

+//Calculation

+//S indicates sigma here

+S_AB=P_AB/A_1 //Stress at section 1 in psi (T)

+S_BC=P_BC/A_2 //Stress at section 2 in psi (C)

+S_CD=P_CD/A_3 //Stress at section 3 in psi (C)

+

+//Result

+printf("\n The stress at the three sections is given as")

+printf("\n Stress at section 1= %0.0f psi/in^2 section 2= %0.0f psi/in^2 section 3= %0.3f psi/in^2",S_AB,S_BC,S_CD)

diff --git a/3705/CH1/EX1.2/Ex1_2.sce b/3705/CH1/EX1.2/Ex1_2.sce
new file mode 100644
index 000000000..f083daa4e
--- /dev/null
+++ b/3705/CH1/EX1.2/Ex1_2.sce
@@ -0,0 +1,29 @@
+
+clear//

+

+//Variable Declaration

+Ay=40 //Vertical Reaction at A in kN

+Hy=60 //Vertical Reaction at H in kN

+Hx=0 //Horizontal Reaction at H in kN

+y=3 //Height in m

+x=5 //Distance in m

+p=4 //Panel distance in m

+A=900 //Area of the member in mm^2

+P_C=30 //Force at point C in kN

+

+//Calculation

+//Part 1

+//Applying summation of forces in the x and y direction and equating to zero

+P_AB=(-Ay)*(x*y**-1) //Force in member AB in kN

+P_AC=-(p*x**-1*P_AB) //Force in member AC in kN

+//Using stress=force/area

+S_AC=(P_AC/A)*10**3 //Stress in member AC in MPa (T)

+

+//Part 2

+//Sum of moments about point E to zero

+P_BD=(Ay*p*2-(P_C*p))*y**-1 //Force in memeber AB in kN (C)

+S_BD=(P_BD/A)*10**3 //Stress in member in MPa (C)

+

+//Result

+printf("\n The Stress in member AC is %0.1f MPa (T)",S_AC)

+printf("\n The Stress in member BD is %0.1f MPa (C)",S_BD)

diff --git a/3705/CH1/EX1.3/Ex1_3.sce b/3705/CH1/EX1.3/Ex1_3.sce
new file mode 100644
index 000000000..e3daf96df
--- /dev/null
+++ b/3705/CH1/EX1.3/Ex1_3.sce
@@ -0,0 +1,28 @@
+
+clear//

+//

+

+//Variable Declaration

+A_AB=800 //Area of member AB in m^2

+A_AC=400 //Area of member AC in m^2

+W_AB=110 //Safe value of stress in Pa for AB

+W_AC=120 //Safe value of stress in Pa for AC

+theta1=60*3.14*180**-1 //Angle in radians

+theta2=40*3.14*180**-1 //Angle in radians 

+

+//Calculations

+//Applying sum of forces 

+//Solving by matrix method putting W as 1

+A =[-cos(theta1),cos(theta2);sin(theta1),sin(theta2)]

+

+B = [1;1]

+C=inv(A)

+D=C

+

+//Using newtons third law

+//Two values of W hence the change in the notation

+W1=(W_AB*A_AB)*D(2,2)**-1 //Weight W in N

+W2=(W_AC*A_AC)*D(1,2)**-1 //Weight W in N

+

+//Result

+printf("\n The maximum value of W allowable is %0.1f kN",W2*1000**-1)

diff --git a/3705/CH1/EX1.4/Ex1_4.sce b/3705/CH1/EX1.4/Ex1_4.sce
new file mode 100644
index 000000000..3b1200aad
--- /dev/null
+++ b/3705/CH1/EX1.4/Ex1_4.sce
@@ -0,0 +1,17 @@
+
+clear//

+

+//Variable Declaration

+d=3*4**-1 //Rivet diameter in inches

+t=7*8**-1 //Thickness of the plate in inches

+tau=14000 //Shear stress limit in psi

+sigma_b=18000 //Normal stress limit in psi

+

+//Calculations

+//Design Shear Stress in Rivets

+V=tau*(d**2*(%pi/4))*4 //Shear force maximum allowable in lb

+//Design for bearing stress in plate

+Pb=sigma_b*t*d*4 //lb

+

+//Result

+printf("\n The maximum load that the joint can carry is %0.0f lb",V)