diff options
Diffstat (limited to '3776/CH6/EX6.16/Ex6_16.sce')
-rw-r--r-- | 3776/CH6/EX6.16/Ex6_16.sce | 20 |
1 files changed, 20 insertions, 0 deletions
diff --git a/3776/CH6/EX6.16/Ex6_16.sce b/3776/CH6/EX6.16/Ex6_16.sce new file mode 100644 index 000000000..2f20f6648 --- /dev/null +++ b/3776/CH6/EX6.16/Ex6_16.sce @@ -0,0 +1,20 @@ +clear +// +M = 10 //KN.m - The moment applied +I_max = 23.95*(10**6) //mm4 - I_z The moment of inertia +I_min = 2.53*(10**6) //mm4 - I_y The moment of inertia +o = 14.34 // degress the principle axis rotated +//Coponents of M in Y,Z direction +M_z = M*(10**6)*cos((%pi/180)*(o)) +M_y = M*(10**6)*sin((%pi/180)*(o)) +//tanb = I_z /I_y *tan14.34 +b = atan((I_max*tan((%pi/180)*(o))/I_min )) +B = (180/%pi)*(b) +y_p = 122.9 // mm - principle axis Y cordinate +z_p = -26.95 //mm - principle axis z cordinate +stress_B = - M_z*y_p/I_max + M_y*z_p/I_min //MPa - Maximum tensile stress +y_f = -65.97 // mm - principle axis Y cordinate +z_f = 41.93 //mm - principle axis z cordinate +stress_f = - M_z*y_f/I_max + M_y*z_f/I_min //MPa - Maximum compressive stress +printf("\n The maximum tensile stress %0.2f MPa",stress_B) +printf("\n The maximum compressive stress %0.2f MPa",stress_f) |