From b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b Mon Sep 17 00:00:00 2001
From: priyanka
Date: Wed, 24 Jun 2015 15:03:17 +0530
Subject: initial commit / add all books

---
 1898/CH6/EX6.18/Ex6_18.sce | 46 ++++++++++++++++++++++++++++++++++++++++++++++
 1 file changed, 46 insertions(+)
 create mode 100755 1898/CH6/EX6.18/Ex6_18.sce

(limited to '1898/CH6/EX6.18')

diff --git a/1898/CH6/EX6.18/Ex6_18.sce b/1898/CH6/EX6.18/Ex6_18.sce
new file mode 100755
index 000000000..c5938e163
--- /dev/null
+++ b/1898/CH6/EX6.18/Ex6_18.sce
@@ -0,0 +1,46 @@
+clear all; clc;
+
+disp("Scilab Code Ex 6.18 : ")
+
+//Given:
+M = 12; //kNm
+l_bc = 0.2; //m
+l_be = 0.4; //m
+
+//Internal Moment Components:
+My = (-4/5)*M;
+Mz = (3/5)*M;
+
+Iy = (1/12)*(l_be*l_bc^3);
+Iz = (1/12)*(l_bc*l_be^3); 
+
+//Bending Stress:
+sigma_B = (-Mz*1000*(l_be/2))/Iz + (My*1000*(-l_bc/2))/Iy;
+sigma_B = sigma_B/10^6;
+sigma_C = (-Mz*1000*(l_be/2))/Iz + (My*1000*(l_bc/2))/Iy;
+sigma_C = sigma_C/10^6;
+sigma_D = (-Mz*1000*(-l_be/2))/Iz + (My*1000*(l_bc/2))/Iy;
+sigma_D = sigma_D/10^6;
+sigma_E = (-Mz*1000*(-l_be/2))/Iz + (My*1000*(-l_bc/2))/Iy;
+sigma_E = sigma_E/10^6;
+
+//Orientation of Nuetral Axis:
+z = (0.45)/(sigma_E + sigma_B);
+
+//theta = -atan(4/3);
+tanA = (Iz/Iy)*(-4/3);
+alpha = atan(tanA);
+alpha = alpha*(180/%pi);
+
+
+//Display:
+
+    
+ printf("\n\nThe normal stress at B = %1.2f MPa',sigma_B);
+ printf("\nThe normal stress at C = %1.2f MPa',sigma_C);
+ printf("\nThe normal stress at D = %1.2f MPa',sigma_D);
+ printf("\nThe normal stress at E = %1.2f MPa',sigma_E);
+ printf("\nThe orientation of the nuetral axis = %1.1f degrees',alpha);
+ 
+ //------------------------------------------------------------------------END---------------------------------------------------------------------------------------
+
-- 
cgit