diff options
Diffstat (limited to '3776/CH8/EX8.4/Ex8_4.sce')
| -rw-r--r-- | 3776/CH8/EX8.4/Ex8_4.sce | 19 |
1 files changed, 9 insertions, 10 deletions
diff --git a/3776/CH8/EX8.4/Ex8_4.sce b/3776/CH8/EX8.4/Ex8_4.sce index 75a4498ff..fb3b4c6ff 100644 --- a/3776/CH8/EX8.4/Ex8_4.sce +++ b/3776/CH8/EX8.4/Ex8_4.sce @@ -3,18 +3,17 @@ clear // S_x = 3.0 //MPa _ the noraml stress in x direction S_y = 1.0 //MPa _ the noraml stress in Y direction -c = (S_x + S_y)/2 //MPa - The centre of the mohr circle +c = (S_x + S_y)/2 //MPa - The centre of the mohr circle point_x = 1 //The x coordinate of a point on mohr circle point_y = 3 //The y coordinate of a point on mohr circle -//Caliculations +//calculations -Radius = ((point_x-c)**2 + point_y**2**0.5) // The radius of the mohr circle -//22.5 degrees line is drawn -o = 22.5 //degrees -a = 71.5 - 2*o //Degrees, from diagram -stress_n = c + Radius*sin((180/%pi)*(o)) //MPa The normal stress on the plane -stress_t = Radius*cos((180/%pi)*(o)) //MPa The tangential stress on the plane +Radius = ((point_x)**2 + point_y**2)**0.5 // The radius of the mohr circle +//22.5 degrees line is drawn +o = 22.5 //degrees +a = 71.57 - 2*o //Degrees, from diagram +stress_n = c + Radius*sin((180/%pi)*(o)) //MPa The normal stress on the plane +ang = sind((-a)) +stress_t = Radius*ang //MPa The tangential stress on the plane printf("\n The normal stress on the 22 1/2 plane %0.2f MPa",stress_n) printf("\n The tangential stress on the 22 1/2 plane %0.2f MPa",stress_t) -printf("\n answer varies due to rounding off errors") - |