blob: 43f31cf7dbb825272309f8a18382e17330ef7954 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
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)
|
