1 files changed, 44 insertions, 0 deletions
diff --git a/764/CH4/EX4.6.b/solution4_6.sce b/764/CH4/EX4.6.b/solution4_6.sce
new file mode 100755
index 000000000..e25066af7
--- /dev/null
+++ b/764/CH4/EX4.6.b/solution4_6.sce
@@ -0,0 +1,44 @@
+
+//Function to print the designation
+function[] = desig(z)
+    printf('\nISLB%d is suitable for this application\n',z)
+endfunction
+ 
+
+//Obtain path of solution file
+path = get_absolute_file_path('solution4_6.sce')
+//Obtain path of data file
+datapath = path + filesep() + 'data4_6.sci'
+//Clear all
+clc
+//Execute the data file
+exec(datapath)
+//Load supported by each beam P (kN)
+P = (W/2)
+//Distance of centre of gravity of the load from the rigid support d (mm)
+d  = (l/3)
+//Calculate the bending moment Mb (N-mm)
+Mb = (P * 1000)*d
+//Calculation of ratio Ixx/y (mm3)
+ratio = Mb/sigmab
+//Selection of beam
+//There are 5 standard beams provided 
+for i = 1:1:5 
+    x = Ixx(i)/(h(i)/2)
+//Break the loop if value of x exceeds value of ratio
+    if (x > ratio)
+        break
+    end
+end
+//Print the suitable beam designation obtained
+if (i == 1)
+    desig(150)
+elseif (i == 2)
+    desig(175)
+elseif (i == 3)
+    desig(200)
+elseif (i == 4)
+    desig(250)
+else
+    desig(250)
+end