1 files changed, 39 insertions, 0 deletions
diff --git a/1898/CH13/EX13.7/Ex13_7.sce b/1898/CH13/EX13.7/Ex13_7.sce
new file mode 100755
index 000000000..2ee337dac
--- /dev/null
+++ b/1898/CH13/EX13.7/Ex13_7.sce
@@ -0,0 +1,39 @@
+clear all; clc;
+
+disp("Scilab Code Ex 13.7 : ")
+
+//Given:
+d = 30; //mm
+r = d/2;
+L = 600; //mm
+sigma_pl = 150;//MPa
+
+//Calculations:
+I = (%pi/4)*(r^4);
+A = %pi*r^2;
+r_gyr = sqrt(I/A);
+K = 1;
+sl_ratio = (K*L)/(r_gyr);
+flag1 = 0;
+
+//Assuming the critical stress is elastic:
+E = 150/0.001;
+sigma_cr1 = (%pi^2*E)/(sl_ratio^2); //Pcr = (%pi^2*EI)/(l^2)
+
+
+if(sigma_cr1 > sigma_pl)
+    Et = (270 - 150)/(0.002 - 0.001);
+    sigma_cr2 = (%pi^2*Et)/(sl_ratio^2); //Pcr = (%pi^2*EI)/(l^2)
+    
+   if(sigma_cr2>150 & sigma_cr2<270) 
+       Pcr = sigma_cr2*A;
+       Pcr = Pcr/1000; //in kN
+        printf('\n\nThe critical load when used as a pin supported column = %1.0fkN',Pcr);
+   
+   end
+    
+       
+end
+
+
+