diff options
Diffstat (limited to '3776/CH6/EX6.5/Ex6_5.sce')
-rw-r--r-- | 3776/CH6/EX6.5/Ex6_5.sce | 30 |
1 files changed, 30 insertions, 0 deletions
diff --git a/3776/CH6/EX6.5/Ex6_5.sce b/3776/CH6/EX6.5/Ex6_5.sce new file mode 100644 index 000000000..551491f53 --- /dev/null +++ b/3776/CH6/EX6.5/Ex6_5.sce @@ -0,0 +1,30 @@ +clear +//Given +//We will divide this into three parts +F = 8 //k - force applied +d = 16 //inch -distance +l_1 = 1 //in +l_2 = 3 //in +b_1 = 4 //in +b_2 = 1 //in +A_1 = l_1* b_1 //in2 - area of part_1 +y_1 = 0.5 //in com distance from ab +A_2 =l_2*b_2 //in2 - area of part_1 +y_2 = 2.5 //in com distance from ab +A_3 = l_2*b_2 //in2 - area of part_1 +y_3 = 2.5 //in com distance from ab + +y_net = (A_1*y_1 +A_2*y_2 + A_3*y_3)/(A_1+A_2+A_3) //in - The com of the whole system +c_max = (4-y_net) //in - The maximum distace from com to end +c_min = y_net //in - the minimum distance from com to end +I_1 = b_1*(l_1**3)/12 + A_1*((y_1-y_net)**2) //Parallel axis theorm +I_2 = b_2*(l_2**3)/12 + A_2*((y_2-y_net)**2) +I_3 = b_2*(l_2**3)/12 + A_2*((y_2-y_net)**2) +I_net = I_1 + I_2 + I_3 //in^4 - the total moment of inertia +M_c = F*d*c_max +stress_cmax = M_c/I_net //ksi - The maximum compressive stress + +M_t= F*d*c_min +stress_tmax = M_t/I_net //ksi - The maximum tensile stress +printf("\n The maximum tensile stress %0.3f ksi",stress_tmax ) +printf("\n The maximum compressive stress %0.1f ksi",stress_cmax) |