initial commit / add all books

2 files changed, 156 insertions, 0 deletions

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