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/CH6/EX6.10/Ex6_10.sce | 40 ++++++++++++++++++++++++++++++++++++++++
 1 file changed, 40 insertions(+)
 create mode 100644 3776/CH6/EX6.10/Ex6_10.sce

(limited to '3776/CH6/EX6.10/Ex6_10.sce')

diff --git a/3776/CH6/EX6.10/Ex6_10.sce b/3776/CH6/EX6.10/Ex6_10.sce
new file mode 100644
index 000000000..43f31cf7d
--- /dev/null
+++ b/3776/CH6/EX6.10/Ex6_10.sce
@@ -0,0 +1,40 @@
+clear
+//Given 
+l = 50.0    //mm - the length of the beam 
+b =  50.0   //mm - the width of the beam
+M  = 2083   //Nm
+A = l*b     //mm2 - The area 
+//straight beam 
+I = b*(l**3)/12.0 //mm4 - The moment of inertia of the beam
+c_1= l/2          // the distance where the stress is maximum 
+c_2 = -l/2        // the distance where the stress is maximum 
+s_1 = I/c_1
+s_2 = I/c_2
+stress_max_1 = M*(10**3)/s_1 //MPa - the maximum strss recorded in the  crossection
+stress_max_2 = M*(10**3)/s_2 //MPa - the maximum strss recorded in the  crossection 
+printf("\n The maximum stress upward in straight case is %0.3f MPa",stress_max_1)
+printf("\n The maximum stress downward in straight case is %0.3f MPa",stress_max_2)
+
+//curved beam 
+//
+r = 250.0                 //mm Radius of beam curved 
+r_0 = r - l/2             // inner radius 
+r_1 = r + l/2             // outer radius
+R = l/(log(r_1/r_0)) //mm 
+e = r - R 
+stressr_max_1 = M*(10**3)*(R-r_0)/(r_0*A*e)
+stressr_max_2 = M*(10**3)*(R-r_1)/(r_1*A*e)
+printf("\n The maximum stress upward in curved case is %0.3f MPa",stressr_max_1)
+printf("\n The maximum stress downward in curved case is %0.3f MPa",stressr_max_2)
+
+//curved beam _2 
+//
+r = 75.0      //mm Radius of beam curved 
+r_0 = r - l/2 // inner radius 
+r_1 = r + l/2 // outer radius
+R = l/(log(r_1/r_0)) //mm 
+e = r - R 
+stressr_max_1 = M*(10**3)*(R-r_0)/(r_0*A*e)
+stressr_max_2 = M*(10**3)*(R-r_1)/(r_1*A*e)
+printf("\n The maximum stress upward in curved case2 is %0.3f MPa",stressr_max_1)
+printf("\n The maximum stress downward in curved case2 is %0.3f MPa",stressr_max_2)
-- 
cgit