From 7f60ea012dd2524dae921a2a35adbf7ef21f2bb6 Mon Sep 17 00:00:00 2001
From: prashantsinalkar
Date: Tue, 10 Oct 2017 12:27:19 +0530
Subject: initial commit / add all books

---
 3776/CH9/EX9.4/Ex9_4.sce | 28 ++++++++++++++++++++++++++++
 1 file changed, 28 insertions(+)
 create mode 100644 3776/CH9/EX9.4/Ex9_4.sce

(limited to '3776/CH9/EX9.4/Ex9_4.sce')

diff --git a/3776/CH9/EX9.4/Ex9_4.sce b/3776/CH9/EX9.4/Ex9_4.sce
new file mode 100644
index 000000000..f645233f8
--- /dev/null
+++ b/3776/CH9/EX9.4/Ex9_4.sce
@@ -0,0 +1,28 @@
+clear
+//Given 
+//
+b = 40.0    //mm - The width of the beam crossection
+h = 300.0   //mm - The length of the beam crossection 
+V = 40.0    //KN - The shear stress in teh crossection
+M = 10.0    //KN-m - The bending moment on K----K crossection 
+c = h/2     //mm -The position at which maximum stress occurs on the crossection
+I = b*(h**3)/12                //mmm4 - the moment of inertia 
+//Caliculations 
+
+stress_max_1 = M*c*(10**6)/I   //The maximum stress occurs at the end
+stress_max_2 = -M*c*(10**6)/I  //The maximum stress occurs at the end
+y = 140                        //mm The point of interest, the distance of element from com
+n = y/(c)                      // The ratio of the distances from nuetral axis to the elements
+stress_L_1 = n*stress_max_1    //The normal stress on elements L--L
+stress_L_2 = -n*stress_max_1   //The normal stress on elements L--L
+x = 10       //mm The length of the element
+A = b*x      //mm3 The area of the element 
+y_1 = y+x/2  // the com of element from com of whole system
+stress_xy = V*A*y_1*(10**3)/(I*b) //MPa - The shear stress on the element 
+//stresses acting in plane 30 degrees 
+o = 60 //degrees - the plane angle
+stress_theta = stress_L_1/2 +  stress_L_1*(cos((%pi/180)*(o)))/2 - stress_xy*(sin((%pi/180)*(o))) //MPa by direct application of equations
+stress_shear =  -stress_L_1*(sin((%pi/180)*(o)))/2 - stress_xy*(cos((%pi/180)*(o)))               //MPa Shear stress
+  
+printf("\n a)The principle stresses are  %0.2f MPa %0.2f MPa",stress_max_1,stress_max_2)
+printf("\n b)The stresses on inclines plane  %0.2f MPa noraml,  %0.2f MPa shear ",stress_theta,stress_shear)
-- 
cgit