1 files changed, 29 insertions, 0 deletions
diff --git a/764/CH4/EX4.19.b/solution4_19.sce b/764/CH4/EX4.19.b/solution4_19.sce
new file mode 100755
index 000000000..d7fa73adc
--- /dev/null
+++ b/764/CH4/EX4.19.b/solution4_19.sce
@@ -0,0 +1,29 @@
+
+//Obtain path of solution file
+path = get_absolute_file_path('solution4_19.sce')
+//Obtain path of data file
+datapath = path + filesep() + 'data4_19.sci'
+//Clear all
+clc
+//Execute the data file
+exec(datapath)
+//Calculate the permissible stresses for steel parts sigmat (N/mm2)
+sigmat = Syt/fs
+//Calculate the eccentricity e (mm)
+Rn = (((bi + bo)/2)*h)/((((bi * Ro) - (bo * Ri))/h)*log(Ro/Ri)-(bi - bo))
+R = Ri + ((h*(bi + (2 * bo)))/(3 * (bi + bo)))
+e = R - Rn
+hi = Rn - Ri
+A = (h * (bi + bo))/2
+//Assume the value load P to be 1N 
+P = 1
+//Calculate the bending moment Mb (N-mm)
+Mb = R * P
+//Calculate bending stress at the inner fibre B (N/mm2)
+B = (Mb * hi)/(A * e * Ri)
+//Calculate the direct tensile stress T (N/mm2)
+T = P/A
+//Calculate the actual load carrying capacity P (N)
+P = sigmat/(B + T)
+//Print results
+printf('\nThe load carrying capacity(P) = %f N\n',P)