diff options
Diffstat (limited to '3776/CH6/EX6.10/Ex6_10.sce')
| -rw-r--r-- | 3776/CH6/EX6.10/Ex6_10.sce | 40 |
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) |