1 files changed, 36 insertions, 0 deletions
diff --git a/764/CH4/EX4.13.b/solution4_13.sce b/764/CH4/EX4.13.b/solution4_13.sce
new file mode 100755
index 000000000..72a4a0fe2
--- /dev/null
+++ b/764/CH4/EX4.13.b/solution4_13.sce
@@ -0,0 +1,36 @@
+
+//Obtain path of solution file
+path = get_absolute_file_path('solution4_13.sce')
+//Obtain path of data file
+datapath = path + filesep() + 'data4_13.sci'
+//Clear all
+clc
+//Execute the data file
+exec(datapath)
+//Calculate the yield strength in shear Ssy (N/mm2)
+Ssy = (50/100)*Syt
+//Calculate the permissible shear stress tau (N/mm2)
+tau = Ssy/fs
+//Calculate the bending moment Mb (N-mm)
+Mb = (P * 1000)*v
+//Calculate the torsional moment Mt (N-mm)
+Mt = (P * 1000)*h
+//Assume value of diameter to be 1mm d
+d = 1
+//Calculate the value of y (mm)
+y = d/2
+//Calculate the second moment of area I (mm4)
+I = (%pi/64)*(d^4)
+//Calculate the polar moment of area J (mm4)
+J = (%pi/32)*(d^4)
+//Calculate the bending shear stress B (N/mm2)
+B = (Mb * y)/I
+//Calculate the torsional shear stress T (N/mm2)
+T = (Mt * (d/2))/J
+//Refer Fig.4.40(c) for Mohr's Circle
+taumax = sqrt(((B/2)^2) + (T^2))
+//Calculate the true value of diameter d (mm)
+d = ((taumax/tau)^(1/3))
+//Print results
+printf('The value of diameter(d) = %f mm\n',d)
+