1 files changed, 38 insertions, 0 deletions
diff --git a/1898/CH6/EX6.16/Ex6_16.sce b/1898/CH6/EX6.16/Ex6_16.sce
new file mode 100755
index 000000000..22a9b104b
--- /dev/null
+++ b/1898/CH6/EX6.16/Ex6_16.sce
@@ -0,0 +1,38 @@
+clear all; clc;
+
+disp("Scilab Code Ex 6.16 : ")
+
+//Given:
+t1 = 15/1000; //m
+t2 = 20/1000; //m
+l = 250/1000; //m
+b = 200/1000; //m
+P = 2.4; //kN
+l_a = 2; //m
+l_b = 1; //m
+
+//Internal Moment:
+y1 = b/2;
+y2 = t2/2;
+A = (2*t1*b)+(t2*l);
+y_bar = ((2*y1*t1*b)+(y2*t2*l))/A;
+
+M = (P*l_a)+(1*y_bar);
+
+//Section Property:
+I1 = (1/12)*(l*t2^3) + (l*t2*(y_bar - y2)^2);
+I2 = (1/12)*(t1*b^3) + (t1*b*(y1 - y_bar)^2);
+I =I1+ 2*I2;
+
+//Maximum Bending Stress:
+c = b - y_bar;
+sigma_max = (M*c)/(I*1000);
+
+//Display:
+    
+  printf('\n\nThe maximum bending stress at section a-a = %1.1f MPa',sigma_max);
+
+     
+//-----------------------------------------------------------------END--------------------------------------------------------------------------
+
+