From b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b Mon Sep 17 00:00:00 2001
From: priyanka
Date: Wed, 24 Jun 2015 15:03:17 +0530
Subject: initial commit / add all books

---
 764/CH4/EX4.2.b/result4_2.txt   |  51 +++++++++++++++++++
 764/CH4/EX4.2.b/solution4_2.sce | 105 ++++++++++++++++++++++++++++++++++++++++
 2 files changed, 156 insertions(+)
 create mode 100755 764/CH4/EX4.2.b/result4_2.txt
 create mode 100755 764/CH4/EX4.2.b/solution4_2.sce

(limited to '764/CH4/EX4.2.b')

diff --git a/764/CH4/EX4.2.b/result4_2.txt b/764/CH4/EX4.2.b/result4_2.txt
new file mode 100755
index 000000000..61640e334
--- /dev/null
+++ b/764/CH4/EX4.2.b/result4_2.txt
@@ -0,0 +1,51 @@
+-->//(Design against Static Load) Example 4.2
+ 
+-->//Refer section 4.10 on page 90
+ 
+-->//Axial tensile force acting on acting on each rod P (kN)
+ 
+-->P = 50
+ P  =
+ 
+    50.  
+ 
+-->//Tensile yield strength of the material Syt (N/mm2)
+ 
+-->Syt = 400
+ Syt  =
+ 
+    400.  
+ 
+-->//Factor of safety for rod f1
+ 
+-->f1 = 6
+ f1  =
+ 
+    6.  
+ 
+-->//Factor of safety for cotter f2
+ 
+-->f2 = 4
+ f2  =
+ 
+    4.  
+ 
+
+Diameter of the rods(d) = 31.000000 mm
+
+Thickness of cotter(t) = 10.000000 mm
+
+Diameter of spigot(d2) = 40.000000 mm
+
+Outer diameter of socket(d1) = 55.000000 mm
+
+Diameter of spigot collar(d3) = 47.000000 mm
+
+Diameter of socket collar(d4) = 80.000000 mm
+
+Width of cotter(b) = 50.000000 mm
+
+Thickness of spigot collar(t1) = 15.000000 mm
+
+Design of cotter joint is safe
+ 
\ No newline at end of file
diff --git a/764/CH4/EX4.2.b/solution4_2.sce b/764/CH4/EX4.2.b/solution4_2.sce
new file mode 100755
index 000000000..0f4a52c27
--- /dev/null
+++ b/764/CH4/EX4.2.b/solution4_2.sce
@@ -0,0 +1,105 @@
+
+//Function to calculate roots of a quadratic equation
+function[r] = quadratic(a,b,c)
+    //Calculate discriminant D
+    D = (b^2)-(4 * a * c)
+    r1 = ((-1 * b)+ sqrt(D))/(2 * a)
+    r2 = ((-1 * b)- sqrt(D))/(2 * a)
+    if(r1 > 0)
+        r = r1
+    else
+        r = r2
+    end
+endfunction
+
+//Function to round-up a value such that it is divisible by 5
+function[v] = round_five(w)
+    v = ceil(w)
+    rem = pmodulo(v,5)
+    v = v + (5 - rem)
+endfunction
+
+//Obtain path of solution file
+path = get_absolute_file_path('solution4_2.sce')
+//Obtain path of data file
+datapath = path + filesep() + 'data4_2.sci'
+//Clear all
+clc
+//Execute the data file
+exec(datapath)
+//Refer Fig.4.17 on page 93
+//For rod
+//Permissible tensile stress sr1 (N/mm2)
+sr1 = Syt/f1
+//Permissible compressive stress sr2 (N/mm2)
+sr2 = 2 * sr1
+//Yield strength of the material in shear Ssy (N/mm2)
+Ssy = (50/100)*Syt  
+//Permissible shear stress tau1 (N/mm2)
+tau1 = Ssy/f1
+//For cotter
+//Permissible tensile stress sc1 (N/mm2)
+sc1 = Syt/f2
+//Permissible shear stress tau2 (N/mm2)
+tau2 = Ssy/f2
+//Calculate diameter of the rods d (mm)
+d = sqrt((4 * P * 1000)/(%pi * sr1))
+//Round up d
+d = ceil(d)
+//Calculate thickness of cotter t (mm)
+t = 0.31 * d
+//Round up t
+t = ceil(t)
+//Calculate the diameter of spigot d2 (mm)
+//a,b,c are the coefficients of the resulting quadratic equation
+//ax^2 + bx + c = 0  x=d2
+a = (%pi/4)*sr1
+b = -1* t * sr1
+c = -1 * P * 1000
+//Call the declared functions
+d2 = quadratic(a, b, c)
+d2 = round_five(d2)
+//Calculate outer daimeter of the socket d1 (mm)
+//a,b,c are the coefficients of the resulting quadratic equation
+//ax^2 + bx + c = 0  x=d1
+a = (%pi/4)*(sr1)
+b = (-1)*(t * sr1)
+c = (((-1 * (%pi/4) * (d2^2))+(t * d2))*sr1)-(P * 1000)
+//Call the declared functions
+d1 = quadratic(a, b, c)
+d1 = round_five(d1)
+//Calculate diameters of spigot collar d3 and socket collar d4 (mm)
+d3 = 1.5 * d
+d3 = ceil(d3)
+d4 = 2.4 * d
+d4 = round_five(d4)
+//Calculate dimensions a1 and c1 (mm)
+a1 = 0.75 * d
+a1 = ceil(a1)
+c1 = 0.75 * d
+c1 = ceil(c1)
+//Calculate width of cotter b (mm)
+b = (P * 1000)/(2 * tau2 * t)
+//Calculate thickness of spigot collar t1 mm
+t1 = 0.45 * d
+t1 = round_five(t1)
+//Print results
+printf('\nDiameter of the rods(d) = %f mm\n',d)
+printf('\nThickness of cotter(t) = %f mm\n',t)
+printf('\nDiameter of spigot(d2) = %f mm\n',d2)
+printf('\nOuter diameter of socket(d1) = %f mm\n',d1)
+printf('\nDiameter of spigot collar(d3) = %f mm\n',d3)
+printf('\nDiameter of socket collar(d4) = %f mm\n',d4)
+printf('\nWidth of cotter(b) = %f mm\n',b)
+printf('\nThickness of spigot collar(t1) = %f mm\n',t1)
+//Check for crushing and shear stresses in spigot end
+sigc = (P * 1000)/(t * d2)
+tauc = (P * 1000)/(2 * a1 * d2)
+//Check for crushing and shear stresses in socket end
+sigs = (P * 1000)/((d4 - d2) * t)
+taus = (P * 1000)/(2 * (d4 - d2) * c1)
+if ((sigc < sr2) & (tauc < tau1) & (sigs < sr2) & (taus < tau1))
+    printf('\nDesign of cotter joint is safe\n')
+else
+    printf('\nDesign of cotter joint is not safe\n')
+end
-- 
cgit