summaryrefslogtreecommitdiff
path: root/3776/CH6/EX6.10/Ex6_10.sce
diff options
context:
space:
mode:
authorprashantsinalkar2017-10-10 12:27:19 +0530
committerprashantsinalkar2017-10-10 12:27:19 +0530
commit7f60ea012dd2524dae921a2a35adbf7ef21f2bb6 (patch)
treedbb9e3ddb5fc829e7c5c7e6be99b2c4ba356132c /3776/CH6/EX6.10/Ex6_10.sce
parentb1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b (diff)
downloadScilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.tar.gz
Scilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.tar.bz2
Scilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.zip
initial commit / add all books
Diffstat (limited to '3776/CH6/EX6.10/Ex6_10.sce')
-rw-r--r--3776/CH6/EX6.10/Ex6_10.sce40
1 files changed, 40 insertions, 0 deletions
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)